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A general framework for non-abelian symmetries is presented for matrix-product and tensor- 
network states in the presence of local orthonormal basis sets. The two crucial ingredients, the 
Clebsch-Gordan algebra for multiplet spaces as well as the Wigner-Eckart theorem for operators, 
are accounted for in a natural, well-organized, and computationally straightforward way. The unify- 
ing tensor-representation for quantum symmetry spaces, dubbed QSpace, is particularly suitable to 
deal with standard renormalization group algorithms such as the numerical renormalization group 
(NRG), the density matrix renormalization group (DMRG), or also more general tensor networks 
such as the multi-scale entanglement renormalization ansatz (MERA). In this paper, the focus is 
on the application of the non-abelian framework within the NRG. A detailed analysis is given for 
a fully screened spin-3/2 three-channel Anderson impurity model in the presence of conservation 
of total spin, particle-hole symmetry, and SU(3) channel symmetry. The same system is analyzed 
using several alternative symmetry scenarios. This includes the more traditional symmetry setting 
SU(2) spin ® SU(2)® 3 arge , the larger symmetry SU(2) spin ® U(l) c har g c ® SU(3) c hannci, and their much 
larger enveloping symplectic symmetry SU(2) sp i n ® Sp(6). These three symmetry settings are com- 



gether with simple self-contained numerical procedures to obtain Clebsch-Gordan coefficients and 
' | ' irreducible operators sets. The symmetries considered can consist of abelian symmetries together 

D with arbitrary non-abelian symmetries with compact, i.e. finite-dimensional, semi-simple Lie alge- 
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I. INTRODUCTION 

Numerical methods for strongly correlated quantum- 
many-body systems are confronted with exponentially 
large Hilbert spaces. With a limited number of ex- 
act analytical solution at hand and with perturbative 
treatments for low-energy or ground-state physics of- 
ten insufficient, a certain systematic treatment with re- 
spect the Hilbert space is required. Besides quantum 
Monte Carlo approaches, that explore quantum systems 
stochastically, 1 a systematic state space decimation is 
provided by renormalization group (RG) techniques such 
as the density matrix renormalization group (DMRG) 2 



or the numerical renormalization group (NRG), 3 both 
highly efficient for quasi-one-dimensional systems, and 
since non-perturbative, considered essentially exact. 

Both methods, DMRG as well as NRG, are based on 
the same algebraic structure of matrix product states 
(MPS). 4,0 Initially introduced for one-dimensional sys- 
tems with MPS owing its name to this case, a wide range 
of activity has emerged within recent years to generalize 
MPS to tensor-networks for two- or higher-dimensional 
systems. 6-9 While clearly appealing from the point of 
area laws for entanglement-entropy, 10-12 tensor network 
states (TNS) often share the same disadvantage as linear 
systems with periodic boundary conditions within the 
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DMRG, namely that state spaces become intrinsically 
non-orthogonal. Therefore also the unique association of 
symmetry labels with each index in a tensor is compro- 
mised. This, however, can be circumvented by introduc- 
ing an emerging extra-dimension, which is at the basis of 
the recently developed multi-scale entanglement renor- 
malization ansatz (MERA). 13,14 Nevertheless, the tradi- 
tional DMRG approach applied to 2D systems 15 with 
open or cylindrical boundary conditions yet with long- 
range interactions has continued to provide a highly com- 
petitive, extremely well-controlled, even though numeri- 
cally expensive approach. 

Within both, traditional DMRG as well as NRG, state 
spaces of entire blocks are built iteratively by adding and 
merging one site at a time. Clearly, the single index de- 
scribing an effective basis for the entire block or site can 
be chosen orthogonal. Moreover, the basis states can 
be labeled in terms of the symmetries of the underly- 
ing Hamiltonian. Operators written as matrix elements 
in this very same basis therefore also share the same 
well-defined partitioning in terms of symmetry sectors. 
By grouping symmetry state spaces together, the Hamil- 
tonian becomes block-diagonal, while general operators 
usually obey well-defined selection rules between symme- 
try sectors. Consequently, the sparsity of these operators 
due to symmetry can be efficiently and exactly included 
in the numerical description, such that usually only a few 
dense data blocks with non-zero matrix elements remain, 
given the symmetry constraints. While this well repre- 
sents the advantage of implementing generic abelian or 
point symmetries in a calculation, the presence of non- 
abelian symmetries offers yet another strong simplifica- 
tion: many of the non-zero matrix elements are actually 
not independent of each other, bearing in mind, for ex- 
ample, the Wigner-Eckart theorem. Therefore going be- 
yond abelian symmetries, non-abelian symmetries allow 
to significantly compress the non-zero blocks in terms of 
multiplet spaces, 16,17 while also reducing their number. 
With the Clebsch-Gordan coefficient spaces factorizing," 
they can be split off systematically in terms of a tensor- 
product and dealt with separately. 

MPS is optimal for one-dimensional systems. When 
exploring systems that are not strictly one dimensional 
but acquire width, such as ladders of several rungs in 
DMRG or multi-channel models in NRG, the price to 
be paid for orthonormal state spaces is that one must 
represent the system as a one-dimensional MPS never- 
theless. This introduces longer-range interactions to the 
mapped f D system, with the effect that the typically re- 
quired dimensions of the state spaces to be kept in a cal- 
culation, grow roughly exponentially with system width. 
The number of symmetries then that (i) are available and 
(ii) are also be exploited in practice, decides whether or 
not a calculation is feasible. Abelian symmetries such 
as particle (charge Q) or spin (S z ) conservation are usu- 
ally implemented in DMRG calculations. However, only 
very few groups have implemented non-abelian symme- 
tries, and these are also constrained to SU(2) symme- 



tries only, 16 due to its complexity in the actual imple- 
mentation. General treatment of non-abelian symmetries 
within the MERA, on the other hand, is currently under- 
development. 9,1 8 NRG, in contrast, had been set up in- 
cluding non-abelian SU(2) spin symmetry from its very 
beginning, 3 dictated by limited numerical resources. So 
far, however, only a very few isolated attempts includ- 
ing more complex non-abelian settings exist within the 
NRG, 19 while to our knowledge there exists no general 
realization yet of arbitrary non-abelian symmetries in ei- 
ther method. 

This paper focuses on the systematic description and 
implementation of non-abelian symmetries of a given 
Hamiltonian within the generalized MPS framework. 
This naturally also does include the description of abelian 
symmetries where necessary, as they can be trivially writ- 
ten in terms of Clebsch-Gordan coefficients. While the 
focus within non-abelian symmetries belongs to SU(A) 
and the symplectic group Sp(2n), the generalization to 
other non-abelian symmetries or also point groups is 
straightforward once their particular Clebsch-Gordan co- 
efficients are worked out. In contrast to the well-known 
SU(2) then, general non-abelian symmetries, such as 
SU(iV > 3), represent a significant increase in algorith- 
mic complexity, in that they can and routinely do exhibit 
inner and outer multiplicity. The latter, for example, 
implies that in the decomposition of the tensor-product 
of two irreducible representations (IREPs) into a direct 
sum of IREPs, the same IREP may occur multiple times. 
Nevertheless, this can be dealt with properly on the al- 
gorithmic level, as will be shown in detail in this paper. 

While the presented non-abelian framework for gen- 
eral tensors is straightforwardly applicable to traditional 
DMRG as well as NRG, the paper focuses on the appli- 
cation within the NRG. Detailed results are presented 
for a fully screened spin-3/2 Anderson impurity model 
with SU(3) channel-symmetry [i.e. see Hamiltonian in 
Eq. (25)]. This model has been suggested as the effective 
microscopic Kondo model for iron impurities in gold or 
silver 20 , historically the first system where Kondo physics 
was observed experimentally. 21,22 Being a true three- 
channel system, this cannot be trivially rotated into a 
simpler configuration of fewer relevant channels. The re- 
sult is an extremely challenging calculation within the 
NRG that requires non-abelian symmetries for fully con- 
verged numerical results for reasonable coarse-graining 
of the continuous bath. The non-abelian symmetries 
present in the model considered are (i) particle-hole sym- 
metry in each of the three channels, SU(2)^ , (ii) to- 
tal spin symmetry, SU(2) sp j n , and (iii) channel symme- 
try, SU(3) c hannci- The non-abelian particle-hole SU(2) 
symmetry, however, does not commute with the chan- 
nel SU(3) symmetry, while the plain abelian charge U(l) 
symmetry does commute. Overall, this suggests a larger 
enveloping symmetry, which turns out to be the symplec- 
tic symmetry Sp(6) [for an introduction, see App. A 10]. 
With this, the following symmetry scenarios are consid- 
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ered and compared in detail, 

SU(2) spin ®SU(2)f h 3 rgo , 

SU(2) spin <g> U(l) c hargo ® SU(3) ch anncl , and 

SU(2) spin (g Sp(6). 

While the first setting represents a more traditional setup 
based on multiple sets of plain SU(2) symmetries only, 
the second setting already includes the larger channel 
SU(3) symmetry. Both of these symmetries do not cap- 
ture the full symmetry of the model, which finally is 
achieved by using the enveloping Sp(6) symmetry. 

Due to the internal two-dimensional structure of the 
SU(3) symmetry based on the fact that SU(3) has two 
commuting generators, i.e. is of rank 2, its multiplets 
have significantly larger internal dimension, in practice, 
up to over a hundred. Therefore despite the reduction 
of the particle-hole symmetry to a plain abelian symme- 
try, the second setting with the SU(3) channel symmetry 
allows to outperform the more traditional setup based 
on SU(2) symmetries only. Similarly, with Sp(6) a rank- 
3 symmetry, multiplets then easily reach dimensions of 
several thousands there, which allows to reduce multiplet 
spaces significantly further still. A detailed analysis of 
this is provided in this paper, with a more general self- 
contained introduction to non-abelian symmetries con- 
sidered given in the appendix [cf. App. C3]. 

From an NRG point of view, 23 a few essential steps 
are required. These are (i) the evaluation of relevant op- 
erator matrix elements required to construct the Hamil- 
tonian, (ii) the generic setup of an iteration, adding one 
site to the so-called Wilson chain, and finally, for thermo- 
dynamical properties (iii) also the treatment of the full 
thermal density matrix. 23 All of these steps are simple in 
principle, yet come with the essential challenge to have a 
flexible transparent framework for the treatment of non- 
abelian symmetries in practice. In this paper, such a 
framework is presented in terms of generalized contrac- 
tions of tensors in the presence of symmetry spaces, in- 
troduced as QSpaces below. 

The paper is thus organized as follows. Section II 
describes the MPS implementation of non-abelian sym- 
metries in terms of QSpaces. Section III describes the 
implications for calculating correlation functions in the 
presence of irreducible operator sets. Section IV gives 
a short review of the NRG together with specialties re- 
lated to non-abelian symmetries, such as calculating re- 
duced density matrices. This section also introduces the 
model Hamiltonian of a fully symmetric 3-channel An- 
derson model. Section V then presents explicit NRG re- 
sults, followed by summary and outlook. Finally, also an 
extended Appendix has been added to the paper. The 
latter is intended to provide a more general pedagogical 
self-contained introduction to non-abelian symmetries as 
they occur in fermionic lattice models, together with their 
actual implementation in practice in terms of QSpaces. 



II. MPS IMPLEMENTATION OF 
NON-ABELIAN SYMMETRIES 

Consider some Hamiltonian H that is invariant under 
a set of ns symmetries, 

«s 

S = (g)S\ (1) 

A=l 

that is, [H,SqI — 0, where a identifies the generator S x 
for the simple (non-abelian) symmetry S x . To be spe- 
cific, for example, S = SU(2) spin ® SU(2) c har g e = S 1 ®^ 2 
with A € {1,2} would stand for the combination of spin 
and charge SU(2) symmetry, respectively. The tensor- 
product notation in Eq. (1) indicates that the symmetries 
act independently of each other, that is [S%, S x ,] = for 
A + A'. 

Given the symmetries as in Eq. (1), this allows to or- 
ganize the complete basis of eigenstates of H in terms 
of the symmetry eigenbasis. Every state then belongs to 
a well-defined irreducible multiplet q x for each symme- 
try iS A . The multiplet itself has an internal state space 
structure that is described by the additional quantum 
labels q z . For example, in the case of S x = SU(2), q x 
(q*) corresponds to the spin multiplet S (the S z label), 
respectively. 

Thus all states in a given vector space can be catego- 
rized using the hierarchical label structure 

\qn;q z ) (state-space label structure), (2) 

where 

(i) q = (q 1 , q 2 , . . . , q ns ), to be referred to as q-labels 
(quantum labels), references the irreducible rep- 
resentations (IREPs) for each symmetry S x , A = 
1, . . . , ng. All states in given Hilbert space with the 
same q-labels are blocked together, to be referred to 
as symmetry block q. 

(ii) Given a symmetry block q then, the multiplet index 
n^g) identifies a specific multiplet within this space. 
It is therefore a plain index associated with given 
symmetry space q. Together with the q-labels, this 
forms the multiplet level which is considered the 
topmost conceptual level. Using the composite no- 
tation (qn) to identify an arbitrary multiplet, the 
subscript q to the multiplet index ntg) is considered 
implicit and hence is dropped, for simplicity. 

(iii) Finally, the set of labels q z = (ql,q 2 , . . . , q z s ) , to 
be referred to as z-labels, resolves the internal struc- 
ture of each multiplet in q. That is, for each IREP 
q x , referring to the symmetry S x in q, q x labels 
its internal IREP space. As such, the z-labels are 
entirely defined by the symmetries considered. By 
construction, the eigenstates of the Hamiltonian H 
are fully degenerate in the z-labels. 



5 



Here the symmetry labels q and q z describe the combined 
record of labels derived from all symmetries considered. 
In practice, states can mostly be treated on the higher 
multiplet level, while the lower level in terms of the z- 
labels is split off and taken care of by Clebsch-Gordan 
algebra and the coefficient spaces derived from it. 

When non-abelian symmetries are broken, they are of- 
ten reduced to their abelian subalgebra. This can be eas- 
ily implemented, nevertheless, consistent with the pre- 
sented framework. In particular, in the abelian case, 
the non-abelian multiplet labels q are absent, while the 
abelian q z quantum numbers remain. Therefore the q z 
labels can be promoted to the status of q-labels, q := q z . 
As a consequence, the concept of the actual q z labels 
becomes irrelevant (therefore subsequently, the q z la- 
bel space may simply be set to zero, q z :— 0). The 
corresponding Clebsch Gordan coefficients are all triv- 
ial scalars, i.e. equal to 1. Yet these "Clebsch Gordan 
coefficients for abelian symmetries" do maintain an im- 
portant role, in that they take care of the proper addi- 
tion rules that come with abelian symmetries, resulting 

Given the MPS background of NRG or DMRG, states 
spaces are generated iteratively, in terms of a product- 
space of a given effective state space with a newly added 
local site. Operators, on the other hand, are typically 
represented in local state spaces, and starting from there, 
they can be written in terms of matrix elements in the 
effective global state spaces. With this in mind, the im- 
plementation of non-abelian symmetries within the MPS 
framework therefore is based on the following two ba- 
sic observations with respect to state space and operator 
representations, respectively. 

(1) State spaces: consider two distinct state spaces, 
\Qn\Q z ) and \ql\q z ) that, for example, represent 
a large effective state space and a small new local 
state space, respectively. Assuming that both state 
spaces all well-categorized in terms of IREPs, then 
their tensor-product space can also be decomposed 
into a direct sum of new combined IREPs \Qh; Q z ) 
using Clebsch-Gordan coefficients (CGCs), 

\Qn;Q*)= E Eft)Ii' c Sl 

Qn;Qz ql;qz 

x\Qn;Q z )\ql;q z ). (3) 

Note that the Clebsch Gordan coefficients given by 

q]q z = (QQz'iQQzlQQz) (i) f ull V define the inter- 
nal multiplet space as specified by the Lie algebra, 
and (ii) determine the splitting, i.e. which output 
multiplets Q occur for given multiplets Q and q. 
On the multiplet level, on the other hand, where 

(i^.)^ combines the multiplets Qn and ql into 
WW 

the multiplet Qh consistent with the splitting pro- 
vided by the CGCs, the coefficients (^lOBj may 

WW 

encode an arbitrary unitary transformation within 



the h output space for each Q. The r.h.s. of Eq. (3) 
demonstrates, that the CGC spaces clearly factor- 

ize from the multiplet space A^j- as a tensor prod- 

ww 

uct. 

(2) Operators: the matrix elements of a specific irre- 
ducible operator set (IROP) F q , i.e. an IROP that 
transforms according to multiplets q for given sym- 
metries [cf. App. Eq. (A3b), or also Sec. A 7] within 
some symmetry space \Qn;Q z ) can be written us- 
ing the Wigner-Eckart theorem as 

(Q'n'-Q' z \Fl\Qn;Q z ) = ' Cfely., (4) 

with C [ %j Q , again the Clebsch-Gordan coefficients 
as in Eq. (3). On the multiplet level, the reduced 

matrix elements (-Fqq/)J^/ = (QV||F«||Qn) refer 
to the single irreducible operator set labeled by 
q, which is indicated by the superscript [1]. The 
Wigner-Eckart theorem thus allows to compactify 
the operator matrix elements on the l.h.s. of Eq. (4) 
as the tensor-product of reduced matrix elements 
and CGCs, as shown on the r.h.s. of Eq. (4). 

Therefore in both cases above, i.e. in all tensor objects 
relevant for a numerical calculation, the CGC spaces fac- 
torize. This allows to strongly compress their size, and 
thus to drastically improve on overall numerical perfor- 
mance. Moreover, note that in both cases, Eq. (3) as 
well as Eq. (4) the underlying structure comprises ten- 
sors of rank-3 throughout. This rank-3 structure holds 
for both, the reduced multiplet space as well as the CGC 
spaces. Therefore, in either case, the final data structure 
of either state space decomposition as well as reduced 
operator sets is exactly the same. It is implemented, in 
practice, in terms of what will be referred to as QSpace 
for general tensors of arbitrary rank. 

A. General quantum space representation 
(QSpaces) 

The generic representation, used in practice to describe 
all symmetry related tensors B, is given by a listing of 
the following type, 



q-labels 


reduced space ||£?|j 


CGC spaces 


{Qh 




{C}i 


{Qh 


B 2 


{C} 2 



By notational convention, an actual operator B will be 
written with a hat, while its representation in terms of 
matrix elements in a specific basis will be written with- 
out the hat, hence the corresponding QSpace is referred 
to as QSpace B. Many explicit examples of QSpaces 
are introduced and discussed in detail in the appendix 
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[Sec. C]. As an up-front illustration, consider, for exam- 
ple, the general Hamiltonian of a single spinful fermionic 
site in the presence of SU(2) symmetry in the spin (S) and 
charge sector (C), which can be written as the QSpace 
[see Eq. (C8)] 





' (S;C) (S';C>) 


\\H\\ 


CGC spaces 


-1 


|;0 §;0 
0;| 0;i 


hi 

|,0 

h i 

0.5 


l( 2 ) 1. 
1. 1(2) 



With every non-zero block listed as an individual row, 
one can see that the only two reduced matrix elements 
|| free to choose without compromising the SU(2)®q = 
SU(2) sp i„ <g) SU(2) c hargo symmetry are the parameters 
(numbers) h±/2 t o and ft<o,i/2- By definition, the Hamil- 
tonian is a scalar operator, therefore it is the only op- 
erator within its IROP, hence can be written as plain 
rank-2 QSpace (the third dimension for this IROP would 
be a singleton dimension, hence can be dropped). Be- 
ing a scalar operator, the Hamiltonian is block diagonal, 
which is reflected in equal symmetry sectors (C; S) and 
(C';S') in each row for first and second dimension, re- 
spectively. Moreover, in given case, the corresponding 
Clebsch-Gordan coefficient (CGC) spaces also result in 
trivial identities, with l' 2 ) the two-dimensional identity. 
Note that the full set of CGC spaces in each row needs to 
be interpreted as appearing in a tensor product with the 
multiplet space, here the reduced matrix elements h 1 / 2: o 
or ^0,1/2 [e.g. see Eq. (8) below]. 

In general, the representation of a tensor B of arbi- 
trary rank-r in the QSpace in Eq. (5) [with Eq. (6) an 
example for a rank-2 QSpace], only lists the non-zero, i.e. 
relevant symmetry combinations. Having r tensor dimen- 
sions, each of its r indices refers to its specific state space 
\qn; q z )i = \{q)inu with i = l,...,r, and hence car- 

ries its own label structure as in Eq. (2). The q-labels 
(q)i = {q x }i already represent the combined set of rtg 
IREP labels from all symmetries S x for the state space 
at tensor dimension i. In general, by convention, the 
internal order of the q-labels (q)i w.r.t. A is fixed and 
follows the order of symmetries used in Eq. (1). 

For a certain row v of the QSpace listing in Eq. (5) 
then, the set of r q-labels arc grouped into 

{Qh={{q)i,---,{q)rh- (7a) 

The reduced matrix elements are stored in the dense 
rank-r tensor B v indexed by u, with % = 1, . . . ,r. This 
is a plain tensor, with the multiplet spaces possibly al- 
ready rotated by arbitrary unitary transformations and 
truncated. This is also reflected in the fact that the in- 
dices rii are plain indices, i.e. carry no further internal 
structure. Finally, for every one of the A = 1 . . . ng sym- 
metries included, the corresponding CGC space is stored 
in the sparse tensors „, each of which is also of rank 
r. These CGC spaces are grouped into {C}„ in the last 
column, 

{C}„ = {C 1;1 ,,...,C„ s; „}. (7b) 



As the q-labels {Q}v also define the z-labels, there is no 
explicit need to store the z-labels (q z )i,v The internal 
running indices in C\„, however, are uniquely associated 
with the z-labels. Note also the different index setting: in 
contrast toEq. (7a), which contains aset of r q-labels, i.e. 
one for every dimension of the rank-r tensor B, Eq. (7b) 
contains a set of ng rank-r CGC spaces, i.e. one for every 
symmetry. 

In addition to the QSpace listing in Eq. (5), also the 
type and order of symmetries considered is stored with 
a QSpace, cf. Eq. (1), even though this is usually the 
same throughout an entire calculation. Moreover, note 
that the row or record index v in Eq. (5) is purely for 
convenience without any specific meaning, as the order 
of records in a QSpace can be chosen arbitrarily. Nev- 
ertheless, it is required to refer to a specific entry in a 
QSpace. 

For a given record v in the QSpace in Eq. (5) then, the 
reduced space and the CGC spaces are to be interpreted 
as an overall tensor-product, 

"S 

B„® {C}„ = B„<8> (0C s;i/ ), (8) 

8=1 

while, of course, this is never explicitly done in practice. 
Yet Eq. (8) demonstrates the single most important mo- 
tivation to implement non-abelian symmetries in a nu- 
merical computation. By splitting off the CGC spaces 
in terms of a tensor product, block dimensions can be 
strongly reduced for larger calculations with several sym- 
metries present. For the models analyzed in this paper, 
for example, this was typically an average dimensional 
reduction from plain abelian symmetries by a factor of 
10 up to several hundreds. Considering that both NRG 
and DMRG scale like 0(D 3 ) with D the typical dimen- 
sion of data blocks, this is an enormous gain in efficiency 
The factorized CGC spaces, on the other hand, can be 
dealt with independently, as will be explained in detail 
later. Assuming that usually the dimensions of the re- 
duced states spaces B v still exceed by far the typical di- 
mensions encountered for the CGC spaces, the latter bear 
little numerical overhead. Only for larger-rank symme- 
tries, such as the symmetry Sp(6) discussed later, mul- 
tiplet dimensions can become large themselves such that 
one needs to pay more attention to an efficient treatment 
of their corresponding CGC spaces [see App. C3b]. 

For QSpaces where the CGC spaces in Eq. (5) exactly 
correspond to the standard Clebsch-Gordan coefficients 
for each symmetry, one may argue that actually similar 
to the z-labels, it is not explicitly necessary to store the 
CGC spaces altogether, since these are known. This is 
true, indeed, for these particular cases, and CGC spaces 
may simply be referenced then. Nevertheless, the explicit 
storage of the CGC spaces with a QSpace as in Eq. (5) has 
practical value. When combining QSpaces through con- 
tractions, i.e. sum over shared indices, for example, quite 
frequently intermediate objects can arise that do have 
rank different, in particular also larger than 3 [e.g. see 
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the intermediate objects indicated by the dashed boxes 
marked by X in Fig. 4]. These then elude a description 
in terms of standard rank-3 CGCs. In this case, the ac- 
tual CGC spaces for intermediate QSpace are important, 
and even though they do not necessarily resemble the 
interpretation of the original standard rank-3 spaces of 
standard Clebsch-Gordan algebra anymore, these spaces 
will be referred to as CGC spaces nevertheless, owing to 
their origin. 

For specific algorithms such as NRG and DMRG then, 
one typically ends up with simple scalar operators such 
as the Hamiltonian or a density matrix, apart from inter- 
mediate steps where complex CGC structures can arise. 
Therefore the full combined contractions on the CGC 
level [e.g. see Fig. 4] could be replaced by analytical ex- 
pressions or sum rules for Clebsch-Gordan coefficients. 
Then, again, the explicit storage of the CGS spaces with 
the QSpaces would not be strictly necessary and could 
be avoided altogether. Note, however, that these afore- 
mentioned contractions must be implemented within the 
code dependent on the specific contraction and relying on 
explicit analytical expressions which, depending on the 
symmetry, may or may not be known. For the QSpace 
as outlined in this paper, on the other hand, no spe- 
cial treatment is required for specific contractions, and 
no explicit knowledge of possibly symmetry dependent 
CGC sum rules is required. The QSpace approach solely 
relies on the correct construction of the standard CGC 
spaces to start with, with the subsequent sums over CGC 
spaces performed explicitly numerically and not analyt- 
ically through exactly the same contraction as on the 
reduced multiplet level, as discussed in more detail later. 

Finally, the explicit inclusion of the CGC spaces al- 
lows to build in strong consistency checks in the actual 
numerical implementation. Imagine that the Hamilto- 
nian is built by a sequence of complex contractions. The 
Hamiltonian eventually must be a scalar operator, i.e. it 
is block diagonal in the symmetries and the CGC spaces 
reduce to plain identities. This can simply be checked at 
the end of the calculation, which thus provides a strong 
check of whether the symmetries have been implemented 
correctly or not. At the stage of intermediate contraction, 
however, the CGC spaces guarantee the correct splitting 
and weight distribution between different emerging sym- 
metry sectors. 



B. A-tensors and Operators 

Consider the prototypical MPS scenario as in Eq. (3) 
that takes some previously constructed state space \i) = 
\Qn; Q z ) and adds a new local state space \a) = \ql; q z ), 
e.g. a new physical site. The state spaces are thus com- 
bined in a product-space described in terms of the IREPs 
\j) = \Qn; Qz}- Here the states i, a, and j are introduced 
as notational shorthand for better readability. The prod- 
uct space then is spanned by \io~) = |cr)|i). The order of 
states in the latter product emphasizes that state \o~) is 



typically added after and thus onto the existing state \i), 
which is of particular importance for fermionic systems. 
In general, the combined states 24 




are described in terms of linear superpositions of the 
product space \ia) given by the coefficients A^\ hence- 
forth called A-tensor (rank-3) or A ^-matrices (rank-2). 
Without truncation, A^' denotes a full unitary matrix 
U(ia),j where the round bracket indicates that the in- 
dices i and a have been fused, i.e. combined into an ef- 
fective single index. The presence of symmetry and the 
proper categorization of state spaces, however, imposes 
certain constraints on this unitary matrix, as pointed out 
already with Eq. (3). In particular, the fully determined 

CGC spaces factorize from the A-tensor, allowing 

an arbitrary rotation in the reduced multiplet space 

only. For the specific case then, that the reduced multi- 
plet spaces are identical to partitions of identity matrices 
with a clear one-to-one correspondence still of input and 
output multiplets, the corresponding A-tensor will be re- 
ferred to as the identity A-tensor [see Fig. 2 later; for 
explicit examples, see App. Eq. (C4) or Eq. (C6)]. An 
identity A-tensor therefore represents the full state space 
still without any state space truncation, and is unique up 
to permutations in the combined output space. Its ex- 
plicit construction is a convenient starting point, in prac- 
tice, when merging new local state spaces with existing 
effective state spaces. 

The entire construction of an A-tensor can be encoded 
compactly in terms of a rank-3 QSpace. Both coefficient 
spaces m Eq. (3), C^ 2 l as well as A^l, directly enter 

the QSpace description in Eq. (5). A schematic pictorial 
representation of an A-tensor is given in Fig. 1. There 
the states i (j) represent the open composite index to 
the left (right), respectively, while a refers to the open 
composite index at the bottom. 

As already argued with Eq. (4), an irreducible oper- 
ator shares exactly the same underlying CGC structure 
as an A-tensor. Thus also its representation in terms 
of a QSpace is completely analogous. Consider an IROP 
set F q = {-Fg }) which transforms according to IREP 
q. Here, the composite index a = (ql;q z ), for short, 
identifies the specific operators in the IROP set. As al- 
ready indicated by the superscript [1] in Eq. (4), its as- 
sociated multiplet index I has the trivial range 1 = 1, 
since, by definition, the IROP represents a single IREP 
on the operator level. With the states \i) = \Q'n';Q' z ) 
and \j) = \Qn; Q z ) now representing the same state space 
within which the operator acts, with usually many mul- 
tiplets and different symmetries, the operator represen- 
tation of the IROP F q in the states (i\ and \j) is evalu- 
ated using the Wigner-Eckart theorem in Eq. (4) . Similar 
to the A-tensor earlier, the resulting factorization of the 
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Figure 1. (Color online) Schematic depiction of a rank- 
3 QSpace as an example for a basic building block for an 
MPS or a tensor network, where lines (boxes) represent in- 
dices (data spaces), respectively. Every index is assumed to 
refer to a state space with similar physical background, hence 
refers to the same global symmetries as in Eq. (1), and has 
the generic composite structure \qn;q z ) as in Eq. (2), where 
q z specifies the states within the CGC spaces. The rank- 
3 QSpace depicted can be interpreted in two entirely differ- 
ent ways while sharing exactly the same underlying algebraic 
structure. These are (i) the state space decomposition into 
IREPs and (ii) operator representation for a given 1ROP in 
a given basis (see text). For the general interpretation of the 
QSpace depicted, consider for simplicity, a single row v in 
Eq. (5). The set {Q}„ defines the q-labels for all tensor di- 
mensions (here a total of three). With the q-labels fixed, the 
corresponding multiplet index n indexes the typically large 
reduced multiplet space, indicated by the thick black lines 
for each tensor dimension. The corresponding reduced rank-3 
multiplet space A v is depicted by the large gray box in the 
background. Moreover, with the q-labels fixed, this fixes the 
IREPs for every tensor dimension and every symmetry. The 
resulting sparse CGC spaces are indicated by the small boxes 
around the center, with one box for every symmetry, such as, 
for example, abelian particle conservation, non-abelian spin 
SU(2), non-abelian channel SU(3), or other. By construction, 
all CGC spaces share the same rank as the underlying QSpace. 
Therefore each CGC space also has three lines attached, one 
for every tensor dimension. In general, the CGC spaces re- 
fer to finite multiplet dimensions for non-abelian symmetries, 
while for simpler symmetries, such as abelian symmetries, the 
CGC spaces actually become trivial, i.e. scalars. These, nev- 
ertheless, are also interpreted as having the same rank as the 
QSpace using singleton dimensions throughout. 



CGC spaces Cq q, together with the remaining multi- 
plet space Fqq, of reduced matrix elements directly enter 
the QSpace description in Eq. (5). 

So even though an operator is usually considered a 
rank-2 object, the fact that an IROP consists of an oper- 
ator set indexed by a, adds a third index to the QSpace. 
In contrast to the state interpretation of a for the A- 
tensor above, however, here the "index" a has a different 
interpretation in that it points to a specific operator in 



the IROP set. By convention, the operator index a will 
always be listed as third tensor dimension in its QSpace 
representation. Given the three-dimensional representa- 
tion of a general IROP, therefore its entire construction 
mimics the construction of an yl-tensor in terms of a 
QSpace. As a consequence, Fig. 1 exactly also resem- 
bles the QSpace structure of an IROP. The states i (j) 
used for the calculation of the matrix element represent 
the open index to the left (right), respectively, while the 
operator index a refers to the open index at the bottom. 

Scalar operators, finally, such as the Hamiltonian of 
the system or density matrices, represent a special case, 
since there the IROP set contains just a single operator. 
Therefore the third index, i.e. the operator index, be- 
comes a singleton and hence can simply be dropped [e.g. 
see Eq. (6)]. Scalar operators therefore are represented 
by rank-2 QSpaces. They are block-diagonal in their sym- 
metries, and their CGC spaces are all equal to identity 
matrices, with an example already given in Eq. (6). 



C. Multiplicity 

For general non-abelian symmetries, frequently inner 
and outer multiplicity occur. 25 ' 26 Both are absent in 
SU(2), yet do occur on a regular basis for SU(iV > 3). In- 
ner multiplicity describes the situation where for a given 
IREP, several states may share exactly the same z-labels. 
Let ml denote the number of times a specific z-label oc- 
curs within IREP q. Then the presence of inner multi- 
plicity implies m\ > 1 for at least one z-label. Within 
such degenerate subspaces an arbitrary rotation is al- 
lowed in principle. For global consistency, therefore the 
CGC spaces must adopt a well-defined internal conven- 
tion on how to deal with inner multiplicity. This issue, 
however, is entirely contained within the CGC algebra, 
which is explored in more detail in the App. A [e.g. see 
discussion following Eq. (A19), and App. Bl]. On the 
level of a QSpace, it is of no further importance other- 
wise. Essentially the only implication of inner multiplic- 
ity is q z -> (q z ,a z ) with a z = 1, . . . , m% [cf. Eq. (A19)], 
where m\ depends on the multiplet q. With this minor 
adjustment, it is assumed throughout that the z-labels 
fully identify the internal multiplet space. Note that, in 
practice, the extra label a z is never included explicitly. 
What is important, however, is a consistent internal mul- 
tiplet ordering that respects multiplicity [see App. B 1]. 

Outer multiplicity, on the other hand, describes the 
situation where in the state space decomposition of a 
product-space of two IREPs, qi and q 2 , the same output 
IREP q may appear multiple times, the number of which 

is specified by M^ 1 '® 1 [cf. App. Eqs. (A35-A38) and 
discussion]. Therefore outer multiplicity primarily also 
enters at the level of Clebsch-Gordan coefficients, as it is 
based on pure symmetry considerations. In contrast to 
inner multiplicity, however, outer multiplicity also affects 
the reduced multiplet space, as will be elaborated upon 
in what follows. 
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In the absence of outer multiplicity [i.e. m] 91 ' 92 ' < 1 
for all gi, (?2, and q of the symmetry, an example being 
SU(iV < 2)], all rows in the QSpace in Eq. (5) must have 
unique {Q} u . If this is not the case, then the rows can be 
made unique by combining the rows with the same {Q}. 
Assume, for example, {Q} v = {Q}v' with v ^ v': clearly, 
the {<5}'s are already the same. Having the same sym- 
metry labels, this refers to the same set of IREPs, hence 
also the CGC spaces of these records must be identical, 
up to a possible global normalization factor which can 
be associated with the multiplet space, instead. Further- 
more, given {Q}y — {Q} U ' , the A v and A u > data blocks 
do live in exactly the same vector spaces for each indi- 
vidual tensor dimension! Therefore A„ and A v * can be 
simply added up [here multiple contributions with the 
same {Q} are considered additive, consistent with gen- 
eral conventions regarding sparse tensors; otherwise, say 
having given the same matrix element twice with differ- 
ent values, would immediately lead to contradictions]. 

In the presence of outer multiplicity, on the other hand, 
the uniqueness of the q-labels {Q}„ in the QSpace in 
Eq. (5) has to be relaxed. The reason for this is as follows. 
Since outer multiplicity derives from the Clebsch-Gordan 
algebra as in Eq. (A38), the CGC spaces 



(~<[<lz\ 



(QQz-,qqz\aQ,Q z ) (10) 
[different 



where M [ 9- q] 



acquire an additional label a = 1, . 

from the a z used with inner multiplicity 

indicates the outer multiplicity in Q, given the product 
space of the IREPs Q and q. In terms of a QSpace object, 
one may therefore be tempted to enlarge the CGC space 
from rank-3 to rank-4, with the dimension of the 4th in- 
dex being equal to mJ^' 9 ' . This strategy alone, however, 

does not capture the full picture since outer multiplic- 
ity also enlarges and thus effects the multiplet space A v 
of an A-tensor. By definition, outer multiplicity means 
that different multiplets with the same q can emerge. 
The only way they can be distinguished is through their 
Clebsch-Gordan coefficients. Therefore rather than en- 



larging the CGC space in a QSpace, M~ 



[QM 



records with 



the same {Q} u are allowed, instead. These records have 
CGC spaces of the same rank-3 dimensions, which, how- 
ever, are clearly distinguishable, as they are orthogonal 
to each other [cf. appendix Eq. (A39)]. The M [ 9' q] sets of 
Clebsch-Gordan coefficients arising from outer multiplic- 
ity are thus spread over M^' q ^ records within a QSpace 
object. 

The situation in the multiplet space for an identity A- 
tensor is depicted schematically in Fig. 2. In the absence 
of outer multiplicity, each symmetry combination (q, q')i 
can only contribute at most once to a given symmetry 
space q" and gets its space allocated, as depicted, for 
example, for {q 1 q')^ in Fig. 2(a), having only one non-zero 
block (shaded block) within the q" output multiplet. The 
symmetry combinations (q, q')\ and (q, q')%, on the other 
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Figure 2. Effect of outer multiplicity on multiplet space (A v ) 
in terms of an identity A-tensor- Panel (a) Schematic depic- 
tion of the state space decomposition of two input multiplet 
spaces with unique symmetry combinations (q, q') into com- 
bined multiplets q" (rows and columns, respectively). State 
spaces of the same symmetry are grouped into blocks sepa- 
rated by solid lines (horizontally and vertically). For simplic- 
ity, an identity A-tensor is depicted, for which the individ- 
ual sectors in q" can be uniquely associated with the (q, q') 
they originate from. Hence each column, separated by solid 
lines, has exactly one shaded block considered non-zero, with 
all-zero blocks shown in white. Here vertical thin lines indi- 
cate sub-blocks that originate from different (q,q'), yet are 
eventually combined in the same block as they belong to the 
same symmetry q" (separated by thick lines). Now, in the 
presence of outer multiplicity a specific (q,q') can contribute 
to the same q" several times, as depicted schematically by 
the spaces Mi and M2 for the rows (q, q')i and (q,q')2, re- 



spectively, both showing a multiplicity of M q . 



2. Panel 



(b) depicts the enlarged multiplet space for the output mul- 
tiplet q" of panel (a) in order to accommodate the additional 
multiplets arising from outer multiplicity. Being an identity 
A-tensor, the entire block shown in panel (b) represents an 
identity matrix (in contrast to an arbitrary A-tensor, which 
may have an arbitrary unitary matrix in its place). The ver- 
tical lineup of (q, q') sectors is arbitrary, making the identity 
A-tensor unique up to permutations. The identity matrix 
shown in the panel is sliced into horizontal blocks as indi- 
cated, each of which is associated with its own unique CGC 
space [not shown] as derived from the Lie algebra of the sym- 
metry under consideration. Each of these slices then directly 
enters as a reduced multiplet space A v in a separate row in 
the QSpace as in Eq. (5). 



hand, show outer multiplicity, in that they result twice 
in the same multiplet q", i.e. M [q ; q]l = M [q ; qh = 2. 

For simplicity, in the absence of truncation and with- 
out any further unitary rotation, the tensor-product on 
the multiplet level can be represented as an identity A- 
tensor with a clear one-to-one correspondence of input to 
output multiplets. This is depicted in Fig. 2(b) in terms 
of an identity matrix in the reduced multiplet space. The 
identity matrix in panel (b) then is sliced horizontally 
into blocks for each {q,q') that contributes to q". In the 
presence of outer multiplicity, the state space for q" needs 
to be enlarged to accommodate the additional multiplets. 
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The slicing (horizontal solid lines) then also proceeds for 
every output multiplet resulting from outer multiplicity, 

as indicated in panel (b). As a result, M^,f ' slices are 
associated with exactly the same Q = {q,q',q"}, distin- 
guishable only through their Clebsch-Gordan coefficients. 
These slices directly enter as A v in separate rows in a 
QSpace as in Eq. (5). 

In summary, outer multiplicity requires an adaptation 
of the multiplet space, which is naturally incorporated 
into a QSpace by allowing multiplet entries with the same 
{Q} v labels yet with clearly distinguishable CGC spaces. 
That is, specific records are also considered to refer to 
different state spaces if their CGC spaces are not exact 
copies (up to a global factor that can be incorporated 
into the multiplet data) but rather orthogonal to each 
other [see App. Eq. (A39)]. In practice, this is checked 
within a small numerical threshold (~ 10 -12 ) accounting 
for numerical double precision noise. The great advan- 
tage of this prescription is that then multiplicities fall 
completely in line with the rest of the MPS algorithm 
without any specific further treatment. 

Finally, it is important to notice that the same con- 
cept of relaxing the uniqueness of the {Q} v labels ac- 
tually also can become relevant for symmetries that do 
not have intrinsic outer multiplicity in its actual sense. 
Yet, in fact, through contractions intermediate objects 
can arise of rank larger than three [e.g. see the QSpaces 
indicated by the dashed boxes marked by X in Fig. 4], 
where records in a QSpace with the same {Q}„ labels can 
also have incompatible CGC spaces, in the sense that they 
are not the same up to overall factors. In this case, also 
the uniqueness of the {Q}„ must be relaxed temporar- 
ily. For simplicity, this will also be referred to as outer 
multiplicity. 

D. Contractions 

The contraction of QSpaces will be introduced in the 
following in terms of a simple example, namely the or- 
thonormalization condition on the combined state space 
in a tensor-product space. Putting symmetry labels aside 
for the sake of the argument, the A-tensor Ay = (icr\j) 
in Eq. (9) combines the state spaces \ia) into a combined 
(possibly truncated) orthonormal state space \j). This 
directly leads to the standard orthogonality relation for 
an A-tensor, 

£4?4? = %'> (ii) 

ia 

which is a simple example for the simultaneous contrac- 
tion of two tensors w.r.t. to two common indices, here 
i and a. By construction, it is completely analogous in 
structure to the orthogonality condition of CGCs as in 
App. Sec. A39. Including symmetries, the contraction in 
Eq. (11) is depicted in terms of QSpaces in Fig. 3. Over- 
all, indices are represented by lines, and lines connecting 



two blocks such as the indices i and a are summed over, 
i.e. contracted. In practice, contraction of QSpaces as 
defined in Eq. (5) happens at several levels, since state 
indices are always composite indices that refer to a sym- 
metry basis of the type \qn;q z ). This implies for a con- 
traction J2i=i' °f two QSpace objects with respect to 
some common state space i and i' , that (i) the q-labels 
qi and qy of the QSpaces as in Eq. (5) must be matched 
for the indices i and z', respectively. For a given spe- 
cific match of rows v and v' then, (ii) this is followed 
by the contraction of the corresponding reduced multi- 
plet spaces, and (iii) by exactly the same contraction of 
the CGC spaces, one for each symmetry. Note that the 
rank of a QSpace and its index order are always shared 
by the multiplet space and CGC spaces for consistency. 
Hence the overall contraction of the QSpaces is directly 
reflected in the elementary contraction of the plain nu- 
merical tensors A„ and {C}„. That is, the contraction 
pattern depicted schematically in Fig. 3, drawn in terms 
of boxes with connecting lines, is exactly the same on 
all levels of the contraction. By collecting the remaining 
non-contracted q-labels, this forms a new entry v" in the 
resulting QSpace, with the (tensor) index order of the re- 
sulting tensor dimensions being the same for all {Q}v>>, 
A v » , and {C}„» for consistency. 

Finally, the resulting QSpace is made unique in the 
{Q} u n labels as far as outer multiplicity permits. Records 
can only be combined, i.e. summed over, iff the CGC 
spaces for given records are all the same up to global 
factors which can be absorbed into the multiplet data, 
instead (see Sec. II C). Outer multiplicity plays no spe- 
cial role with contractions otherwise. Note that indepen- 
dent of whether or not outer multiplicity is present, when 
specifying a subset of tensor dimensions within {Q} v for 
contraction, the resulting QSpace will, in general, always 
have many contributions to the same {Q} U ". For com- 
parison, consider the completely analogous case of reg- 
ular square matrices of dimension D > 1: a matrix el- 
ement (M)ij is uniquely identified in the overall index 
while for example, the index i is not sufficient as 
it refers to an entire row of matrix elements. Moreover, 
when two matrices Mi and Mi are multiplied together, 

[B) Sj , = (MiM 2 ) jf = ^(MiMM^, , (12) 

i=i' 

the common index space (second index of Mi and first 
index of Mi) is summed over, i.e. contracted. Every 
match i = i' results in a contribution. In particular, 
for some given j and f , all D matches i = i' contribute 
and are summed up to the same output space (J,j'). In 
the case of QSpaces the situation is exactly analogous. All 
matches i = i' in the q-labels g, and q^ for the contracted 
index must be included. The only real consequence of 
outer multiplicity is that in the resulting QSpace B in 
Eq. (12) not necessarily all records with the same {Q} v 
labels can be merged by adding them together. In the 
specific case of the contraction in Eq. (11), however, the 
resulting QSpace is simply the identity, and as such a 



11 



A 




Figure 3. (Color online) Contraction of (i) an A-tensor or 
(ii) an irreducible operator into a scalar. All indices speci- 
fied are composite indices of the type \qn;q z ). An yl-tensor 
describes a (truncated) basis transformation of the product- 
space of the new local space \a) with an effective previously 
constructed basis \i), resulting in the combined state space 
\j'} = Yli'v' <T ') K')i w hh the corresponding bra-space 

(j\ = ^2 i<7 A^*(i\(a\ depicted in the lower part of the fig- 
ure. The result is the scalar identity operator, reflecting the 
orthonormality condition Eq. (11). An entirely different inter- 
pretation of the same contraction pattern can be given when 
the A-tensor is replaced by an IROP F a . The contraction 
then describes Eq. (13b) and yields a scalar operator, with its 
generic QSpace representation schematically depicted to the 
right. 

scalar operator with unique {Q}„. 



E. Scalar operators 

Given the definition of an A-tensor in Eq. (9), the con- 
traction of the two QSpaces A and A* in Fig. 3 leads 
to the identity operator = Y^j m the possi- 

bly truncated combined space C [cf. Eq. (11)]. Clearly, 
this also provides a strong check on the numerical imple- 
mentation of the symmetries. In particular, i( c ) repre- 
sents a (trivial) example of a scalar operator, that can 
be described as rank-2 QSpace. The CGC spaces are 
all identity matrices (up to overall factors that can be 
associated with the multiplet space), and therefore the 
lines, that usually connect to the CGC spaces within a 
QSpace, can be directly connected through from j to j' 
on the r.h.s. of Fig. 3, with the CGC spaces themselves 
no longer shown. In given case, due the orthonormality 
condition in Eq. (11), also the reduced multiplet space 
is given by identity matrices. This actually also would 
allow to connect through the thick black line on the r.h.s. 
of Fig. 3, and thus also to skip the large remaining block 
on the r.h.s. for the reduced multiplet space altogether. 

Figure 3, however, allows yet an entirely different in- 
terpretation. Remember that an irreducible operator set 
F q has a completely analogous structure and interpreta- 
tion in terms of its internal CGC spaces when compared 
to an A-tensor (cf. Fig. 1). Therefore it must hold that 



the scalar-product-like contraction, 

F^F-F^J2 P K P U ( 13a ) 

9* 

also results in a scalar operator (note that through the 
Wigner-Eckart theorem, by convention, the state space 
associated with the right index of the operator F q is com- 
bined with the multiplet space q; cf. App. A 7). With 
a = (ql;q z ) and the further sum through the opera- 
tor (matrix) multiplication, Eq. (13a) shares exactly the 
same contraction pattern as discussed in Fig. 3 in the 
context of the orthonormality of A-tensors earlier. Here 
the resulting scalar operator, however, can have arbi- 
trary positive hermitian matrices in its multiplet space 
still, represented by the large gray box on the r.h.s. of 
Fig. 3. The reduction of Eq. (13a) to a scalar operator 
is also intuitively clear, given that the Hamiltonian it- 
self is typically constructed in terms of scalar operators 
of exactly this type [see, for example, App. Eq. (A48) 
or Eq. (A60) given the Hamiltonian in Eq. (A46)]. The 
notation in Eq. (13a) emphasizes that in the scalar prod- 
uct the same irreducible operator set F q must be taken, 
considering that the IROP F q is different from the IROP 
(Ft)«. Nevertheless, since (F£ )t ~ (F^l^, up to pos- 
sible signs originating from the definition of the CGC 
algebra [e.g. compare the QSpaces in App. Tbls. C3a 
and C3b and accompanying discussion], these signs are 
irrelevant in the scalar contraction. Hence it follows that 
also 

f 2 ^. f =Y,( p iy tp i ( 13b ) 

is a scalar operator, yet different from Eq. (13a), as in- 
dicated by the tilde on F 2 . Similarly, note that if the 
A-tensor had been contracted on the right instead of the 
left index in Fig. 3, this also would have yielded a scalar 
operator, namely a reduced density matrix up to normal- 
ization (e.g. Fig. 5 below using p). = 1). 



F. Operator matrix elements 

The typical calculation of matrix elements of opera- 
tors for iterative methods such as NRG or DMRG is de- 
picted schematically in Fig. 4. While the complex many 
body states are generated iteratively and described by A- 
tensors [cf. Eq. (9)], an elementary irreducible operator 
set f q , on the other hand, usually operates locally within 
the state space \<r) of a specific site. Therefore, the oper- 
ator is described initially in terms of the matrix elements 
f^l, = (a\f q \a'). The link to the many body states is 
given through the yl-tensor that connects given site to a 

generated effective state space \j) = J2ia ^lj'l cr )l*)- 
The matrix elements of an IROP in the combined space 
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Figure 4. (Color online) Typical evaluation of matrix el- 
ements given an A-tensor. The nested dashed boxes 
indicate the sequential order of contractions prior to the fi- 
nal contraction. In panel (a), the local IROP set f acts 
within the state space \a) of a given site. Its local matrix 
elements, {a\f q \a'), are assumed to be known and described 
in terms of the local rank-3 QSpace /. The local IROP set 
is mapped into the larger effective space linked through the 
A-tensor, \j) = ^2 icr AZ \cr}\i). The overall result is the rank- 
3 QSpace F on the r.h.s., i.e. the desired matrix elements 
F\ q ), = {j\f q \j'}. Panel (b) depicts a typical scalar nearest- 
neighbor contribution to a Hamiltonian H = [f q ]J, ■ [f q ]k+i of 
two consecutive sites, say k and k + 1 using their respective 
A-tensors. This contraction already uses an effective descrip- 
tion of the local operator fl at site k in terms of the QSpace 
F\ obtained from the A-tensor Aik) a t site k as in panel (a) 
form the prior iteration. Using the A-tensor An,+i) of site 
k + 1, the overall contraction can be completed as indicated. 

\j) then become, 

ia.i'cr' 

** J crcr' 

= [E A ^(J2f^%.,, d4) 

o cr' 

It is exactly this procedure that is depicted in Fig. 4(a). 
The matrix elements are calculated in a two-stage pro- 
cess. The sum in the round brackets of Eq. (14) (con- 
traction of a') is carried out first, leading to the tempo- 
rary rank-4 tensor with open indices (i, j', <r, q) [box X in 
Fig. 4(a)]. This rank-4 tensor then is contracted simulta- 
neously in the indices i and a with the A* tensor, provid- 



Figure 5. (Color online) Backward update of density ma- 
trix pfe given in the effective basis \j) of a system up to and 
including site k (right index) by tracing out the local state 
space \ffk) (middle index) given the basis transformation Ak 
that introduced site k. The result is the reduced density ma- 
trix pk-i in the effective basis \i) of the system up to and 
including site k — 1. 



ing the final result shown on the r.h.s. of Fig. 4(a). Quite 
generally, for contractions including several blocks as in 
Fig. 4, these are always done sequentially, adding one 
block at a time. This is explicitly indicated in Fig. 4 by 
the (nested) dashed boxes, with the final contraction con- 
necting the remaining tensor to the outer-most dashed 
box. Every individual contraction then follows the multi- 
stage process over composite indices as described earlier 
in Sec. II D. 

The so obtained effective description of an op- 
erator f q acting on site fc using Ak can be used then to 
describe, for example, the typical scalar nearest-neighbor 
contribution [f 9 ] k ■ [f q ]k+i to the Hamiltonian including 
site k + 1. This operation is shown in Fig. 4(b). In par- 
ticular, one may use the identity A-tensor A 1 ^^ for site 
fc + 1, such that the resulting Hamiltonian is constructed 
in the full tensor-product space |cr)fc+i|«)fe of the system 
up to and including site fc + 1. Here \i)k describes the 
effective space up to and including site fc, whereas | cr) fc+i 
describes the new local state space of site fc + 1. This ex- 
actly corresponds the two-stage prescription used within 
the NRG (and similarly also for the DMRG) to build the 
Hamiltonian for the next iteration: (i) the tensor-product 
space including the newly added site must be mapped 
into proper symmetry spaces. This is taken care of by the 
construction of the identity A-tensor ^^i- (ii) The new 
Hamiltonian is built using this identity A-tensor through 
contractions as shown in Fig. 4(b) [note that while the 
presence of outer multiplicity in QSpace / is typically in- 
herited by QSpace F through the basis transformation as 
in Fig. 4(a), the internal contraction over the IROP set 
index q in Fig. 4(b) eventually leads to a scalar contri- 
bution to the Hamiltonian, as discussed with Eq. (13b)]. 
After diagonalization and state space truncation in the 
combined state space, the part of the resulting unitary 
matrix describing the kept states can be contracted onto 
^fe+i' yielding the actual final Ak+i- 
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G. Density matrix and backward update 

Consider the density matrix p k = Y,jj>(pk)jj'\j){?\ 
given in the basis |j)(k)i which is assumed to include all 
sites of a system up to and including site k. With the lo- 
cal state space of the last site k described by \au)> tracing 
out this last site from the density matrix p k corresponds 
to contracting the A^-tensor that connected site k to the 
system's MPS, 

Pk-i = tr(p fc ) 

Ok 

ij,i'j',<? 

= Y^(A^p k A^"%\i)(if\. (15) 
a' 

Equation (15) leads to a density matrix pk-i, which now 
is written in the many-body basis which includes 

all sites up to and including site k — 1. This backward up- 
date is a well-known operation within the NRG. 17,23 ' 27 ' 28 
Its graphical depiction is given in Fig. 5 [note that the 
sum over i and i' in Eq. (15) connects to state spaces 
that are not yet contracted; hence these correspond to 
open indices in Fig. 5] . 

The backward update of the density matrix in Eq. (15) 
preserves its properties as a density matrix and as a scalar 
operator. The former directly follows from the realization 
that the orthonormality condition Eq. (11) with the A- 
tensor in the last line of Eq. (15) is exactly equivalent to 
a complete positive map. Moreover, by tracing out part 
of a system such as a site that has been added through 
a tensor product space and that itself can be fully cat- 
egorized using given symmetries, this procedure cannot 
break symmetries by itself. This is to say, that the partial 
trace in Eq. (15) preserves the property of a scalar opera- 
tor. However, the trace over CGC spaces adds important 
weight factors to the reduced multiplet spaces, which are 
crucial, for example, to preserve the overall trace of the 
density matrix during back-propagation. While the con- 
traction in Fig. 5 can be easily performed, in practice, 
without the explicit knowledge of these weights, their 
determination is straightforward and instructive, never- 
theless, as will be shown in the following. 

The contraction in Fig. 5 clearly also holds for the CGC 
spaces of every symmetry individually. Therefore it is 
sufficient to focus on one specific symmetry. Let i con- 
tain several multiplets qi, and consider, for simplicity, the 
special case where the local state space a contains one 
specific multiplet q a only. In addition, also the reduced 
density matrix p k is chosen such that it only picks one 
very specific multiplet qj. Focusing on the Clebsch Gor- 

dan coefficients C q ^ q j z = (qiqiz] q<rq<rz\qjqjz) for chosen 
symmetry then, which properly combine the irreducible 
multiplets qi and q a into the multiplet qj , the contraction 



in Eq. (15) with respect to the fixed qj is given by 

E {<ii'qi'z;q<7q<jz\q 1 q 1 >z){qiqiz]q<jq<7z\q q z)* ■ S qjtqj , z 
= ^2(qi>qz>z;q<rqaz\(^2 q«<iaz) 

where the &q jz q-, z in the first line comes from the as- 
sumption that the initial p k is a scalar. The last iden- 
tity follows from the fact that also f>k~i shall be a scalar 
operator. Alternatively, the last equality can also be un- 
derstood as a general intrinsic completeness property of 
Clebsch-Gordan coefficients. Either way, the remaining 
factor fq t q. in the last line must be independent of the 
z-labels. The factor f q . q . then, in a sense, reflects the 
weight of how the IREP qi together with the traced over 
IREPs q a contributes to the final total qj. If, for exam- 
ple, for fixed qi and the known set of q a some final total 
qj cannot be reached, then it holds f qiqj — for this case. 

From the scalar property of Pk-i, Eq. (16) can be fur- 
ther constrained to some specific qi — . Also summing 
over q iz = qi> z then, the second line in Eq. (16) becomes 
equal to tr(^ g \QjQjx){QjQjz\') = d qj , i.e. the internal 
multiplet dimension of the IREP qj . Together with the 
last line in Eq. (16), it follows, 

/•* = t- ( 17 ) 

as demonstrated, for example, for SU(2) in [17]. Note 
that Eq. (17) holds in general for arbitrary symmetries, 
and also in the presence of outer multiplicity. This fol- 
lows by recalling that one of the main assumptions that 
entered Eq. (16) was to pick one specific multiplet qj. 
This single IREP, however, may equally well also have 
been any of the multiplets resulting from outer multi- 
plicity, say multiplet qj — > qj <a , which nevertheless again 
leads to Eq. (17). 

H. Generalized A-tensor for tensor networks 

The prototypical A-tensor as defined in Eq. (9) com- 
bines two physically distinct state spaces in terms of their 
tensor-product space. One may be interested, however, 
in the case where three or more state spaces need to 
be combined into the description of a single combined 
state space, while also respecting symmetries. This sit- 
uation, for example, occurs regularly in the context of 
tensor network states. 6-9 Let m be the number of states 
spaces to be combined. Then this requires the general- 
ized Clebsch-Gordan coefficients (qiqiz', ■ ■ ■ ', q m qm,z\qqz) ■ 
Once known, in principle, they can be combined com- 
pactly into a generalized A-tensor of rank m + 1. The 
question is, how to obtain this generalized A-tensor in a 
simple manner in practice. 

For this, the QSpace structure introduced in this paper 
proves very useful. In particular, a generalized A-tensor 
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Figure 6. (Color online) Generalized A-tensor that com- 
bines multiple state spaces, i.e. four effective state spaces 
. . . , \h) together with one local degree of freedom \a). 
Here it is assumed that all input state spaces describe proper 
orthonormal state spaces that act in different spaces, such 
that they can be combined into a simple product space. The 
index j, finally, represents the common global state space. 
In particular, it can be used to truncate the global Hilbert 
space to the state space of interest (within the DMRG, this 
may simply be the ground state, where the index j, being 
a singleton dimension, simply may be skipped then). While 
the general Clebsch-Gordan coefficients for the entire object 
may not be easily available (object to the right), the overall 
A-tensor can be built iteratively by adding one state space at 
a time (object to the left), starting, say, from Ai which links 
the two state spaces i\ and 12 into the combined state space j'2 
and hence allows to employ Clebsch-Gordan coefficients in the 
usual manner. The state space j'2 can then be combined with 
state spaces 13, and so forth. Contraction of the intermediate 
indices j'2, ■ • ■ ,J*4, finally, leads to the generalized A-tensor to 
the right. 



can be obtained based on the iterative pairwise addition 
of individual state spaces, which is a well-established pro- 
cedure at every step. The situation is depicted schemat- 
ically in Fig. 6. To be specific, Fig. 6 considers four 
effective state spaces \i a ) = \q a q a .z) with a — 1,...,4, 
together with a local state space |cr) = l^gs^), thus hav- 
ing m = 5. This specific setting may correspond, for 
example, to the situation in a tensor network state that 
describes a two-dimensional system which, from the point 
of view of a specific site with state space <t, has four ef- 
fective states spaces to the top, bottom, left, and right, 
respectively. Note, however, that here the state spaces 
\i a ) with a — 1, . . . , m are assumed to be physically dif- 
ferent, in principle, such that the tensor product space is 
a well-defined concept. 

Starting with state spaces \i\) and |i2 in Fig. 6, their 
state space can be combined in terms of and identity A- 
tensor A\ d in the usual fashion using standard Clebsch- 
Gordon coefficients. The resulting combined state space 
lj'2) then can be combined with state space ii in terms of 
another identity A-tensor Alf, thus yielding |j 3 ), again 
using standard Clebsch-Gordan algebra. The procedure 



is repeated, until at the last step, for example, the local 
state space a is added, resulting in the full combined state 
space \j) properly categorized in terms of symmetries. 
The iteratively generated m — 1 identity A-tensors A™, 
on the other hand, can be contracted into a single tensor 
of rank m + 1 by contracting the intermediate indices 
J2j ■ ■ • jim-i- This then results in the desired generalized 
A-tensor, shown at the r.h.s. of Fig. 6. Furthermore, in 
the context of DMRG or tensor network states, one is 
typically interested in a single state such as the ground 
state of the system. In this case, the full combined state 
space |j) is truncated to a single state, such that the 
extra index \ j) becomes a singleton and can be dropped, 
for simplicity. 

The construction of a generalized A-tensors, as out- 
lined above, emphasizes the versatility of the QSpaces 
introduced in this paper. In particular, by explicitly in- 
cluding the CGC spaces in Eq. (5), generalized Clebsch- 
Gordan coefficients can be easily obtained in terms of a 
generalized A-tensor, which itself is constructed through 
an iterative transparent procedure. 



III. CORRELATIONS FUNCTIONS 

Correlations functions are usually calculated with re- 
spect to operators whose transformation under given 
symmetries is known. This is specifically so, as these op- 
erators often naturally derive from the same fundamental 
building blocks, that also enter the Hamiltonian. In the 
presence of non-abelian symmetries, a single specific op- 
erator then that is not a scalar operator, is usually part 
of a larger irreducible operator set. In practice, thus also 
its correlation function is calculated w.r.t. the full IROP, 
for simplicity, as will be explained in the following. 

Consider, for example, the retarded Greens function 



cm = (4114). 



(18) 



that, in the time domain, creates a particle of preserved 
flavor cr, and destroys it some time later. Clearly, a par- 
ticle with the same flavor must be destroyed later, other- 
wise the Greens function is zero, i. e. the Greens function 
is diagonal with respect to symmetries. Now in the pres- 
ence of symmetries, it must be possible to write the oper- 
ators d\ as part of an irreducible operator set, e.g. some 
spinor (IROP) ^ q that transforms according to IREP q 
with internal dimension d q (in the case of plain abelian 
symmetries, it typically holds d q — 1, i.e. the operator d\ 
is the only member of the IROP). Thus the calculation 
of the very specific correlation function with respect to 
specific elements d a and d\ above can replaced by the 
Greens function 



(19) 



To be clear, if d q > 1, this includes the scalar product of 
the spinor components, and thus one is actually calculat- 
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ing the same Greens functions as in Eq. (18) d q times, 

d q 

^G [a) {Lu) = j q G i , q (Lu). (20) 

with G( CT ) independent of cr within its multiplet, as im- 
plied by the round brackets. This apparent overhead, 
however, only affects the CGC space, so this is negligi- 
ble numerical overhead, yet makes the calculation con- 
ceptually simple. Specifically, when calculating matrix 
elements and their contribution to the Greens function, 
eventually all indices can be fully contracted, so there is 
no need for a special treatment of a specific z-label that 
represents a peculiar d\. Moreover, given the discussion 
of scalar operators in Sec. II E earlier, one realizes that 
the scalar product •Ify 1 ■ %j) of the IROP i\) q yields a scalar 
operator. 

In the following, two explicit prototypical examples for 
correlation functions are given that are used explicitly for 
the numerical results presented in this paper. The first 
example is the spin-spin correlation function or magnetic 
susceptibility \d (w) defined at some site d. In the pres- 
ence of spin SU(2) symmetry, 

Xd( w ) = (Sx,d\\Sx,d}oj = (Sy,d\\S y ,d)u> — (S z ,d\\Sz,d}ui 

= l($4$*) u - (21) 

Clearly, the local operator S\ = Sd • is a scalar op- 
erator, with the corresponding spinor Su) given by [cf. 
App. Eq. (All)] 

•■(:!::) 

The second example is the spectral function for a sin- 
gle spinful channel in the presence of spin and particle- 
hole SU(2) symmetry. The spinor is given by [cf. App. 
Eq. (A59)], " 



\-&tJ 

In the evaluation of the correlation function, 

d q =i 

= (c t ||c|) w + (cj||cj.) w + (ejej)^ + (cjpf}w 

the signs, including s = ±1, of the individual compo- 
nents are irrelevant. Given the spin symmetry and the 



fact, that in the presence of particle-hole symmetry spec- 
tral functions are symmetric with respect to lo = 0, and 
in general G B>B t(u) = (B\\B^) U = G B ^ B (-u), it fol- 
lows that all four contributions above describe exactly 
the same function, indeed, and therefore 

G^cj) = {cjcl) u = i(ft\\j)) u . 



IV. THE NUMERICAL RENORMALIZATION 
GROUP 

The non-abelian setup described above is straightfor- 
wardly applicable to the NRG. 3,29 Before doing so in de- 
tail, here a brief reminder of the essentials of NRG is 
given, followed by the introduction of the model Hamil- 
tonian to be analyzed. By construction, the NRG deals 
with so-called quantum impurity models - an arbitrary 
small quantum system (the impurity) that is in contact 
with a macroscopic non-interacting typically fermionic 
bath. Each part is simple to solve exactly on its own. 
The combination of both, specifically in presence of in- 
teractions at the location of the impurity, however, gives 
rise to strongly-correlated quantum-many-body effects. 

The systematic approach introduced by Wilson 3 was a 
logarithmic discretization in energy space of the contin- 
uum of the bath (coarse graining) , followed by an exact 
mapping onto a semi-infinite so-called Wilson-chain, with 
the intact impurity space coupled only to the very first 
site of this chain. Given the half-bandwidth W ;= 1 
of the bath, the discretization parameter A > 1, typi- 
cally A > 1.7, defines the logarithmic discretization in 
terms of the intervals ±[A~ m , A~( m+1 )] with m > an 
integer, and energies taken relative to the Fermi energy 
£t = 0. Each of these intervals is then described by a sin- 
gle effective fermionic state, with its coupling and exact 
energy position chosen consistently w.r.t. the hybridiza- 
tion of the original continuum model. 30 ' 31 The resulting 
discretized model is then mapped onto the semi-infinite 
Wilson-chain (Lanczos tridiagonalization). 32 Hereby, the 
logarithmic discretization of the non-interacting bath 
translates to an effective tight-binding chain, with the 
hopping tk ~ Kr k l 2 between sites k and k + 1, decay- 
ing exponentially in the discretization parameter A. The 
latter then justifies the essential renormalization group 
ansatz of the NRG in terms of energy scale separation 
- large energies are considered first, with approximate 
eigenstates at large energies discarded and considered 
unimportant for the description of the still following 
lower energy scales. Thus each site of the Wilson chain 
corresponds to an energy shell with a characteristic en- 
ergy scale Wfc = f (A + l)A~ fc / 2 . Here the constant a of 
order 1. is chosen such, tk/uJk — > 1 for large fc. 28 

In practice, when considering the system up to site 
k, the Hamiltonian of the rest of the system is ignored, 
equivalent to assuming degeneracy in the state space of 
the remainder of the system. With the full Hamil- 
tonian H including the Wilson chain up to site fc, its 
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eigenstates \s) k , H k \s) k — E%\s)k, and with \e)k an arbi- 
trary state of the rest of the system following site k, then 
the essential spirit of NRG after coarse graining of the 
bath can be condensed in the following approximation, 33 

H\se) k ~E k s \se) k , (22) 

expressing energy scale separation, with \se)k = |s)fe ® 
|e)fe. The energies E* are usually taken relative to the 
ground state energy Eq of iteration k, and rescaled by the 
energy scale ui k . All of this will be referred to as rescaled 
energies, and has the advantage that independent of the 
Wilson shell k, energies are always of order 1. 

In this paper, the state space truncation at a given 
NRG iteration is energy-based, i.e. all states with E* < 
Ek are kept, typically with Ek — 5 ... 7 in rescaled 
energies. The number of kept states Nk thus changes 
dynamically. 28 ' 29 ' 31 

A. Full density matrix 

Within the NRG, 3 a complete many-body basis set 
can be formulated from the state space discarded at 
every iteration. 33 Initially introduced for explicit time- 
dependence of quantum quenches, they actually can also 
be used to improve on existing calculations for thermody- 
namical quantities and expectation values, 34 with a clean 
extension to arbitrary temperatures using the full den- 
sity matrix (FDM). 23 The density matrix p = e~ 0tl /Z 
with (3 = l/ksT, fee the Boltzmann constant and T the 
temperature, obviously commutes with the Hamiltonian 
and is a scalar operator in itself. Within the FDM-NRG 
approach, 23 the density matrix 

P=\ E e ^ f>E ^\l n ^1z;e) kk (qn;q z ;e\ 1 (23) 

k\qn;q z ;e 

can be constructed straightforwardly in terms of a 
QSpace for every Wilson shell k. Here s = (qn) stands 
for the multiplet label for a given shell k. Note that the 
symmetry of the states e is irrelevant here, as this space 
is fully traced over. Given the usual practice of NRG to 
rescale and shift energies at every iteration, all of this, of 
course, must be undone before entering Eq. (23) [given a 
general thermal density matrix, of course, all energies in 
Eq. (23) must be (i) at the same energy scale, i.e. non- 
rescaled, and (ii) specified with respect to a common en- 
ergy reference, e.g. the overall ground state energy of the 
Wilson chain]. 23 

By construction, all eigenenergies E* n are degenerate, 
i.e. do not depend on the z-labels. With the reduced den- 
sity matrix being a scalar operator, therefore the CGC 
spaces in the QSpaces describing Eq. (23) are all propor- 
tional to identity matrices, leading to the overall normal- 
ization 

Z=^^- fe e-K,, (24) 

k\qn 
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Figure 7. Schematic depiction of the fully screened Kondo- 
Anderson hybrid model [Eq. (25) with m = 3] in the NRG 
setup of a Wilson chain. Three d-levels with onsite Hund's in- 
teraction of strength Jh couple uniformly to their respective 
channel with hybridization T. The semi-infinite Wilson chain 
for each channel represents a tight-binding chain with expo- 
nentially decaying couplings, that interacts with the other 
channels through the impurity only. For a given NRG iter- 
ation, all terms in the Hamiltonian of the same energy scale 
must be included simultaneously, leading to an extended Wil- 
son site [dashed boxes] of three spinfull fermionic levels with 
a state space of 3 4 = 64 states each. 

where d q is the internal dimension of multiplet q, and 
jN-k reflects the degeneracy w.r.t. the rest of the Wilson 
chain of final length N, with d w the state space dimension 
of a Wilson site. 23 

B. Model: symmetric three-channel system 

The historically first physical system where Kondo 
physics was observed was that of Fe impurities in Au. 21,22 
The effective microscopic model for this material, how- 
ever, is far from trivial. It was argued only very recently 
in an extended study 20 that the physics of the five d- 
orbitals of substitutional Fe in Ag or Au is dominated 
by 3-fold degenerate triplet space i 2g , with the doublet 
space e g split-off by crystal fields and thus playing a mi- 
nor role. Together with the effective spin 3/2 of the iron 
impurity, this then results in an SU(3) symmetric fully 
screened 3-channel Kondo model. 

The actual model analyzed 20 is depicted schematically 
in Fig. 7. It consists of m = 3 spinful d-levels compris- 
ing the impurity, that are interacting through the Hund's 
coupling of strength Jjj. Each of these impurity levels is 
coupled to its own spinful bath channel with uniform hy- 
bridization r. This leads to the Kondo- Anderson hybrid 
Hamiltonian, 

m— 3 

H=H d+ j2J2 [vf( d -v + H - c -) + £ p4v] . 

1=1 pa 

(25a) 

where all energies will be given in context in units of the 
half-bandwidth W := 1. The impurity is described by 

H d = - J H S • S, (25b) 

with the impurity spin 

m m 
i i—1 a a' 
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given in terms of the Pauli matrices r a with a € {x, y, z}. 
Here d| CT [cj pcr ] creates a particle with spin a G {t, |} on 
d-level z at energy Ed — [in bath channel i at energy 
£p], respectively. For J# > T, an effective spin-3/2 forms 
at the impurity, leading to a symmetric fully-screened 
spin-3/2 system. The resulting Kondo temperature Tj< 
decays exponentially with Ju/T, with Xk quickly becom- 
ing the smallest energy scale in the system. In practice, 
choosing Jh — 2T/(m + j) leads to comparable Kondo 
temperatures Tk for different m. Compared to the stan- 
dard Kondo Hamiltonian with S • s coupling of the dot 
spin S with the lead spin s, the Hamiltonian in Eq. (25) 
in terms of T and Jh also allows for charge-fluctuations, 
while the model maintains particle-hole symmetry. 

In particular, the Anderson-like model in Eq. (25) 
has the advantage that the impurity self-energy S(cj) 
can be evaluated within the NRG in a simple fashion. 
From a more technical point of view, this allows the 
straightforward calculation of an improved spectral func- 
tion from the self-energy. 35 The impurity Greens function 
[cf. Eq. (20)] 

= G{ i<T) (u>)-inG'( i(7) (u;), (26) 

consisting of real and imaginary part, respectively, is con- 
structed within the NRG framework, as usual, from the 
spectral function A {ia) {u) = - AlmG (ier ) (ui) = G" ia) (w). 
Subsequently, the real part G'^s (ui) is obtained through 

Kramers-Kronig transform of A( ia }(uj). 29 The calculation 
of the additional correlation function F/ icr \ (oj) then, 

F {ia){u) = ([d(i<r),Hd]\\ d[ itr ))^ (27) 

obtained similarly from its spectral part F'u a \(w) = 
— —ImF(w), allows to evaluate the self-energy at 
the impurity 35 

S(itr),J H = tt^-, (28) 

Note that, the commutator of the IROP d{i a ) with the 
scalar Hamiltonian in Eq. (27) again leads to an IROP 
w.r.t. the same IREP q. Moreover, by symmetry, both 
G( i(T ) and Fu a ^ are independent of (icr), as indicated by 
the subscript bracket, and hence will be skipped alto- 
gether in the following, for simplicity. 

1. Kondo limit from numerical perspective 

While the procedure to obtain the self-energy is 
straightforward for an Anderson-like model, there is no 
simple way to do so for the plain Kondo-like model with 
S-s interaction. 35 However, from the NRG point of view, 
the transition from one to the other is straightforward. 
That is, knowing that the Kondo temperature Xk decays 



exponentially with Jh/T, both, Jh as well as T can be 
taken much larger than the bandwidth W := 1 of the 
model, while keeping their ratio constant, 

J H , r » 1, ^ = const. (29) 

This is a well-known procedure in the analytical 
Schrieffer- Wolff transformation for the Anderson model 
into a Kondo model. 36 But, of course, exactly the same 
strategy can also be pursued here within the NRG [see 
Fig. 8 later]. For the local density of states at the 
impurity this leads to a well-separated nearly discrete 
contribution to the spectral function at \lo\ 3> 1 far 
outside the bandwidth. For the spectral range within 
the bandwidth, the actual spectral function for the 
Kondo-model emerges. In particular, this procedure al- 
lows to fully eliminate the free-orbital (FO) regime with 
strong charge-fluctuations in the Anderson-like model 
right within the first truncation step. From a numeri- 
cal point of view, this is desirable as the FO regime is 
typically the most expensive one. For example, for the 
model discussed here using the symmetries below, using 
energy-based truncation indicates that about a factor of 
5 ... 10 more multiplets are required for the FO regime as 
compared to the local moment (LM) or strong coupling 
(SC) regime at later NRG iterations [cf. Fig. 9]. Never- 
theless, by maintaining an Anderson-like description, the 
impurity self-energy remains easily accessible numerically 
within the NRG, even though essentially the correlation 
functions for the Kondo model are calculated. 



V. NRG RESULTS 

The model in Eq. (25) is a true three-channel system, 
in that it is not possible to simply decouple a certain uni- 
tary superposition of bath channels. Furthermore, within 
an NRG iteration, a site from each channel must be in- 
cluded as they have the same coupling strength, i.e. en- 
ergy scale, as schematically depicted in Fig. 7 [dashed 
boxes] . 

The non-abelian symmetries present in the system are, 

• total spin symmetry: SU(2) sp j n , 

• particle-hole symmetry in each of the three chan- 
nels: SU(2)f h 3 rgc , and 

• channel symmetry: SU(3) cha nnci- 

The latter symmetry SU(3) c hannci, however, does not 
commute with particle-hole symmetry, while it does 
commute with the total charge U(l) c harge, *-e. the 
abelian subalgebra of particle-hole symmetry [cf. App. 
Eq. (A72), and subsequent discussion]. Having non- 
commutative symmetries, however, directly suggests a 
larger enveloping symmetry, which in the present case is 
the symplectic symmetry Sp(6) [i.e. Sp(2m) with m = 3, 
cf. App. A 10]. 
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This allows to consider the following symmetry set- 
tings, 

SU(2) spin ® SU(2)f h 3 argo ee SU(2)f c 4 , (30a) 
SU(2) spin ® U(l) ch argo ® SU(3) channc i, and (30b) 
SU(2) spin ® Sp(6). (30c) 

All of these symmetry settings have been implemented, in 
practice, and applied within the NRG framework, with 
results presented in the following. The first setting in 
(30a) represents a more traditional NRG scenario based 
on a set of plain SU(2) symmetries. The second set- 
ting (30b) includes SU(3) c h a nnci together with the simple 
abelian symmetry U(l) c harge for total charge, while the 
last setting (30c) represents the actual full symmetry of 
the model. 

Even though the second setting in (30b) actually in- 
cludes an abelian component in terms of charge, it nev- 
ertheless represents a stronger symmetry as compared to 
the first setting (30a). Since SU(3) c h a nnci is a rank-2 sym- 
metry with two commuting z-operators, i.e. generators 
of the Cartan subalgebra, it possesses a £wo-dimensional 
multiplet representation. This results in much larger 
multiplets for setting (30b), with the bare SU(3) multi- 
plet dimensions easily reaching up to 100 {e.g. d q < 125 
for the NRG run underlying Fig. 12, cf. App. C3). As 
a consequence, this allows, on average, smaller multiplet 
spaces and thus better numerical performance. 

The first two symmetry settings in (30) emphasize dif- 
ferent symmetry aspects, yet allow to break certain sym- 
metries which, nevertheless, are present in the model 
Hamiltonian in Eq. (25). The first symmetry setup (30a) 
strongly emphasizes particle-hole symmetry, while it does 
not use the symmetric coupling of the levels to their 
respective channels. The channel symmetry can thus 
be broken without reducing the symmetry setting (30a). 
The second symmetry setting (30b), on the other hand, 
emphasizes the channel symmetry (uniform T), while it 
allows to break the particle-hole symmetry. Hence, in 
principle, a uniform level-shift could be applied to the d- 
levels within this setting. Only the third symmetry (30c) 
captures the full symmetry of the model, as it combines 
channel symmetry with particle-hole symmetry into the 
enveloping symmetry Sp(6). This is a rank-3 symme- 
try with multiplet dimensions now easily reaching up to 
a several thousands (e.g. see Tbl. VIII for actual mul- 
tiplets generated in a full NRG run). A more detailed 
general discussion and comparison of all of above sym- 
metry setups in terms of their overall multiplet structure 
and representation of a site with three spinful levels (i.e. 
a Wilson site) is given in App. C 3. 



A. Spectral functions 

The spinor i\) q to be used for fermionic hopping term 
as well as for the calculation of spectral functions can 
be represented for all symmetry settings by IROPs with 



a well-defined multiplet label q [cf. Eq. (20)]. For the 
first symmetry setting in (30), SU(2)fg, the IROP for 
the calculation of spectral functions involves three 4- 
component spinors, ipf^ for short, one for every chan- 
nel i = l,...,(m = 3). The corresponding IROP la- 
bels are q 1 = (|, |, 0, 0), q 2 = (|, 0, ~,0), and q 3 = 
(|,0, 0, |), respectively. The number of components in 
the spinor derives from the two participating SU(2) mul- 
tiplets (S, C) = I for spin and one specific channel, thus 

having 2x2 = 4 operators in one specific IROP ijj\ 4 \ 
indeed. With this, the hopping in the Hamiltonian, for 

example, is given by hu,k+x = Ya=i ' $f+i,i' Notc 
that this excludes the hermitian conjugate part, as this is 
already fully incorporated through the particle-hole sym- 
metry [see App. Eq. (A73b)]. Furthermore, note that 
particle-hole symmetry gives rise to intrinsic even-odd 
alternations for the spinors along a chain [see App. A 9 b 
for a detailed discussion]. 

In contrast, the second symmetry setting in (30), 
SU(2) spin <g) U(l) chargc ® SU(3) ch annci, leads to a single 
6-componcnt spinor, tjj^ for short. Its IROP multi- 
plet label is given by q = (§,— §,0 1). This com- 
bines a 2-dimensional SU(2) spin multiplet S — | and 
an abelian 1-dimensional multiplet C z = — \ with the 
3-dimensional SU (3) channel multiplet T = (0 1), result- 
ing in the 2x1x3 = 6 operators in the multiplet. 
For comparison, here the hopping term in the Hamil- 
tonian in Eq. (25) is reduced to a total of two terms, 

h k , k+1 = (ViP • 4+a + H.c.) [see App. Eq. (A73a)]. 

Finally, for the third symmetry setting in (30), 
SU(2) sp i n ® Sp(6), again a single spinor ij} q is obtained, 
but now with 12 components, written as ^ 12 ^ for short. 
Its IROP label is given by q = (i, 1 0), which combines 
the 2-dimensional SU(2) sp i n multiplet S — \ with the 
6-dimensional Sp(6) multiplet (1 0), i.e. defining repre- 
sentation of Sp(6) . Overall, this again recovers the 2x6 = 
12 components of the spinor, indeed. For comparison, 
now the hopping term in the Hamiltonian in Eq. (25) is 

reduced to the single term hu,h+i = ($k 2 ^ ' The 
scalar contraction of the spinor fo 12 ^ with itself recovers 
the original 12 terms in the fermionic hopping structure 
between two sites in the Hamiltonian in Eq. (25). Since 
particle-hole symmetry is part of Sp(6), this again im- 
plies that (i) the hermitian conjugate is already taken 
care of in the hopping term, and (ii) that Sp(6) again 
gives rise to the same intrinsic even-odd alternations for 
the spinors along a chain, exactly analogous to what has 
already been encountered for standard particle-hole sym- 
metry. 

The correlation functions calculated for the model in 
Eq. (25) are presented in Fig. 8, with the model param- 
eters indicated at the bottom left of panel (a). Panel 

(a) shows the spectral data on a linear scale, while panel 

(b) shows the same data vs. log(|cj|) which therefore al- 
lows a logarithmic zoom into the low energy regime. The 
legend shown with panel (b) also applies to panel (a). 
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Figure 8. (Color online) SU(2) S pi n ®U(l) c h arge <g>SU(3) ch armci 
analysis of the symmetric three-channel Anderson model 
[Eq. (25) with model parameters specified in the lower left of 
panel (a)]; the same data as in panel (a) is shown vs. log(|cj|) 
in panel (b) to zoom into the Kondo peak at small frequencies 
with the legend for both panels shown in panel (b). The spec- 
tral data Anrg and the auxiliary F{^ RG are shown together 
with the derived self-energy S"(w) and the improved spectral 
function A- lmp (see text). A zoom around uj — is show in 
the right inset of panel (a), with the left (right) axis belonging 
to A(oj) [S"(w) and F"(uj)\, respectively. The spectral data 
for A(uS) and E" '(to) is symmetric around ui = and strictly 
positive, while F"(oj) is antisymmetric. In panel (b) therefore 
the ui < branch of F"{ui) has been plotted in dashed lines, 
same color otherwise. The left inset to panel (a) shows the 
spin-spin spectral data x"(^), with the resulting Tk = l/4xo 
indicated in panel (b) and the left insets of both panels by 
the vertical dashed line. 



The data in Fig. 8 was obtained using the symmetry set- 
ting in (30b) including SU(3) channel- Note that having 
chosen an energy-based NRG truncation with = 7, 
the spectral data for the other two symmetry settings 
is identical, hence not shown. While the calculation is 
somewhat more involved for the more traditional setup 
(30a), it becomes significantly more compact still when 
finally including Sp(6) as in (30c). Their individual nu- 



merical efficiency will be discussed with Fig. 9 below. 

In Fig. 8, the spectral function obtained from the NRG 
is plotted as j4nrg( w ) = G'^ rg (uj). The spectral data 
satisfies the Friedel sum-rule to an excellent approxima- 
tion, in that lim w _ s .o( 7 rr^NRG(w)) = 1 [see right inset to 
panel (a) for a zoom around to = 0]. The self-energy £(w) 
was obtained by calculating the additional correlation 
function F NRG (u) [Eq. (27), to be used in Eq. (28)]. The 
imaginary part S"(cj) = — ^Iml](u;), plotted in Fig. 8, 
clearly approaches zero in a smooth parabolic fashion at 
the Fermi energy, i.e. lim^^o S(w) = lim^^o — 
0, as expected for a system who's low-energy behavior 
corresponds to that of a Fermi liquid. This is seen more 
clearly still in the zoom around cj = in right inset of 
Fig. 8(a), with the self-energy data associated with the 
right axis. The self-energy E(w) / Jh sharply drops within 
M ^5 Tr from order 1. accurately down to about 10~ 4 
which is considered the NRG resolution limit. 

The improved spectral function Ai mp (uj) derived from 
the self-energy 35 is also shown in Fig. 8 [dashed red 
(black) line]. Within the Kondo regime, the result closely 
follows the original AnrgCw), as demonstrated in the 
zoom in the right inset of Fig. 8(a) or also in panel (b). 
As expected from the self-energy treatment, 35 the im- 
proved spectral function Ai mp (uj) allows clearly sharper 
resolution for structures at finite frequencies, specifically 
so for larger A. This can be observed, for example, for 
the hybridization side peaks in Fig. 8(b) at the energy of 
the Hund's coupling Jh- Having chosen Jh much larger 
than the bandwidth [with the bandwidth indicated by 
the vertical dotted line in panel (b)], these hybridization 
side peaks essentially correspond to very narrow, nearly 
discrete peaks that are much overbroadened through the 
standard log-Gauss broadening of the NRG. 23 ' 29 In prin- 
ciple, these side peaks could be narrowed significantly 
further by an adaptive broadening scheme. 3 ' For the pur- 
poses of this paper, however, this was irrelevant. 

The dynamically generated exponentially small Kondo 
temperature Tk for the system can be determined by 
taking the full-with-at-half-maximum (FWHM) of the 
Kondo peak in the spectral function. However, with 
NRG somewhat sensitive to broadening of the underly- 
ing discrete data 29 (see also supplementary material in 
Ref. 23), Tk is simply determined therefore through the 
static magnetic susceptibility xo ='■ 1/4Tk, 35 where xo 
is obtained from the impurity spin-spin correlation func- 
tion xM = (S( z ) jd \\S( z ) id ) u = x'(w)-wrx"(w) evaluated 
at uj = 0, with S( z ),d the total spin at the impurity [cf. 
Eq. (21)] . The resulting spin-spin spectral function ~x!'(oS) 
is shown in the left inset to Fig. 8(a), together with the 
resulting Tk = 4.4 • 10~ 7 (in units of bandwidth). As 
expected, x"( w ) shows a pronounced maximum around 
Tk- The value for Tk is also indicated by the vertical 
dashed line in panel (b). 

The NRG data presented in Fig. 8 clearly suggests con- 
verged data, even without necessarily having to resort to 
self-energy to get the low energy physics correct. 20 The 
convergence is also supported by the analysis of the dis- 
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carded weight 28 which, inspired by DMRG, analyzes the 
decay of the eigenspectrum of site-specific reduced den- 
sity matrices built from the ground state space a few it- 
erations later. For given NRG run, the discarded weight 
is estimated as e^ =5 % = 3 • 10~ n . This suggests good 
convergence, in agreement with Ref. 28. If, for example, 
an energy truncation of Ek = 5, . . . , 6 had been used, 
instead, NRG intrinsic parameter dependent deviations 
of up to ten percents can still be seen w.r.t. to the Friedel 
sum-rule or the agreement of ^4nrg( w ) with Ai mp (u>). 



B. Detailed comparison of symmetry settings 

An NRG specific technical comparison of the symme- 
try settings in Eq. (30) for the calculation in Fig. 8 is 
presented in Fig. 9. The underlying truncation had been 
energy-based in all calculations (Ek = 7), thus leading 
to a fair comparison in terms of accuracy. With this, 
the physical properties, and in particular the energy flow 
diagram 3 ' 35 in Fig. 9(a), show perfect agreement using 
either symmetry setting. Having sufficiently many states 
implies that for symmetries, that are not explicitly and 
thus exactly included in the QSpace setup, their unin- 
tended breaking due to numerical double precision noise 
does not play role. 

Figure 9(b) shows the number of kept multiplets for 
each iteration k. Having chosen Jjj and T much larger 
than the bandwidth [cf. Fig. 8(a)], the free-orbital regime 
is absent, with the transition from the local moment to 
the strong coupling regime given by the energy scale of 
Xk [vertical dashed line at k ~ 22]. As expected from 
physical grounds, also the local moment regime (k < 22) 
requires a larger state space (multiplet) dimension still 
for the same accuracy, i.e. the same £k ; as compared to 
the strong coupling regime (k > 22). 

With the state space truncation based on the energy 
cutoff Ex = 7, the actual Hilbert state space dimension, 
i.e. when including the internal CGC space dimensions, 
is exactly the same for all symmetry settings. In particu- 
lar, the maximum total Hilbert state space dimension per 
iteration that was diagonalized exactly for either symme- 
try setting was A* ot < 4,369,024 or AT* < 68,266 w.r.t. 
kept space only. These state spaces could be strongly re- 
duced to the effective and manageable multiplet dimen- 
sion as indicated at the top of Fig. 9(b), with Wilson 
shell specific multiplet dimensions plotted in the panel. 

Figure 9(c) analyzes the actual reduction in multiplet 
space due to presence of the CGS spaces in terms of the 
ratio of the actual Hilbert space dimension A t * ot relative 
to the total multiplet dimension A to t for each site along 
the Wilson chain. Depending on the symmetry setting, 
on average, the treatment of non-abelian symmetries al- 
lows to reduce the Hilbert space dimension by at least 
a factor of 16, 20, or 300 for the symmetries in (30), 
respectively. This demonstrates an enormous numerical 
gain, considering that the numerical cost of NRG roughly 
scales like C(A 3 ot ). Note that it is exactly through the di- 
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Figure 9. (Color online) Comparison of the efficiency of the 
symmetry settings as outlined in Eq. (30) for the calculation of 
the spectral data in Fig. 8 for the 3-channel model in Eq. (25). 
For a fair comparison, all calculations were performed using 
the same energy-based truncation with Ek = 7 for the same 
discretization A = 4 as in Fig. 8. The vertical dashed lines 
in all panels indicates the energy scale of Tk. Panel (a) com- 
pares the energy flow diagrams resulting for even iterations 
from the individual NRG runs, indicating perfect consistency 
for all symmetry settings. Panel (b) shows the number of kept 
multiplets for each iteration. For each symmetry setting, at 
the top of the panel the maximum dimension in the multi- 
plet space over the entire NRG run is specified for kept (Nk) 
and total (JV t ot), *-e. kept and discarded space, respectively. 
The inset shows the ratio qu = Nk/N^ of the multiplets 
that needed to be kept for the symmetry settings in Eq. (30a) 
and Eq. (30b) relative to the case when the full Sp(6) is in- 
cluded [Eq. (30c)]. Panel (c) shows the ratio N* ot /N to t of the 
actual Hilbert-space dimension (7V t * ot ) at a given iteration, 
which includes the internal multiplet dimensions, relative to 
the dimension of the multiplet space (-/V to t). The inset shows 
the ratio Xk that describes the increase in the number of mul- 
tiplets when adding a new site prior to truncation. 



mensional reduction to multiplet spaces, that above NRG 
calculations had been feasible in practice, and this within 
a few hours of runtime. In contrast, the plain abelian 
setting simply would not have been able to deal with the 
underlying Hilbert state space dimension using state of 
the art workstations [cf. App. Tbl. VI]. 

Within the kept space, the multiplet dimension of the 
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first two settings in (30a) and (30b) relative the setting 
including the Sp(6) are shown in the inset to Fig. 9(b) 
[(a) and (b), respectively]. This clearly demonstrates a 
further reduction by a factor of about 5 ... 8 when in- 
cluding the full Sp(6) symmetry. From the same inset, 
it is also clear that the symmetry setting in Eq. (30b) 
including SU(3) c h a nnci allows, on average, a 40% further 
reduction of the number of multiplets in the simulation 
as compared to the SU(2)fp setting. 

Furthermore, the inset to Fig. 9(c) shows the ratio 

Xk = ^tot +1 /Ny?' which indicates the increase in to- 
tal number of multiplets when adding a new site prior 
to truncation. While this factor shows a clear reduction 
from the actual dimension of the local Hilbert space of a 
Wilson site of 4 3 = 64 states, the ratio Xk is somewhat 
larger than what one may naively expect, considering 
that, depending on the symmetry, a Wilson site reduces 
to a total number of 4, 10 or 13 multiplets [see App. 
Tbl. V, Tbl. IV, and Tbl. Ill, respectively]. On the other 
hand, given non-abelian symmetries, the combination of 
two multiplets typically leads to clearly more than just 
one overall multiplet. In this sense, the major gain of 
using non-abelian symmetries is given by the state space 
reduction demonstrated in Fig. 9(b). For the first two 
symmetry settings in (30), the multiplet space increase 
by adding a new site in terms of a product space reduces 
the original abelian factor of 64 only modestly down to 
about 38. Only when using of the full Sp(6), this leads 
to a significant further reduction of the ratio Xk down 
to about 20, which thus becomes nearly comparable in 
numerical cost to a two-channel calculation with abelian 
symmetries, where a Wilson sites adds 4 2 = 16 states to 
the system. 

The SU(3) representations that are explicitly gen- 
erated in the calculation of Fig. 8 using SU(2) spin <E> 
U(l) 

charge 

<& SU(3) c hannci are listed in Fig. 12 [App. C3]. 
The largest Clebsch-Gordan space that is split off with 
respect to the SU(3) sector only is the (4,4) represen- 
tation with an internal multiplet dimension of 125. In 
other words, by explicitly accounting for SU(3) symme- 
tries, in the present case, a 125-fold degeneracy in the 
Hamiltonian had been reduced to a single multiplet, with 
the SU(3) symmetry space taken care of separately with 
minor computational overhead. Nevertheless, the eigen- 
states in the SU(3) setting still show significant degenera- 
cies. These can be entirely removed only by using the full 
Sp(6) symmetry, which allows to remove original degen- 
eracies in the Hamiltonian of several thousands. Note 
that on top of above symmetries, the spin SU(2) mul- 
tiplets present yet another independent multiplet space 
that enters as a tensor product, thus enlarging the overall 
symmetry space still further. 

In terms of overall runtime on a state-of-the-art 8-core 
workstation, this translated to about 6 hours of runtime 
for the SU(2)f^ symmetries, as compared to about 4.5 
hours of runtime when including SU(3) c h a nnci- Using the 
full symmetry as in (30c), on the other hand, took about 
24 hrs. While significantly more efficient in terms of stor- 



age requirements [cf. Tbl. VII] thus facilitating calcula- 
tions on standard workstations, the huge CGC spaces in 
the last setting must be dealt with carefully. As can be 
seen from Fig. 9(b), the total number of kept multiplets 
hardly reaches 400, while the Sp(6) multiplets are fully 
comparable in terms of dimensionality, with some mul- 
tiplets even much larger internally than the actual num- 
ber of multiplets considered [cf. App. Tbl. VIII]. While 
the sparse algebra had been optimized by ourselves to 
also make use of the parallel shared memory capacity [cf. 
App. C3b], in contrast, the full multiplet spaces had ac- 
cess to the highly optimized shared BLAS libraries. The 
latter benefitted the first two symmetry settings (30a) 
and (30b) in terms of overall runtime. However, there is 
clearly room for further improvement in dealing with the 
sparse algebra for larger rank symmetries as in (30c). 



VI. SUMMARY AND OUTLOOK 

A generic and transparent framework has been pre- 
sented for the implementation of non-abelian symme- 
tries in tensor-networks in terms of QSpaces. For this, 
it was assumed that all participating state spaces are 
strictly orthonormal and can be assigned proper well- 
defined symmetry labels. Therefore the presented frame- 
work is straightforwardly applicable to the traditional 
DMRG as well as to the NRG. The latter was demon- 
strated in detail in this paper for an SU(3) symmetric 
3-channel problem which, in the presence of particle-hole 
symmetry, can be enlarged into the symplectic symmetry 
Sp(6). By reducing the actual state space to the reduced 
multiplet space, while factorizing the Clebsch Gordan 
coefficient space, this allows an efficient description of 
all relevant tensors. While the explicit Clebsch Gordan 
algebra bears little overhead for combinations of lower 
rank symmetries, the average internal multiplet dimen- 
sions grow quickly with increasing rank r of a symmetry. 
In practice, one may roughly estimate that the typical 
internal multiplet dimension grows like O(10 r ), for ex- 
ample, having r = 0,1,2,3 for abelian, SU(2), SU(3), 
and Sp(6), respectively. Starting with r — 3, an effi- 
cient sparse scheme on all CGC spaces becomes crucial. 
For symmetries with rank larger than three, finally, it 
appears desirable to develop general strategies and sum 
rules for the contraction of extended complex networks 
of CGC spaces. 

A detailed self-contained general introduction to non- 
abelian symmetries is given in the appendix, together 
with many explicit examples that arise in practice 
(App. A and C). Several further highlights explained in 
detail in the appendix (App. B) are: (i) a straightforward 
numerical recipe for the general calculation of Clebsch 
Gordan coefficients based on explicit product space de- 
composition in the presence of multiplicity, (ii) a generic 
recipe for the determination of irreducible operator sets, 
and last but not least, (iii) also a general algorithm to 
get the framework for several symmetries initialized from 
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plain Fock space. The latter does not require any initial 
detailed knowledge of specific symmetry labels other than 
the general action of the underlying generators. These 
are known in second-quantized form and thus also easily 
defined in Fock space. 



VII. ACKNOWLEDGEMENTS 

I want to thank Jan von Delft and Arne Alex for fruit- 
ful discussions on SU(iV) representation theory, and also 
G. Zarand and I. Weymann for helpful discussions on the 
implementation of non-abelian symmetries in context of 
conventional NRG calculations. In addition, I want to 
thank the above people and also Markus Hani for their 
critical reading of the manuscript. This work received 
support from the DFG (SFB-631, De-730/3-2, De-730/4- 
2, SFB-TR12, WE4819/1-1). 



1 W. M. C. Foulkes, L. Mitas, R. J. Needs, and G. Rajagopal, 
Rev. Mod. Phys. 73, 33 (Jan 2001) 

2 S. R. White, Phys. Rev. Lett. 69, 2863 (Nov 1992) 

3 K. G. Wilson, Rev. Mod. Phys. 47, 773 (Oct 1975) 

4 S. Rommer and S. Ostlund, Phys. Rev. B 55, 2164 (1997) 

5 A. Weichselbaum, F. Verstraete, U. Schollwock, J. I. Cirac, 
and J. von Delft, Phys. Rev. B 80, 165117 (2009) 

6 A. W. Sandvik and G. Vidal, Phys. Rev. Lett. 99, 220602 
(2007) 

7 V. Murg, F. Verstraete, and J. I. Cirac, Phys. Rev. A 75, 
033605 (Mar 2007) 

8 J. I. Cirac and F. Verstraete, arXiv:0910.1130vl [cond- 
mat.str-el](2009) 

9 S. Singh, R. N. C. Pfeifer, and G. Vidal, Phys. Rev. A 82, 
050301 (Nov 2010) 

10 J. D. Bekenstein, Phys. Rev. D 7, 2333 (Apr 1973) 

11 M. M. Wolf, F. Verstraete, M. B. Hastings, and J. I. Cirac, 
Phys. Rev. Lett. 100, 070502 (Feb 2008) 

12 J. Eisert, M. Cramer, and M. B. Plenio, Rev. Mod. Phys. 
82, 277 (Feb 2010) 

13 G. Vidal, Phys. Rev. Lett. 99, 220405 (Nov 2007) 

14 G. Vidal, Phys. Rev. Lett. 101, 110501 (Sep 2008) 

15 E. M. Stoudenmire and S. R. White, arXiv:1105.1374v2 
[cond-mat .str-el] (201 1 ) 

16 I. P. McCulloch and M. Gulcsi, EPL (Europhysics Letters) 
57, 852 (2002) 

17 A. I. Toth, C. P. Moca, 6. Legeza, and G. Zarand, Phys. 
Rev. B 78, 245109 (2008) 

18 S. Singh, R. N. C. Pfeifer, and G. Vidal, arXiv:1008.4774vl 
[cond-mat. str-el] (2010) 

19 L. De Leo and M. Fabrizio, Phys. Rev. Lett. 94, 236401 
(Jun 2005) 

20 T. A. Costi, L. Bergqvist, A. Weichselbaum, J. von Delft, 
T. Micklitz, A. Rosch, P. Mavropoulos, P. H. Dederichs, 
F. Mallet, L. Saminadayar, and C. Bauerle, Phys. Rev. 
Lett. 102, 056802 (2009) 

21 W. J. de Haas, J. de Boer, and G. J. van den Berg, Physica 
1, 1115 (1934) 

22 J. Hondo, Progress of Theoretical Physics 32, 37 (1964) 

23 A. Weichselbaum and J. von Delft, Phys. Rev. Lett. 99, 
076402 (2007) 

24 U. Schollwock, Rev. Mod. Phys. 77, 259 (Apr 2005) 

25 J. P. Elliott and P. G. Dawber, Symmetry in physics, Vol. 
1+11 (Oxford University Press, New York, 1979) 

26 A. Alex, M. Kalus, A. Huckleberry, and J. von Delft, J. 
Math. Phys. 52, 023507 (2011), ISSN 00222488 

27 W. Hofstetter, Phys. Rev. Lett. 85, 1508 (Aug 2000) 



28 A. Weichselbaum, Phys. Rev. B 84, 125130 (Sep 2011), 
http://xxx.lanl.gov/abs/1106.0916 

29 R. Bulla, T. Costi, and T. Pruschke, Rev. Mod. Phys. 80, 
395 (2008) 

30 M. Yoshida, M. A. Whitaker, and L. N. Oliveira, Phys. 
Rev. B 41, 9403 (1990) 

31 R. Zitko and T. Pruschke, Phys. Rev. B 79, 085106 (2009) 

32 J. Demmel, J. Dongarra, A. Ruhe, and H. van der Vorst, 
Templates for the solution of algebraic eigenvalue problems: 
a practical guide (Society for Industrial and Applied Math- 
ematics, Philadelphia, PA, USA, 2000) ISBN 0-89871-471- 


33 F. B. Anders and A. Schiller, Phys. Rev. Lett. 95, 196801 
(2005) 

34 R. Peters, T. Pruschke, and F. B. Anders, Phys. Rev. B 
74, 245114 (2006) 

35 R. Bulla, A. C. Hewson, and T. Pruschke, J. Phys.: Con- 
dens. Matter 10, 8365 (1998) 

36 J. R. Schrieffer and P. A. Wolff, Phys. Rev. 149, 491 (Sep 
1966) 

37 A. Freyn and S. Florens, Phys. Rev. B 79, 121102 (Mar 
2009) 

38 C. Pope, Geometry and Group Theory (Online lecture 
notes, Texas A& M University, 2006) 

39 R. Gilmore, Lie groups, Lie algebras, and some of their 
applications (Dover Publications, 2006) 

40 J. J. Sakurai, Modern Quantum Mechanics, Revised ed. 
(Addison- Wesley Publishing Company, Reading, MA, 
1994) 

41 M. A. A. van Leeuwen, A. M. Cohen, B. Lisser, and 
C. Arjeh, LiE, A Package for Lie Group Computations 
(Universite de Poitiers, 1992) ISBN 90-74116-02-7, http: 
//www-math. univ-poitiers . f r/~maavl/LiE/ 



23 



Appendix A: Non-abelian symmetries 101 

The general more pedagogical introduction of non- 
abelian groups in this appendix emerges from a practical 
numerical background of treating quantum many-body 
phenomena. It does not claim to cover non-abelian sym- 
metries in every theoretical detail, yet requires certain 
elementary concepts which will be briefly reviewed. The 
main focus of this appendix then is on practical applica- 
tions in quantum lattice models. Specifically, this targets 
the numerical renormalization group (NRG), 3,29 density 
matrix renormalization group (DMRG) 2 ' 24 or more gen- 
erally tensor networks, 9,18 yet also exact diagonaliza- 
tion, which itself may be formulated in a matrix prod- 
uct state language. The appendix offers a general treat- 
ment of continuous non-abelian symmetries, with modi- 
fications towards abelian, point-groups, or discrete non- 
abelian symmetries straightforward. Overall, this ap- 
pendix should be self-contained, sufficient and hopefully 
helpful to deal with general abelian and non-abelian sym- 
metries in numerical simulations. 

The non-abelian symmetries of concern in this paper 
are continuous symmetries. An element G of the corre- 
sponding Lie group Q can be parameterized by a set of g 
continuous, independent, and real parameters a CT , 25 

g 

G(ai,...,a g ) =exp(i^2a a S c ), (Al) 

(T=l 

with g the dimension of the symmetry group. Infinites- 
imal operations with a a <C 1 then define the set of g 
generators {S a }, the number of which thus also reflects 
the dimension of the group (note that the identity oper- 
ator is a trivial operation which therefore is never part of 
the set of generators) . For unitary symmetries, as consid- 
ered throughout in this paper, the generators in Eq. (Al) 
are hermitian. Furthermore, when dealing with exponen- 
tially large yet /imie-dimensional quantum-many-body 
Hilbcrt spaces, the non-abelian symmetries also must 
have finite-dimensional Lie algebras. 

The commutator relations of the generators in 
Eq. (Al), 

[S a , Sfj] = ^ fa^i/Sv, (A2) 
v 

determine the tensor of the structure constants f^^u, 
which itself fully defines the underlying Lie algebra. The 
tensor / CTM „ is antisymmetric in that by construction 
fo-fiv — —ffj.au-, yet not necessarily fully antisymmetric 
also w.r.t. to the last index v [in principle, it can be 
made fully antisymmetric using the Cartan-Killing met- 
ric, while distinguishes between co- and contravariant in- 
dices in Eq. (A2); 38,39 for simplicity, however, this dis- 
tinction is not made in this paper]. All generators are 
assumed to be connected to each other through above 
commutator relations. That is, if a subgroup of gener- 
ators fully decouples in that it commutes with the rest 



of the generators, then this subgroup forms a symmetry 
of its own. In this sense the group of generators for a 
specific simple symmetry is irreducible. 

A set of matrices {R a }, that obeys exactly the same 
commutator relations as the generators (operators) {S&} 
in Eq. (A2), allows a one-to-one correspondence between 
the matrices {R&} and the generators of the symmetry. 
It is called a matrix representation of the Lie algebra. 
If the carrier space, i.e. the vector space within which 
the matrix representation is defined, is fully explored 
through repeated application of the individual matrices 
of the representation, then this is called an irreducible 
matrix representation, to be denoted as {I a } henceforth. 
It is unique up to an overall similarity transformation. 
Together with its carrier space it refers to an irreducible 
representation (IREP), specified by a unique label q. If, 
on the other hand, part of the carrier space of a matrix 
representation decouples, the representation is called re- 
ducible. This will be discussed in significantly more de- 
tail later in the context of state space decomposition in 
Sees. A 5 and A 6. 

Consider an irreducible matrix representation {!%} for 
IREP q of dimension d q . Its carrier space is spanned 
by the multiplet \q) = (l^^)}, where q z references the 
individual states within the multiplet q (consider, for ex- 
ample, spin multiplets, where \qq z ) = \S,S Z }). The states 
\qq z ) forms an irreducible space w.r.t. the action of the 
generators, in that for an arbitrary symmetry operation 
G as in Eq. (Al), 

G\qq z )=Y,G q q ^\qq l z ), (A3a) 

some linear superposition within the same multiplet 
space arises. The coefficients G q q , form a d q x d q di- 
mensional matrix, which represents the symmetry op- 
eration G within multiplet q, and is given by G q = 
exp (i J2 a 0>ala) f° r some arbitrary but fixed values a a . 

Similar to the multiplet space \qq z ) of dimension d q , 
an irreducible operator (IROP) set F q = {F q z } can be 
defined in a completely analogous manner. While it is 
not constrained to a specific carrier space, the IROP F q 
consists of a set of d q operators that are associated with 
multiplet q. As such, it can be written as a vector of op- 
erators, i.e. a generalized spinor. For a given symmetry 
operation G then, the IROP transforms analogously to 
Eq. (A3a), which for an operator implies 

On the level of infinitesimal operations, \a a \ <C 1, in con- 
trast to the plain action of generators on a ket-state as 
in Eq. (A3a), Eq. (A3b) shows that the transformation 
of an IROP directly translates to commutator relations 
[l.h.s. of Eq. (A3b)] with the generators of the symmetry, 
instead. 

The practical relevance of above general statements 
will be discussed in much detail in what follows, together 



24 



with many examples relevant in actual numerical calcu- 
lations. 



1. Simple example: rotational symmetry 

A simple and well-known example of a non-abelian 
symmetry is the rotational group 0(3) in real space in 
three dimensions. An arbitrary rotation can be written as 
G = e lS with S an arbitrary hermitian, yet fully complex 
three-dimensional matrix (hence no hats). The latter is 
required for G to be real. Consequently, this leaves three 
real parameters (a xi a y , a z ), with S = a x S x +a y S y +a z S z . 
The generators 25 
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(A4) 

represent infinitesimal rotations around the x, y, and z- 
axis, respectively. The 0(3) symmetry therefore has di- 
mension g — 3, and its Lie algebra is defined by, 



[S a 



(A5) 



with <7, /i, v € {x, y, z} and e a ^ v the Levi-Civita tensor, 
having switched to general operator notation [operator S 
(with hat) rather matrix S]. Being generators of 0(3), 
the matrix representation in Eq. (A4) already represents 
a 3-dimensional IREP. As it is the simplest non-trivial 
IREP for 0(3), it is also called its defining representa- 
tion. By combining state spaces that share this symme- 
try then, many other IREPs can be generated, including, 
for example, the (trivial) scalar representation of dimen- 
sion 1. 

With respect to continuous functions f(x,y,z) in 
three-dimensional space, the generators of infinitesimal 
rotations are given by the differential operator described 
by the angular momentum operator L = r x p with 
p ~ V r . By construction, its three components Lj also 
obey exactly the same Lie algebra as the generators in 
Eq. (A5). The same also holds for the spin algebra SU(2) 
in complex space, which describes the symmetry for spin- 
ful particles such as electrons if rotational spin symmetry 
is not broken, i. e. in the absence of an external magnetic 
field. Hence the rotational group 0(3) is isomorphic to 
the spin SU(2). In contrast to 0(3), however, the defining 
representation of SU(2) is two-dimensional [c/. Eq. (A6)], 
and hence also allows half-integer spin multiplets, which 
are entirely absent in 0(3). Having essentially twice as 
many multiplets in SU(2) as compared to 0(3), SU(2) is 
thus called a double cover or 2:1 cover of 0(3). 



2. SU(2) spin algebra 

In this paper, the setup and notation for non-abelian 
symmetries is generalized from SU(2). Therefore the 



symmetry SU(2) will be recapitulated in some more de- 
tail, introducing the semantics used for the general treat- 
ment of non-abelian symmetries. In this sense, the se- 
mantics used in this paper is somewhat more inclined 
towards the physics background, rather than strictly ad- 
hering to the mathematical language of Lie algebras. The 
latter, nevertheless, will be indicated in context. 

Similar to the 0(3) symmetry, an arbitrary unitary 
transformation in two-dimensional complex space is given 
by G = e lS with S an arbitrary two-dimensional hermi- 
tian matrix. This again has three independent real pa- 
rameters (a x , a y , a z ), such that S = a x S x + a y S y + a z S z . 
Here S a = with a € {x,y,z}, is given by the stan- 
dard Pauli spin matrices t ct , 
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For SU(2), this represents the smallest non-trivial ma- 
trix representation, therefore this also becomes its defin- 
ing representation. The commutator relations of the ma- 
trices T a are exactly the same as for 0(3) in Eq. (A5), 
since SU(2) also refers to the same rotational symmetry. 
Therefore, the generators for SU(2) will again also be de- 
noted by the operators {S a } with a € {x,y,z} in what 
follows. 

For a general irreducible representations of SU(2), e.g. 
a spin multiplet, the usual choice of basis is such that 
the z-component of the spin operator, S z , becomes diag- 
onal in its matrix representation S z , while the other two 
operators S x and S y remain non-diagonal (due to their 
non-commuting properties, only one spin component can 
be fully diagonalized, given the freedom of a similarity 
transformation for the whole representation). Using the 
notation \qq z ) = \S, S z ) for general spin multiplets, the 
multiplet label q (q-label) then can take the values q — 
0, |, 1, 1, 2, . . . with the internal multiplet label (z-label) 
spanning the 2q + 1 values q z G {—q, —q + 1,..., +q}- 
The raising and lowering operators (RLOs) are defined 
as 40 

S ± = S x ±iS y , (A7) 

such that S- = (S*+)^, with the commutator relations 

[&A]=±£± (A8a) 
[S+,S-] = [§+,&+] =2S Z . (A8b) 

For spin multiplets \qq z ) then, it holds 40 

S z \qq z ) = q z \qq z ) 

S±\qq z ) = ^q(q + l)-q z (q z ±l) \q,q z ±l). (A9) 

While the operator set {S x , S y , S z } generates the 
SU(2) symmetry group, this set itself does not yet repre- 
sent an irreducible operator (IROP), in that it docs not 
yet transform according to a specific symmetry multi- 
plet. For this, a specific linear superposition of the orig- 
inal operators as in Eq. (A8b) is required. In particular, 
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the transformation of an IROP set under given symme- 
try is completely analogous to the transformation of the 
symmetry eigenstates in Eq. (A9). As indicated with 
Eq. (A3b), the major difference is that the action of a 
generator applied onto a state is simply replaced by 
the commutator of the generator with an operator. For 
example, for an IROP F q given in terms of the set of op- 
erators {Fq x } which transform like the (state) multiplet 
q, it follows for consistency with Eq. (A9), 

[S z ,Fl\ = q z Fl (AlOa) 
[S± , f% } = ^q(q + l)-q z {q z ±\) ■ ±1 ( AlOb) 

This allows, for example, to complete the operator >S 2 
into an irreducible spin operator set as follows. Clearly, 
[S z , S z ] — 0-S z , which implies that the operator S z has z- 
label q z — 0, i.e. Sq = S z with q still unknown. Applying 
the RLOs yields the operators corresponding to q z = ±1, 

[S ± ,S%] =^q(q+l) + 0-S q ±1 . 

= [S ± ,S,]=tS± 

With the further application of RLOs yielding zero, i.e. 
[S + ,S + ] = [§-,§-] = 0, the operator space is thus ex- 
hausted. The irreducible spin operator set therefore has 
three members q z G { — 1, 0, +1}, and thus transforms like 
a spin multiplet q — max(q z ) = 1, 



S 1 = 01 
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Note that the signs and prefactors are crucial for consis- 
tency with the Wigner-Eckart theorem later. 

In above derivation, the z-operator in Eq. (AlOa) al- 
lowed to directly determine the z-label q z . The RLOs 
in Eq. (AlOb), on the other hand, served to explore the 
multiplet space, in that they generated the remaining op- 
erators Fy, with proper well-defined prefactors. In given 
case of spin SU(2), these factors are known [cf. r.h.s. 
of Eq. (AlOb)]. In situations, where they may not be 
known right away, they can nevertheless be determined 
in a straightforward manner. For simplicity, in the ab- 
sence of inner multiplicity for given multiplet, for canon- 
ical raising or lowering operator S± (see Sec. A 3 b) the 
combined application of S± followed by S± onto an op- 
erator of given multiplet q results in the same operator, 
i.e. 

[Sl[S ± ,Fl]] = alFl, 

from which the prefactor a± can be easily determined. 
The analogous situation for a state space multiplet \qq z ) 
is S±S ± \qq z ) = a±[qq z ), with a± > since S±S ± is 
a positive operator; in case of spin SU(2), this exactly 
reflects the prefactor on the r.h.s. of Eq. (AlOb), i.e. 



a± = q(q + 1) — q z (q z ± 1) > 0. Therefore if the appli- 
cation of S± results in a new operator component in the 
multiplet, i.e. a± > 0, then this operator is exactly given 
by 

F q q ,=^[S ±1 Fl\. (A12) 

This already contains the correct normalization and sign, 
with the the latter strictly determined by the outcome of 
the commutator. The z-label q z can be derived directly 
from the structure constants of the underlying Lie alge- 
bra, i.e. Eq. (AlOa). For a more general discussion on 
IROPs and their general decomposition also in the pres- 
ence of inner multiplicity for the IROP multiplet q, see 
Sec. A 7. 



3. Generators and symmetry labels 

Symmetries S within a quantum mechanical frame- 
work are described by a set of generators S a that leave 
the Hamiltonian H of the system invariant. Therefore it 
must hold for all generators of the symmetries considered 
that 

[S tT ,H]=0. (A13) 

Thus by definition, the Hamiltonian is a scalar opera- 
tor. The generators of independent symmetries S and 
S' commute trivially, by definition, as they operate in 
independent symmetry sectors. Therefore, for simplicity, 
a single specific non-abelian symmetry S is considered 
in the following, also referred to as simple non-abelian 
symmetry, a prototypical example being SU(iV). 

Therefore let S be a simple non-abelian symmetry. By 
construction then, its set of generators {S'o-} is fully con- 
nected via the structure constants in Eq. (A2), i.e. is irre- 
ducible but not necessarily an IROP yet [e.g. see previous 
discussion for SU(2) ]. With the symmetry reflected in 
the unitary transformation G — e leS ° with hermitian S a 
[cf. Eq. (Al)], it follows that for infinitesimal e <§C 1, the 
invariance of the Hamiltonian under this unitary trans- 
formation, i.e. UHU^ = H, is trivially equivalent to 
Eq. (A13). 

In order to ensure maximally independent generators, 
all operators in {S a } can be taken orthogonal with re- 
spect to each other and specifically also with respect to 
the identity matrix (which is always excluded from the 
set of generators {S a }). This requires a scalar or inner 
product for matrices, which is provided by 

{A,B)=tx{A^B), (A14) 

together with the resulting Frobenius norm ||A|| 2 = 
(A, A) = tiuVA). For the generators of the symmetry, 
thus one requires 

tY(SlS a ,) = a a 5 aa , (Al5a) 
tr^) = tr(l (t) 5 CT ) = 0, (Al5b) 
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The generators in Eq. (Al5b) are understood as finite- 
dimensional matrix representations of the operators Sa- 
in some specific carrier space, here the defining represen- 
tation. Moreover, the orthogonality w.r.t. to the identity 
in the last equation implies that all generators S a are 
traceless. 



a. Z-operators (Cartan subalgebra) 

For a given simple non-abelian symmetries, it is al- 
ways possible to identify a maximal set of mutually com- 
muting hcrmitian generators which form the so-called 
Cartan subalgebra of the symmetry's Lie algebra. These 
can be fully diagonalized simultaneously (together with 
the Hamiltonian) , and hence can be considered diago- 
nal. They shall be referred to as the z-operators [as they 
generalize the concept of the operator S z for SU(2)], 

[S z ,S z ,]=0 (z-operators). (A16) 

For a given Hamiltonian H then, this implies that every 
eigenstate H\n) = E n \n), in addition, can also be labeled 
with its respective set of symmetry labels \n) — > \qn;q z ), 
leading to 

H\qn;q z ) = E qn \qn-q z ). (A17) 

Here q identifies the multiplet, i.e. a set of states q z that 
are connected in an irreducible manner through all of 
the generators of the symmetry. While the index n origi- 
nally identified all states in given Hilbert space, it is now 
sufficient that it labels the multiplet within the space of 
multiplets that share the same q. The composite index 
(qn) then is referred to as multiplet index. Similarly, also 
the eigenenergies E qn in Eq. (A17) acquire symmetry la- 
bels. These, however, are independent of q z since, by 
construction, the states within a symmetry multiplet are 
degenerate in energy. More generally, with q z entirely de- 
termined by symmetry for a given multiplet q, they can 
easily be generated and thus omitted where convenient. 

Given a specific multiplet qn, the labels q z are equal 
to the eigenvalues of the z-operators, 

S z \qn\q z ) = q z \qn;q z ) (z-labels), (A18) 

which will be referred to as z-labels. If more than one 
z-operator is associated with given symmetry <S, say a 
total of r z-operators, where r thus defines to the rank 
of the symmetry, then the z-label structure associated 
with a multiplet also consist of a collective set of r z- 
labels (note that r needs to be differentiated here from 
the rank r of a tensor or QSpace as used in the main 
text). For example, the symmetry group SU(iV) has rank 
r = N — 1. Therefore the rank of SU(2) is 1, e.g. a single 
label q suffices to identify a state within an SU(2) spin 
multiplet. SU(3), on the other hand, already acquires a 
two-dimensional label structure for q z , and thus also for 

q- 



Note that the z-labels in Eq. (A18) for the states of a 
specific multiplet q may not necessarily be unique, in that 
the same q z may occur multiple times. Let m z describe 
how often a specific z-label occurs within given multiplet 
q. Then the presence of m z > 1 for at least one z-label is 
called inner multiplicity. It is then necessary to introduce 
an extra label a that uniquely identifies the state within 
this degeneracy, 

\Q n i Qz) \q n ] Qz a z)j (inner multiplicity) (A19) 

with a z € {1, . . . , m z }. While inner multiplicity is absent 
for SU(7V < 2), it occurs on a regular basis for SU(./V > 
3). The situation for outer multiplicity is analogous (see 
Sec. A 5). 

The label for the entire multiplet q (to be referred to 
collectively as q-labels) is in principle arbitrary, yet must 
be unique to identify the multiplet. Since for a continu- 
ous symmetry infinitely many IREPs exist, it is natural 
that the q-labels inherit the r-dimensional label structure 
of the z-labels. In particular, it is possible to construct 
a set of r scalar operators, called Casimir operators (see 
Sec. ??), that define a unique set of r constants for each 
multiplet. In practice, however, the q-labels are derived 
from q = m&x{q z }, i.e. by the z-labels corresponding to 
the maximum weight state (see Sec. A 3 c) which in prin- 
ciple can be related to the constants derived from the 
Casimir operators. 25 



b. Raising and lowering operators (roots) 

While for an arbitrary unitary element G of the sym- 
metry hermitian {S a } are required, on the level of gener- 
ators, in principle, arbitrary linearly-independent linear 
superpositions within the space of generators S a can be 
taken. Using such a reorganized set of generators, in- 
stead, this still preserves Eq. (A13), yet alters the struc- 
ture constants j ailv for given symmetry S. This freedom 
is used in the following to define canonical raising and 
lowering operators, which are non-hermitian, in general. 

Consider the action of a generator S a onto a symmetry 
eigenstate \qn;q z ). The z-operators are special, in that 
they are diagonal and hence return the same state, yet 
weighted by the eigenvalue q z . The remaining genera- 
tors, however, are non-diagonal, hence change the state 
and thus explore the multiplet space. In general, these 
generators can be reorganized such that all of them rep- 
resent proper raising or lowering operators (RLOs), with 
the canonical commutator relations, 

["S*,z > 5*0-] = f Z a<jS a = f Z aS a , (A20) 

with no summation over a. The action of these canonical 
RLOs in z-label space, in the literature also referred to as 
root space, then defines the canonical form. By definition, 
the canonical RLOs {S±} of a specific Lie algebra are ex- 
pected to have the property that their application onto a 
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symmetry eigenstate in the multiplet with well-defined z- 
labels will generate another eigenstate of the z-operators, 
yet with raised or lowered, i.e. well-defined different z- 
labels. This is exactly what is expressed through the 
commutator relations in Eq. (A20). In particular, the 
structure constants take the simple form, where a non- 
zero contribution can only arise if the last two indices 
in fzaa' are identical, hence the shortcut notation f za 
in the last term in Eq. (A20). By construction, f zrJ is 
fully antisymmetric. Note that Eq. (A20) also can be in- 
terpreted as an eigenvalue equation for the generators of 
the group. Since the z-operators S z are symmetric, the 
resulting eigenvalue problem is always well-defined with 
real eigenvalues f za . 

As a specific example, Eq. (A20) was already encoun- 
tered for SU(2) in Eq. (A8a). Here it states more gen- 
erally that the commutator of an arbitrary generator S a 
with a z-operator yields the very same operator S a up to 
the scalar structure factor f Z(T . This factor can be zero, 
e.g. when S a refers to another z-operator as in Eq. (A16), 
therefore / zz < = 0. For every z-operator, however, there 
must exist at least one RLO S a with f za =^ 0, since oth- 
erwise the group of generators would be reducible. 

With Eq. ( A20) , the application of a generator S a onto 
a symmetry eigenstate \qn;q z ) yields 

S z ■ S a \qn; q z ) = [S z ,S a ] \gn; q z ) + S a ■ S z \qn;q z )^ 

— J za •Ja 

= (qz + fza)S a \qn-q z ). (A21) 

If S a is an RLO with f za ^ 0, the state S a \qn;q z ) is 
again a symmetry eigenstate, yet with a uniform shift in 
the z-labels, 

q z , = q z + f za . (A22) 

Therefore the action of an RLO S a in root space is 
generic, i.e. independent of the specific multiplet q or 
the state q z under consideration. Nevertheless, the RLO 
may annihilate the state, i.e. S a \qn;q z ) — 0, which is 
essential to obtain a finite-dimensional multiplet space. 
Furthermore, Eq. (A22) allows to pair up raising and low- 
ering operators. That is, if S a is a raising operator, then 
with 

[S z , = -[S z , = -f za (S a )\ (A23) 

the operator (S a y changes the z-labels exactly in the 
opposite direction as S a in Eq. (A22). In this sense, 
(SaY = S-cr represents the corresponding lowering oper- 
ator. The actual definition of what is a raising or lowering 
operator is not entirely unique, as it depends on the spe- 
cific underlying sorting scheme of the z-labels adopted in 
root space. This does not matter, however, as long as 
the sorting is done consistently throughout. 38 ' 39 

In the presence of inner multiplicity a few complica- 
tions arise. Most importantly, an RLO usually will gen- 
erate a superposition in the m z <-fold degenerate state 



space in the resulting q z i , 

m,i 

S a \qn;q z a z ) = ^ \qn; q z -a z ,) (A24) 

a z ,=l 

with some coefficients s q \ ^ s ; 9 /c » , • As a consequence, the 
application of a raising operator So- followed by its com- 
plimentary lowering operator S^ onto a symmetry eigen- 
state, 

SlS a \qn;q z a z ) = ^ s l ^ a } aza ^\qn;q z a zl ) (A25) 

a,,=l 

with some other coefficients s q q z a -}a z a z , , does return to the 
same symmetry labels q z , yet not necessarily to the same 
state. If the resulting state in Eq. (A25) does not repli- 
cate the initial state \qn;q z a z ) up to an overall factor, 
then this allows to explore the other states in the de- 
generate subspace at q z . This is relevant for the decom- 
position of state spaces, where the resulting state as in 
Eq. (A25) needs to be orthonormalized in a consistent 
fashion with respect to the already explored states of the 
multiplet including the state \qn;q z a z ) (see Sec. Bl, for 
more detail on the numerical implementation). 

While all z-operators are required, e.g. for the defini- 
tion of the z-labels, it is usually not required to explicitly 
construct all of the RLOs, as some of these operators can 
be generated through a product of a smaller set of RLOs. 
As will be seen below in the case of SU(JV) or Sp(2n), 
the number of actually required RLOs can always be re- 
duced to the rank of the symmetry, i.e. the number of 
z-operators. This minimal set of RLOs will be referred 
to as simple RLOs, consistent with their general notation 
in the literature as simple roots of the symmetry. In a 
sense, these simple RLOs are the ones that induce the 
smallest shifts in the z-labels. 38,39 Again, their definition 
is not entirely unique, depending on conventions such 
as normalization of generators or what specific sorting 
scheme is applied to the z-labels. The simple RLOs still 
fully generate and connect the state spaces of any IREP. 
The underlying intuitive reason is that an r-dimensional 
z-label structure only requires r linearly independent vec- 
tors to explore its space (for a rigorous proof, see for ex- 
ample Refs. 38 and 39). Therefore given r z-operators 
{Z\, . . . , Z r }, it is sufficient to choose a specific subset of 
r raising operators {Si+, . . . , with the correspond- 

ing lowering operators 5,_ = (Si+Y . This reduction to 
simple RLOs is very useful in practice, yet does not re- 
strict the non-abelian treatment in any way. 

c. Maximum-weight state 

Consider some multiplet q of internal dimension d q for 
a given non-abelian symmetry group S of rank r. Then 
each of the d q states carries a set of r z-labels. When de- 
picted graphically as points in r-dimensional space, this 
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is called the weight diagram for the multiplet [for SU(3), 
for example, a collection of weight diagrams generated 
in an actual NRG run is shown in Fig. 12]. Since the 
z-operators are traceless, the values of the z-labels are 
naturally centered around the origin, i.e. q z = 0. Inner 
multiplicity, if present, decreases as a function of dis- 
tance \q z \ to the origin, such that the outermost points 
in a weight diagram always refer to unique states without 
any remaining multiplicity. By choosing a lexicographic 
ordering in the r z-labels, 25 the maximum weight (MW) 
is defined by 



max 



{«■}• 



(A26) 



The state with q z = qMW is called the maximum-weight 
state. This state is guaranteed to be unique to the multi- 
plet for non-abelian symmetries, 25 ' 38 ' 39 hence can be used 
as label for the entire multiplet, i.e. q — <?mw- While the 
state space \q z \ — m&x(\q z \) will not be unique, in gen- 
eral, since it refers to several states at the circumference 
of the weight diagram, max.{q z } does provide a unique 
set of z-labels. This underlines the importance of lexico- 
graphic ordering. 

As an example, consider the well-known spin SU(2). 
The states within the multiplet q are labeled by \qq z ) 
where q z = —q, . . . , +q identifies each state within the 
multiplet. This results in a one-dimensional weight dia- 
gram, with the multiplet itself labeled by the maximum 
weight states, q = max(q z ), indeed. 

Clearly, the q-labels for a multiplet themselves are also 
not entirely unique and hence depend on convention. In 
particular, if the rank of a symmetry group S is r > 1, 
the order of the z-operators themselves is a priori ar- 
bitrary. Hence there is a certain freedom in the order 
of the z-labels, which in return affects the definition of 
the maximum weight state. Given a certain order in the 
z-operators then, the lexicographic sorting of sets of z- 
labels is typically done in reverse order, i.e. starting with 
the last of the r label for a given q z . Moreover, hav- 
ing identified <?mw, this still leaves the freedom to use a 
linearly independent transform of ^mw &s label for the 
entire multiplet for consistency with literature. For ex- 
ample, for SU(3) [Sp(6)] this is discussed with Eq. (A34) 
[Eq. (A91a)], respectively. 



4. Example SU(iV) 

a. Defimng representation 

The symmetry SU(7V) is defined as the unitary sym- 
metry of an A^-dimensional space. The defining represen- 
tation, i. e. the IREP with smallest non-trivial dimension, 
is therefore given by N x N dimensional matrices. Since 
according to Eq. (A15b) all generators are traceless, only 
N — 1 diagonal z-operators exist, the diagonals of which 
form an ./V-dimensional orthogonal vector space that is 
also orthogonal to the diagonal of the identity matrix. 



The raising (lowering) operators are chosen as N x N 
matrices with a single entry of 1. anywhere in the up- 
per (lower) triangular space, respectively, away from the 
diagonal. For this, let 



fa) = (0,...,0,l (l) ,0,...,0) T , 



(A27a) 



with i e {1, . . . , N} be the ./V-dimensional cartesian col- 
umn basis vectors, and 



E; 



\ e i)( e i 



(A27b) 



the matrices of the related operator basis, which also 
contain just a single entry of 1 in their NxN dimensional 
matrix space, i.e. (Eij)i'j' — £«'<%'• Then the generators 
can be written as follows, 

c SU(iV) _ w J raising operator for i < j , . , 
rj I lowering operator tor % > j 



(A28b) 



These matrices are orthogonal as in Eq. (A15), while the 
(arbitrary) normalization was chosen such that, for con- 
venience, all entries are integers. The choice of generators 
for SV(N) in Eq. (A28) guarantees canonical RLOs, and 
thus simplifies the groups commutator relations w.r.t. z- 
operators exactly the way as indicated in Eq. ( A20) . This 
can be easily seen by observing that for a diagonal op- 
erator of the type (Z)ij = ZiSij, the matrix elements of 
the commutator with an arbitrary operator (S)ij = 
is given by 

\Z , S]ij — Siji^Zi 

that is, existing non-zero matrix elements in S are 
weighted by differences in diagonal elements of Z, while 
there cannot arise any new matrix elements unequal zero 
in [Z, S] as compared to S. Clearly, if S± only has a sin- 
gle non-zero entry as for the operators in Eq. (A28a), it 
follows [Z, S±] = const -Si, in agreement with Eq. (A20). 

From Eq. (A28a) above, a total of |JV(JV-1) different 
raising operators arise. However, not all of these are 
required to fully explore the multiplet space. Consider, 
for example, the subset of r = N — 1 raising operators 



{s s ; 



SV(N) X ^fqSV(N) QSV(N) 



,S£!JU (A29) 



which thus matches the rank r of the symmetry group 
SU(A r ) and thus also the number of z-operators. From 
repeated application of these operators, it is easily seen 
that the remaining raising operators not contained in 
Eq. (A29) can be generated. For example, Sf^^^ is 



generated by S- 



SU(AT) gSU(JV) 



2:) 



. Therefore, above minimal 



set of r raising operators with their hermitian conjugate 
set of lowering operators is sufficient, indeed, to explore 
all multiplet spaces. 
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b. The symmetry SU(3) 

The defining representation for SU(3) is chosen as in 
Eq. (A28), with the z-operators given by, 



(A30) 



Their diagonals can be collected as rows into a matrix, 
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(A31) 



the columns of which give the z-labels (z\,Z2) for the 
three states in the defining representation (see large black 
dots in Fig. 10). This represents the weight diagram for 
the defining representation. 

The canonical commutator relations as in Eq. (A20) 
yield the structure constants = f z<y for z € {1, 2} 

and a € {12,23,31}, 



/ = 



2 -1 -1 
3-3 



(A32) 



The columns in Eq. (A32) thus define the roots, i.e. the 
shift in z-labels when applying either 5 , 12 , S23, or S31, re- 
spectively. These vectors (roots) are depicted in Fig. 10 
by large thick arrows. Clearly, the three points in the 
weight diagram of the defining representation can be con- 
nected by these roots, equivalent to (repeatedly) applying 
raising or lowering operators. 

With the convention, that z-labels are lexicographi- 
cally sorted starting with the last z-label, i.e. sorting 
w.r.t. Z2 first and then z\, the three states in the defining 
representation are already properly sorted from largest 
to smallest [left to right in Eq. (A30)]. Furthermore, 
Eq. (A32) shows that S12 and £23 correspond to posi- 
tive roots, since (2,0) > (0,0) and also (—1,3) > (0,0). 
As their application makes z-labels larger, they represent 
raising operators, indeed, while S31 is a lowering opera- 
tor, all in agreement with Eq. (A28a). The third raising 
operator thus would be Si 3 with root (1,3) which, how- 
ever, is not a simple root and hence can be dropped. 

Finally, SU(3) still contains well-known SU(2) subal- 
gebras. That is, for example, by using S12 as a raising 
operator for the (a;, y) subspace together with its corre- 
sponding z-operator [Si2,S{ 2 ] =: 2Si 12 ' ) = Z\ while keep- 
ing the y component abelian, this shows that every line 
of points in the (zi,z 2 ) plane in Fig. 10 parallel to S12 
must correspond to a proper SU(2) multiplet. The same 
also holds for the two remaining permutations of (x, y, z) 
using S23 or S31 for the SU(2) subspace, instead. These 
SU(2) subalgebras clearly obey the standard commutator 
relations for SU(2). 




Figure 10. (Color online) Root space (z\,Z2) for SU(3). The 
three large black dots depict the weight diagram of the three- 
dimensional defining representation of SU(3). Large arrows 
indicate the shifts in z-labels due to the action of the rais- 
ing operators S12, S23, and the lowering operator S31, while 
dashed lines close to orthogonal to these arrows indicate lines 
of constant Si 2, S23, and S31 (color match with correspond- 
ing arrows). Dark symbols indicate accessible z-labels, while 
light crossed-out symbols are not accessible within SU(3) (see 
text). 



c. Symmetry labels for SU(3) 

The q-labels for a given IREP within SU(3) are derived 
from its maximum- weight labels (?mw = max{(zi, 22)} ■ 
With the z-labels additive through tensor products (see 
latter in the appendix), the z-labels of arbitrary mul- 
tiplets must be integer multiples of the z-labels of the 
defining representation. This immediately excludes the 
z-labels (points) in Fig. 10 that are crossed out. In par- 
ticular, with the columns of Eq. (A31) being linearly de- 
pendent, one may therefore use the columns of 




(A33) 



as basis for the maximum weight labels, for consistency 
with literature. 25,41 Given <7mw, the actual label of the 
multiplet then is determined by 



(gi,g 2 ) = z 1 • <7mw = 




(A34) 



This prescription makes the q-labels independent of 
the specific normalization conventions chosen for the z- 
operators. Furthermore, with z = z ■ (1, 0, — 1) T and the 
vectors in the columns in Eq. (A33) being positive by 
the adopted sorting scheme, this guarantees plain posi- 
tive integers for the multiplet labels q. These labels also 
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lie dense, in the sense that any (qi, 52) with qi > results 
in a valid multiplet. 

The defining representation with q-wiw = (1, 1) has the 
q-label (1,0). Its weight diagram together with many 
further examples for multiplets, as generated, in practice, 
from an actual NRG run using SU(3), are presented in 
Fig. 12. Note, however, that weight diagrams are mainly 
a matter of presentation of multiplets, while in practice a 
listing of z-labels suffices to describe the multiplet space. 



5. Decomposition into irreducible representations 

The generators of a specific symmetry group S repre- 
sent an irreducible finite set of operators {S a }, assumed 
to act in the full Hilbert space of a given physical system. 
Within (small) subspaces of the system, finite dimen- 
sional matrix representations can be constructed that 
obey exactly the same commutator relations as the gen- 
erators in terms of their structure constants in Eq. (A2). 
As such, a given matrix representation {R a } inherits all 
the properties of the generators. In particular, the matrix 
representation has the same number of operators as {S a } 
with a one-to-one correspondence in the symmetry label 
a. Therefore the z-operators as well as the RLOs share 
exactly the same interpretation within the D-dimensional 
carrier space of {R a }- 

Consider some arbitrary matrix representation {R a } 
that may have emerged, for example, from a tensor prod- 
uct space. As it operates in a D-dimensional carrier 
space, all of its matrices share the same dimension DxD. 
Assume a well-defined symmetry eigenstate within this 
space is available, to be called seed state, with a typically 
easy example being a maximum weight state. Then re- 
peated application of RLOs from the set {R a } generates 
a (sub) space which eventually describes a full symmetry 
multiplet, i.e. an IREP. By construction this subspace 
already diagonalizes the z-operators. Thus the z-labels 
are known, which also provides the q-labels for the mul- 
tiplet, e.g. by simply taking the maximum weight labels, 

qi = 9mw- 

If this multiplet q\ with d qi symmetry eigenstates 
spans the entire 7>-dimensional carrier space, then the 
matrix representation {R a } is already irreducible. If only 
a subspace of the 7?-dimensional carrier space was gener- 
ated, i.e. d qi < D, the matrix representation {R a } is re- 
ducible. Multiplet qi then defines a fully separated space, 
given the symmetry operations in {R a }. Combining the 
orthonormal state space of multiplet q\ as columns into 
a matrix Vi, the matrix representation {Ra-} can be cast 
into the space of multiplet qi, R a — > 7^ = V^RaVi, 
which thus constructs the irreducible matrix representa- 
tion iM for IREP qi . 

In case the 73-dimcnsional vector space is not ex- 
hausted yet, above procedure can be repeated with an- 
other seed state within the remainder of the vector space, 
generating further irreducible multiplets (72, <?3, • ■ until 
the D-dimensional vector space is fully exhausted. By 



combining the state spaces of the multiplets thus gen- 
erated, the resulting unitary matrix V = [Vi, V%, . . .] al- 
lows to block-decompose the original matrix representa- 
tion {R a } in terms of its irreducible representations, 

V^R a V = ($M q lM. (A35) 

Q 

where q runs through all IREPs J^. Note that a given 
IREP may be generated multiple times in the decom- 
position, which is indicated by the outer multiplicity 
M q e [0,1,2,...]. The presence of outer multiplicity 
therefore refers to the situation that M q > 1 for at least 
one q in the decomposition. In this case, also inner mul- 
tiplicity may occur, which refers to non-uniqueness of 
z-labels within an irreducible multiplet [cf. Eq. (A19)], 
both of which are specifically relevant, for example, for 
SXJ(N > 2) or Sp(2n > 2). 

As seen from above construction, the matrix represen- 
tation 7^ of IREP q is tightly connected to the symmetry 
multiplet q. In general, 7^ is unique only up to a global 
similarity transformation, as this does not affect commu- 
tator relations. By using its related multiplet state space, 
however, this space (i) can be chosen such that it diago- 
nalizes all z-operators, and (ii) can put into a well-defined 
order as provided, for example, by the lexicographic or- 
dering in the z-labels used to define the maximum weight 
state. Based on this basis, the matrix representation 7^1 
can be determined uniquely. This procedure on obtaining 
unique irreducible matrix representations will be adopted 
throughout. 

The decomposition in Eq. (A35), finally, can be done 
fully numerically along the same lines as already sketched 
above. Particular attention, however, must be paid to 
issues related to inner and outer multiplicity for over- 
all consistency. This will be discussed in more detail in 
Sec. Bl. 



6. Tensor product spaces 

Consider two irreducible matrix representations I^ qi ^ 
and 1^ of some non-abelian symmetry group S, with 
their matrix elements written in the basis of the symme- 
try eigenstates \qiqi z ) and \q2q2z) of the two IREPs q\ 
and 92 j respectively. The two multiplets are assumed to 
live in different spaces, so they can be joined through a 
tensor product, i.e. \qiqi z )\q2<l2z) = QM2z}- Then 

the generators of the symmetry in the combined space 
are defined in an additive fashion, which derives from 
the origin of the generators in infinitesimal symmetry 
operations, cf. Eq. (Al), 

R tot = jfa] g, jfc] + g, j[q 2 ] ( A36 ) 

Note that the additivity of the symmetry generators 
directly also implies the additivity of z-labels for non- 
abelian symmetries in general. And even if the non- 
abelian part of the SU(iV) symmetry is broken, e.g. re- 
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duced to an abelian symmetry with quantum labels q z , 
these are, of course, still additive. 

By construction, the tensor product representation 
{R^} in Eq. (A36) is also a representation of the sym- 
metry, as it obeys the same commutator relations within 
the combined system as the IREPs 1^ within their in- 
dividual space, 

[ijrtot )jR tot] = [ R ki] jR i^ 8 xte] + j>] g [ R M jR M] 

=y f i B lqi] =v f i R lq2] 

= E/->4 ot - ( A37 ) 

The product representation {i?Jj. ot }, however, is typically 
reducible. The resulting decomposition into IREPs is 
done exactly the same way as in Eq. (A35). 

The unitary transformation that rotates the product 
space \qiqiz] Q2l2z) into its combined symmetry multi- 
plets \qq z ; (a)) is given by the Clebsch-Gordan coeffi- 
cients (CGC), 

!<?<&;(«))= kiHz;q2q2z)(qiqiz;q2q2z\qqz;(a)), 

v v ' 

=c lq2 - ] , , 

(A38) 

with the shorthand notation C 1 ' 92 ^ , > for CGCs, con- 
sistcnt with the MPS tensors in the main body of the 
paper [cf. Eq. (10)]. Note that the CGCs implicitly also 
carry the multiplct labels qi, 52, and q. The index a has 
been added to account for possible outer-multiplicity [cf. 
Eq. A35] , in that for input multiplets q\ and qi the same 

output multiplet q can appear M g 91 ' 92 ' times, therefore 
a = 1, . . . , Afg 91,92 ' for a given q. If outer multiplicity 
is absent, the index a can be omitted, hence the round 
brackets around a in Eq. (A38) or Eq. (A39). 

As outer multiplicity also refers to different multiplets 
and hence state spaces, the Clebsch-Gordan coefficients, 
reflecting a unitary transformation, obey the general or- 
thogonality condition, 

E ^Ua^^)=^-i.^). (A39) 

This holds within the same output multiplet q, whereas 
the overlap between different output multiplets is strictly 
zero. While outer multiplicity is intrinsically connected 
to the underlying symmetry and hence to CGCs, in 
addition, this also affects the output multiplet space 
which must accommodate the additional multiplets [e.g. 
see QSpace discussion in the main text; note also that 
Eq. (A39) is completely analogous in structure to the 
orthogonality relation of ^-tensors as in Eq. (11)]. 



7. Irreducible operator sets and Wigner-Eckart 
theorem 

Consider a set of generators {S a } of some symmetry 
group S that a Hamiltonian H commutes with. Then 
all energy eigenstates of the Hamiltonian can be cate- 
gorized with well-defined quantum-labels, as indicated in 
Eq. (A17). In order to maintain an effective book keeping 
of the quantum labels when calculating matrix elements 
of operators, it must be possible to similarly categorize 
the operators themselves. Typically, the operators of in- 
terest are closely related to the Hamiltonian, i.e. consist 
of operators that also appear in the Hamiltonian or are 
composites thereof, such as creation, annihilation, occu- 
pation, spin operators, etc. Since the Hamiltonian can 
be properly constructed within the given symmetry, so 
can be its constituents. 

An irreducible operator set is constructed in a com- 
pletely analogous fashion as an irreducible state space, 
with an explicit example already derived for the spin op- 
erator in Eq. (All) using Eq. (A10). Consider the generic 
setup of a set of generators {S a } including RLOs. Then 
irreducible state multiplets can be generated through it- 
erative application of these operators, 

SM*) = s [ l°l\qq' z ), (A40) 

as in Eq. (A24), while ignoring inner multiplicity for the 
sake of the argument and having dropped the energy mul- 
tiplet index n for simplicity. Given an operator F, on the 
other hand, its transformation according to a symmetry 
is fully reflected in its commutator relations with the gen- 
erators of the symmetry. This is easily motivated through 
infinitesimal symmetry operations as in Eq. (A3b). The 
commutator relations, on the other hand, also emerge 
naturally when analyzing the effect of a generator of the 
symmetry acting onto a symmetry state \qq z ) that al- 
ready also has the operator F applied to it, 

S a -F\qj z ) = [S^FWqq'J+F-S^qq'J. (A41) 

The second term on r.h.s. clearly describes the symmetry 
properties of the state \qq z ), while the first term yields 
the transformation properties of the operator F which 
are independent of the carrier space. This is similar to 
what has already been seen in Eq. (A21) for the combined 
action of two generators. 

Now iff an operator set F q = {F^ } transforms exactly 
the same way as the state space of IREP q in Eq. ( A40) , 
that is 

then the operator set F q is called an irreducible opera- 
tor (IROP) set that transforms like the multiplet q. It 
carries the symmetry labels {q,q z ) the same way as an 
irreducible state multiplet does. 
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a. IROP decomposition 

In the case that a specific member of the IROP set is 
already known, then Eq. (A42) allows to generate the full 
IROP set exactly the same way as a state multiplet can be 
generated. Again, the maximum weight label determined 
the multiplet q that the IROP represents. This was ex- 
actly the procedure adopted, for example, to obtain the 
spin operator in Eq. (All). Furthermore, Eq. (All) also 
serves as a simple demonstration that the space of gen- 
erators itself clearly also can be cast into a single IROP. 
The corresponding multiplet is called the regular repre- 
sentation then. 

More generally, it is instructive to realize that irre- 
ducible operator sets (IROPs) and symmetry multiplets 
(IREPs) can be treated on a nearly equal footing. In 
particular, the notion of proper orthonormalization of 
state spaces can be directly applied also to IROP sets, 
up to a global normalization factor. This is motivated 
by the observation, that given a scalar multiplet |0) for 
which \F q ) = F q q , |0) ^ 0, i.e. does not vanish, then \F q q , ) 
represents the multiplet vector space for IREP q. With 
proper overall normalization of the IROP F q , it follows 

The last equation also holds, if the scalar multiplet |0) 
is replaced by an arbitrary other symmetry eigenstate 
\qq z )- For matrix representations of IROPs and oper- 
ators more generally, this motivates the scalar or inner 
product for two matrices as in Eq. (A14). Thus equipped 
with scalar product and norm for matrices, an IROP de- 
composition can be done exactly the same way as the 
multiplet decompositions for symmetry multiplets start- 
ing from a specific symmetry eigenstate (IROP compo- 
nent). This is important, in particular, in the presence 
of inner multiplicity in the multiplet of an IROP for con- 
sistency with the Wigner-Eckart theorem. 

b. Wigner Eckart theorem 

It follows from Eqs. (A41) and (A42), that the states 
resulting from the IROP {F q ^ z } applied to a multiplet 

\q2q2z), 

Sa ■ F£Jq 2 q 2z ) 

= [S„,F£J \q 2 q 2z ) + F^- ^jcM^l , 

q lz q lz q lz q 2z q 2z 

transforms exactly the same way under given symmetry 
as a tensor product of two state multiplets, 

So- ■ \qq z )i\qqz)2 

= S la \qq z )^ ® \qq z ) 2 + \qq z )i ® S^alq^uh ) 

w»h =< 1 ?'}, Mzh 

q lz q lz q 2z q 2z 



using Eqs. (A36) and (A40). Therefore the action of an 
IROP F qi onto the state space of an IREP q 2 shares ex- 
actly the same algebraic structure in terms of symmetries 
like the product space of the two multiplets qi and q 2 ■ 

This motivates the Wigner-Eckart theorem. With 
the definition of the Clebsch-Gordan coefficients in 
Eq. (A38), it is thus clear that up to scalar factors de- 
pending on the normalization of the operator set, the 
same CGCs also apply for the state space decomposition 
arising out of Fjj* \q 2 q 2z ). In particular, it follows for the 
matrix elements of the operator w.r.t. a given state space, 

(qq z ; {a)\F q l\q 2 q 2z ) = (qq z ; (a)| • {P q q l\q 2 q 2z )) 

= (q;(a)\\F^\q 2 )-C q l± {aV 

(A43) 

where, again, a accounts for possible outer multiplicity. 
(q; (a)\\ A qi \\q 2 ) is called the reduced matrix element. It 
is entirely independent of the z-labels, i.e. the internal 
structure of the IREPs qi, q 2 , and q. 

The first line in Eq. (A43) specifies the adopted con- 
vention for matrix elements given the Wigner-Eckart the- 
orem: the operator is acting to the right ket-state, the 
symmetry labels of which are combined. The result- 
ing object is contracted with the bra-states. This is 
important for consistency, since the IROP Fi is sub- 
tly different from the IROP (F^)". Therefore one must 
be careful with expressing a matrix element through 
(qq z \F\q 2 q 2z ) = (q 2 q 2z \F^\qq z )* . Even though usually 
(Fq z y ~ (F') q _ qz , further signs may bee needed to en- 
sure for consistency within the Clebsch-Gordan coeffi- 
cients [e.g. see discussion around Eq. A49 later]. 

8. Several independent symmetries 

A physical system often exhibits several symmetries. 
Each of the A = 1, .. .,ng symmetries S x is completely 
described by its own set of generators {5*}. As these 
symmetries act independently of each other, this implies 
that their generators must commute, 

=0 for A ^ A'. (A44) 

This allows to assign independent quantum labels (q x q^) 
with respect to each individual symmetry [cf. discussion 
following Eq. (1) in the main paper]. On the multiplet 
level, the symmetries are given compactly by the com- 
bined q-labels, q = (q , q 2 , . . . , q ns ), while similarly their 
z-labels are given by q z = (q\, q 2 , . . . , q z s ). Here the el- 
ementary multiplet labels q x and q^ can already consist 
of a set of labels themselves, the number of which is de- 
termined by the rank r of the respective symmetry <S A 
[cf. Eq. (A18)]. 

When a non-abelian symmetry is broken, it it is re- 
duced to simpler subalgebras. In particular, it may be 
reduced to its abelian core of z-operators (Cartan subal- 
gebra). For example, consider the rotational spin SU(2) 
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symmetry. This symmetry can be broken by applying a 
magnetic field. The system still maintains, however, a 
continuous rotational symmetry around the axis of the 
magnetic field, leaving the q z symmetry intact, while 
the multiplet label q becomes irrelevant. Similarly, if 
particle-hole symmetry (see later) is broken, only the 
abelian quantum number of total charge [i. e. the z-label] 
remains. 

Abelian symmetries therefore fit seamlessly into the 
general non-abelian framework outlined in this paper. 
With the multiplet label q irrelevant, the q z are promoted 
to the status of a q-label, instead, with no z-labels re- 
maining [with all multiplets being one-dimensional, the 
z-labels are no longer required, i.e. can be set to zero, for 
simplicity] . This then allows to write the abelian symme- 
try in terms of trivial scalar Clebsch-Gordan coefficients. 
The latter, nevertheless, are important as they account 
for the proper addition rules w.r.t. the abelian z-labels, 



9. Symmetries in physical systems 

In the following, several examples of symmetries in 
simple physical systems will be given, with the associated 
spinors and irreducible operator sets explained in detail. 
In particular, this concerns fermionic systems with spin 
or particle-hole symmetry. 

For the model Hamiltonians in strongly correlated 
electron systems, correlation through interaction plays 
an important role, while the terms describing interac- 
tion typically preserve certain underlying global symme- 
tries. Since the arguments of demonstrating symmetries 
of a specific Hamiltonian, however, are rather similar, in 
general, it suffices to consider a simple non-interacting 
Hamiltonian. Simple issues related to interactions are 
discussed with Eq. (A58) below. 

For simplicity, therefore much of the following discus- 
sion will be exemplified in terms of the Hamiltonian of a 
plain spinful fermionic tight-binding chain, 

# = E**E (4A+i, CT + H.c.), (A46) 
k <y 

S v ' 

where c ka creates a particle at site k with spin a € {f, l}. 
The Hamiltonian in Eq. (A46) has spin-independent hop- 
ping amplitudes tk, hence possesses spin-SU(2) symme- 
try, SU(2) sp i n in short. Furthermore, it is particle-hole 
symmetric, implying particle-hole SU(2) symmetry, also 
called charge-SU(2) symmetry, or SU(2) c h a rgc in short. 



a. SU(2) spin symmetry 
Using the two-dimensional spinor 

h,u = (y) (A47) 

for each site k, the Hamiltonian in Eq. (A46) can be 
rewritten as 

# = •^,*+i + H - c -). (A48) 

k 

where the sum over a was incorporated in the scalar 
product of the vector of operators in ips.k- Clearly, the 
two-dimensional scalar product is invariant under an ar- 
bitrary unitary two-dimensional transformation U, i.e. 
ft k ip k+1 — (Ut^kY (U^pk+i), thus exhibiting spin-SU(2) 
symmetry. The spinor in Eq. (A47) is defined in a site 
specific manner. When concentrating on a single site, 
therefore the index k can be dropped for convenience. 

The generators of spin-SU(2) symmetry are con- 
structed in terms of the two-dimensional defining rep- 
resentation of {S„} = {S+,S Z ,S-} [cf. Eq. (A7)]. These 
can be written as operators (distinguished by the hat) 
through second quantization in the full Hilbert space, 

S a = ft s S a ip s , 

which up to prefactors leads to the spin IROP S 1 = 
{-^S + ;S Z ;+^S_}, already derived in Eq. (All). 
The raising operator, for example, is given by 

which flips a down-spin to an up-spin for given site. Sim- 
ilarly, the z-operator is given by 

S z = ft (^r z )ip = |(c|c t - cjcj_) = \ (?V h). 

Furthermore, [S+,cj] = [c|c^,c|] = cj shows that the 

spinor ft s already represents an IROP for the q = \ 
multiplet of SU(2) sp i n , 

(ft s )^ = (j) . (A49a) 

This is already properly sorted w.r.t. z-labels, in that 
the second component correspond to the lower q z = — i 

element of the multiplet, since [S z , c|] = (—5) ■ c|. 

In contrast, the IROP for the spinor -05, i.e. without 
the dagger, is similar, yet has subtle differences. In par- 
ticular, with [S+jCj.] = [ctci,c.J = — c,, the role of spin 
within the multiplet is reversed, i.e. q z — ¥ —q z , while also 
an additional sign is acquired, 




(A49b) 
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This extra sign is important in context of the Wigner- 
Eckart theorem in Eq. (A43), where the particular or- 
der of first applying, i.e. combining an operator with the 
ket-state is directly related to the order in the Clebsch- 
Gordan coefficients. This is convention, of course, but 
consistency is crux. 

In terms of the proper IROPs in Eqs. (A49), finally, 
the Hamiltonian in Eq. (A46) can be written in either 
IROP while, however, one must not mix them, 



H = X>([(&0 [1/2] ] f • ^s) [1/2] +H.c.) (A50a) 

k 

= E^((^) [1/2l -[(^) [1/2l ] t +H.c). (A50b) 



The second line is essentially the same as the spinor ex- 
pression in Eq. (A48), yet with the difference, that here 
the underlying IROP structure has been pointed out ex- 
plicitly. 



b. SU(2) particle-hole symmetry for spinful system 



The particle-hole symmetry SU(2) c h arge of the Hamil- 
tonian in Eq. (A46) can be made apparent in a similar 
way as for the spin symmetry above. Consider the spinor 
in the charge sector, 



1pC,ko 



with alternating phases Sk = (~l) fe along the chain in 
Eq. (A46). Again, the Hamiltonian can be written as 
sum over scalar products in the spinors, 




cLC; , ... - r, j- 



i '.i.n- ' Vc.k : 1,(7 — E v ' A + L.T 

= hk ,fc+i) 



k-tr^k+lj—tr) 



suggesting another underlying SU(2) symmetry. Note 
that the alternating sign Sk is crucial to recover the cor- 
rect hopping structure in Eq. (A46). Given the spinor in 
the charge sector, the raising operator becomes 



which, up to a sign, is the same for both spins. It is there- 
fore sufficient in the charge sector to consider a spinor for 
one specific a in ipc,ka only. Therefore, again concentrat- 
ing on a single site and hence dropping the site index k, 
now with fixed a =t, the spinor in the charge sector is 
given by, 



(A51) 




The associated raising operator becomes 



(A52) 



which now creates a pair of particles with opposite spin, 
while the z-operator is 



= i(n t + (A53) 

With n = fU + hi, the z-operator C z counts the to- 
tal charge on given fermionic site relative to half-filling. 
With 



- I/at j 



fit'' 



[C+,c\\ = [scjct.c;] = scl, 



at 



(A54a) 

this allows to construct the q = | IROPs for SU(2) c h a rge> 



Wo) 



1/2] 




(A55a) 
(A55b) 



again associating the lower component with the q z — 
-\ element of the q = (1/2) multiplet [cf. Eqs. (A49)]. 
An irrelevant overall minus sign has been applied to the 
spinor in Eq. (A55a) for later convenience. With this, 
the hopping term in the Hamiltonian in Eq. (A46) can 
be rewritten in terms of the scalar products 

^,*+i=[(^*) [1/a] ] t -(^c > * +1 ) I1/2] 

+ m k f /2] ]' ■ (^, fe+1 ) [1/2] (A56) 

The spinors in the charge sector do mix spin compo- 
nents, which essentially also requires full spin symme- 
try [see later discussion of symplectic group Sp(2m) in 
Sec. A 10]. More importantly, the construction of the 
SU(2) c harge symmetry allows it to fully commute with 
the spin-SU(2) symmetry introduced earlier, 



[S g ,C z ] = \ [h t - hi,h 



1]=0 



-£+£+ 



[S+,C Z ] — l[c^c^,c^c^.} + 2 [c^.Cp Cj_Cj_] 



= 



= 



-etc 

t 4- 



[S+, C + ] = s[c\c v c|c|] = sc\c\ = 



(A57) 



That is, the two symmetries act completely independent 
of each other and thus can coexist simultaneously, writ- 
ten as the overall symmetry SU(2) spin <g) SU(2) c har g e- 

If interactions are present in the system, such as local 
Coulomb interaction Uh^h^, then the particle-hole sym- 
metric regime requires a specific altered onsite energy 
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relative to the chemical potential. With n 2 
h = TVf + n_i, it follows that n^n^ 
and therefore 



h a , and 
i(n-l) 2 + i(n-l), 



Sdh + Uh-fhi = (ed + ^){n — 1) + ^(n — l) + const. 

(A58) 

The first term on the r.h.s. is proportional to the C z oper- 
ator, which thus acts like a magnetic field for SU(2) sp i n . 
Therefore for full particle-hole symmetry to hold, this 
term must be zero, which requires Ed — — ^ . In particu- 
lar, in the absence of interaction, this implies Ed = 0. The 
actual interaction term, i.e. the second term on the r.h.s. 
in Eq. (A58), also resembles C z . Yet it is quadratic, and 
for this it also holds, C z = C 2 = C 2 . Therefore, this term 
can actually be written as C 2 which itself, like spin S 2 
for SU(2) S pi n , represents the Casimir operator for SU(2), 
and thus is compatible with SU(2) spin (g> SU(2) chargo . 

The actual IROP for particle creation and annihila- 
tion given SU(2)fg = SU(2) spin ® SU(2) chargc symmetry 
can be generated using above symmetry operations. This 
generates a four-dimensional spinor. As it turns out, the 
resulting IROP is the combination of the two IROPs gen- 
erated in the spin symmetric case in Eqs. (A49) as well 
as in the particle-hole symmetric case in Eqs. (A55), 



,/,L2'2j 

res 



fsc\\ 



(A59) 



The multiplet labels [|, |] will be derived with Eq. (A61) 
below. The signs for the individual components in above 
IROP have been properly adjusted, considering that the 
raising operator in the charge sector itself, Eq. (A52), 
carries the alternating sign S(m = (— 1) . For example, 

commuting C+ onto the fourth component, yields the 
third component of the spinor ^>cs [c/. Eq. (A54a)], while 
commuting S + onto the third component yields the first 
component, and so on. Again, keeping track of the alter- 
nating sign Sk = (— l) fc is crucial to recover the hopping 
structure in Eq. (A46), 



CS,k 



J CS,k+l 




l>k,k+l, 



(A60) 



The full tight-binding Hamiltonian simply becomes H = 

J2k ^fc('0cs.fc)^V'cs fe+i wnere the hermitian conjugate 
term has been incorporated already in the spinor struc- 
ture. This also reflects the irrelevance of taking the her- 
mitian conjugate version of the IROP in Eq. (A59) as 
this results in essentially the same object after properly 
reordering of its components and taking care of signs. 



With Eq. (A60) being a scalar product in a four- 
dimensional spinor space, one may be tempted to think 
that a plain tight binding chain actually has a non- 
abelian symmetry with a defining representation of di- 
mension 4. This cannot be the symmetry SU(4), how- 
ever, since SU(4) has rank-3 and thus requires three com- 
muting abelian z-operators. The symmetries discussed 
here, however, only have two abelian z-operators, namely 
total spin and total charge. The symmetry that ap- 
pears compatible with this scenario, at second glance, 
is the symplectic symmetry Sp(4) [see Sec. A 10 be- 
low]. Nevertheless, even the latter can be excluded, since 
raising and lowering operators are severely constrained 
by the fact that the creation and annihilation opera- 
tors appear in pairs for the same fermionic particle in 
the IROP of Eq. (A59). Consequently, quadratic oper- 
ators of the type (c CT )^cJ = i&^yCa = are immedi- 
ately excluded. With this, the symmetry of the spinor in 
Eq. (A59) has to remain the product of two symmetries, 
i.e. SU(2)fg = SU(2) spin ® SU(2) chargc . 



Having determined the IROP -0CS in Eq. (A59) by re- 
peated application of RLOs S± and C±, the z-labels for 
each of the four components, on the other hand, can be 
determined through the z-operators C z = | (clcL — c, c{) 

and S z = \ (c*C| — c| c, ) , resulting in the z-labels q z = 
{C z , S z ), respectively. The results are summarized in the 
following table. 



[ z-operator, IROP component ] (C Z ,S Z ) 



[4,4] 

[S z , sc|] 

[4, h] 

[S z , Cj.] 

w 



[C z ,sc 
[S z ,sc 



[4,-a T ] 



+K4) 
+K4) 

-K h) 
+l( h) 

+K4) 

= -|(-2t) 

= -l(-fir) 




(A61) 



These z-labels demonstrate that both the charge and spin 
multiplet contained in ipcs corresponds to a q = g multi- 
plet. The maximum weight state has the z-labels (|, |), 
which thus labels the spinor, as was already indicated in 
Eq. (A59). 

Similarly, the local state space of a fermionic site must 
be organized consistent with the SU(2) sp j n ® SU(2) c h ar ge 
symmetry above. The local state space consists of the 
empty state |0), the singly occupied states |f) and \l), 
and the doubly occupied state s |f!). Note that the sign 
in the last state is crucial, as it is generated by the raising 
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operator C+ acting on the empty state |0). In summary, 



state 


c C z 


s s z 


10) 

it> = 4io) 

|+) = cj|0) 
s\n)=sc\cl\0) 


1 1 

2 2 




2 '2 




2 '2 

1 1 

2 2 

0. 



Therefore the local four-dimensional state space of 
fermionic site is spanned by the two multiplcts, q = 
(C, S) € {(0, |), (|, 0)}. If particle-hole symmetry is bro- 
ken yet particle number still preserved, then 2C Z from 
the middle column describes the total number of parti- 
cles relative to half-filling [cf. Eq. (A53)]. 

With the ordering convention of state labels being 
|C, C z ; S, S z ) and a e {t,4J = {+1, 1}, the non-zero 
matrix elements of the 4-component spinor in Eq. (A59) 
can be calculated. For example, 

+sa = {-a\c a ■ s\ U) = s(U \s ■ sc\ ■ \ - tr)W 
= s(i,±i;0,0^ i+i iaJ0,0;i,^) 

^2 2 '2 2- 1 

= s<i,±i|i, M ;0 ,0)<0,0|i f;i,^><i,0||^||0,i) 

v v " v ' 

= 1 (charge) = ±g (spin) 

v2 

=> (1,011^110,1) = ^. 

The order inversion of the matrix element in the first 
line was used since the spinor ^cs m Eq. (A59) contains 
the creation operator sc\ and not its hermitian conju- 
gate. The overall complex conjugation however, is 
irrelevant since all matrix elements are real, hence the 
notation of putting the asterisk in brackets. 

The second non-zero reduced matrix element can be 
calculated in a similar fashion, 

l = {Q%\a)=s{cr\-scl • |0)« 
= «<0,0;i,f|$ .1+1 iaJ|,^;0,0) 

*-2 2 '22> 

= , s( 0,0|I,±i;I,^)(I,||I,f;0,0)(0,i||^||i,0) 

v v *■ V ' 

=+i (charge) =1 ( s P in ) 

(0,iW|,0)=«>/2 

Overall, this leads to the reduced matrix elements in the 
charge-spin sectors (C,S) £ {(0, |), (|,0)} 

4l 2 ' m = (^ S f )• (A63) 

Note that although the spinor in Eq. (A59) has four com- 
ponents, i.e. is of rank-3, on the reduced multiplct level in 
Eq. (A63) the spinor becomes a two-dimensional object 
as expected from an IROP. The further internal struc- 
ture is entirely taken care of by rank-3 Clebsch Gordan 
coefficients, which have been omitted in Eq. (A63) [for 



[1/2 1/21 

a full description of ip C s ' including Clebsch- Gordan 
coefficients in terms of a QSpace see Tbl. C5]. 

The operator in Eq. (A63) is non-hcrmitian. In the 
context of two-site hopping, however, this nevertheless 
leads to a hermitian term in the Hamiltonian, as required. 
Indicating the local symmetry eigenspace for site k by 
\<j) k , the matrix elements of the hopping term in the 
tensor-product basis \<j k +\) |cfc) (in this order, assuming 
site k + 1 is added after site k) are given by 

(^IK+i|^ fe+1 K+i)K) = i>{ ® WWi, (A64) 

where the ^'s to the right without the hat denote the ma- 
trix elements in the local \a) basis. Note, that fermionic 
signs apply, when tp k is moved, for example, to the left 
of (ctj. +1 |, such that the tensor-product on the r.h.s. of 
Eq. (A64) contains [zi/>]fc+i rather than ipk+x, where i k = 
( — is diagonal in \a k ) and adds signs corresponding 
to the number of particles in \a k ). Note that with the par- 
ticle number being related to C z — \{n — 1), the operator 
z is well-defined in terms of the symmetry labels. It is a 
scalar operator, since (— = (— = (— 1) 4C *, 
hence does not alter the Clebsch-Gordan content of the 
operator ip but rather acts on the multiplet level only. 
For the hopping ipl.'4' k+1 between two nearest-neighbor 
sites, Eq. (A63) finally leads to 

s * 

/ -s k+1 V2] 
\V2 ) 

written as a plain tensor product on the level of the mul- 
tiplet spaces of two fermionic sites. For the sake of the 
argument, the product space here is not yet described in 
terms of proper combined symmetry multiplcts of sites k 
and k + 1 . 

With s k = (— l) fe , H k:k+ i in the last equation clearly 
yields a hermitian object for all iterations. For example, 
for even k, the hopping term is given by 

[k even] _ ( ( \M 

similar in structure to a hermitian object of the type 
t x <£> t x in terms of Pauli matrices, while for odd fc, 

[fc odd]_ ( V2\ ( -V2\ 

similar in structure to the hermitian (ir y ) <g> (ir y ) = —r y ® 
T y . Hence for every even (odd) iteration, one has a t x ®t x 
(t v (g) T y ) structure, respectively, a prescription that is 
periodic with every pair of iterations. This intrinsic even- 
odd behavior is not specifically surprising, considering, 
for example, that two particles are needed to return to 
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the same charge quantum numbers related to particle 
hole symmetry. 

In summary, using Eq. (A60), the hopping in the 
Hamiltonian is given by 

# = (A65) 

k 

In a sense, the net effect of incorporating spin SU(2) 
was to eliminate the spin index, while incorporation of 
particle-hole SU(2) eliminates the hermitian conjugate 
term in the Hamiltonian. Together they reduce the four 
terms initially required for a single hopping in Eq. (A46) 
to the single scalar term r 4>\'4>k+v 

c. Particle-hole SU(2) symmetry for several channels 

The alternating sign in the raising operator Cu.+ = 
scjcj in Eq. (A52) defines the doubly occupied states 
as ||; 4^) = Ck,+ \\ \ = scjcjjO); for even sites, s — 
+1, therefore ||;-^) = cjcjjO). For odd sites, on the 

other hand, ||; = — clcT jO) = cjcj |0). In practice, for 
consistency, usually a certain well-defined fermionic order 
is adopted. Above raising operator Ck,+ thus suggests 
that it may be useful to reverse the fermionic order of 
every other site for the local state space included there. 

Fully reversing the fermionic order of a given site fc 
with several fermionic channels i = 1, . . . , m implies for 
the matrix elements for particle creation or annihilation 
operators, 

c k,itr ~^ ^k,ia = ~ z k c k,icr- 

This transformation is equivalent to a unitary transfor- 
mation local to site fc. Similar to Eq. (A64), z% = (—l) nk 
with n k = J^ia ™k,i<j and h kyla = c\ l(J c kM again takes 
care of fermionic signs for the full multi- level site fc. Be- 
ing a scalar operator, z k is independent of the fermionic 
order. 

Now consider the effect of flipping the fermionic order 
for the odd sites in the tight-binding chain that carry the 
sign s = — 1, assuming particle-hole symmetry in every 
channel. For a specific channel, this (i) takes away the 
sign in the raising operator Ck t +, and (ii) implies, for 
example, for the 4-component spinor in Eq. (A59) for a 
single channel, 









f+zc\\ 




V'CS.fc odd = 


-4 




+4 

\-ZCf) 


= Zk1pCS,k even, 



having intermittently dropped the index k for readabil- 
ity. Therefore, up to the local operator z k which assigns 



fermionic signs to the full Hilbert space of a local site, 
the matrix elements of the spinor for the odd sites are 
exactly the same as the matrix elements of the spinor for 
even sites. Therefore with vb taken as the spinor for even 
sites in the chain, the required spinor for odd sites be- 
comes zip. Together with the additional fermionic signs 
in the nearest-neighbor hopping term as already encoun- 
tered in Eq. (A64), the hopping structure hk, k +i of the 
tight-binding Hamiltonian in Eq. (A65) becomes, 

vj\®[z ■ (ztp)]k+i =ipl® ip k+1 for k even 
(zip)l <£) [z ■ i/j]k+i = (zt/Ol ® ( z ^)k+i for k odd. 

This result generalizes to any number of channels with 
particle-hole symmetry. As such it much simplifies the 
structure and thus the treatment of the two different 
kinds of spinors for even and odd sites, respectively, that 
had been required initially. 



d. Symmetric three-channel system 

Consider the generalization of the spinful one-channel 
setup in Eq. (A46) to a spinful three-channel system, 

m— 3 
k i=l a 

" v ' 

where c\ ia creates a particle at site k in channel i with 
spin a. This model is relevant for the system analyzed 
in the main body of the paper where the specific num- 
ber of three channels, for example, originates from the 
underlying orbital band structure in terms of a partially 
filled cZ-shell. The Hamiltonian in Eq. (A67) can also 
be complemented with interaction terms that are com- 
patible with the symmetries discussed in the following. 
This can include onsite interaction U at half-filling [cf. 
Eq. (A58)], or uniform local Hund's coupling Ju [e.g. see 
Eq. (25b)]. Here, however, the focus of the discussion is 
on symmetries, for which the Hamiltonian in Eq. (A67) 
suffices. 

The Hamiltonian in Eq. (A67) possesses SU(2) spin 
symmetry, SU(2) particle-hole symmetry in each channel, 
and also SU(3) channel symmetry, while not all of these 
symmetries necessarily are independent of, i.e. commute 
with each other. All of these symmetries can be defined 
within the Hilbert space of a local site, hence again focus- 
ing the discussion on a single site k in the following, while 
dropping the site index fc, for simplicity. For each of the 
three channels, the associated spinful fermionic level is 
represented by the four states as in Eq. (A62), leading to 
a total of 4 3 — 64 state for a given site. 

The total spin-SU(2) symmetry of a site is described 
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by the generators 

i 



= e4 

i 



(A68) 



where Si a represents the spin operators for the fermionic 
level i, with S z = ^[S+,S\], as expected for SU(2). 

The particle- hole symmetry exists for every channel i, 
and is described by the SU(2) symmetry, 



C i+ — sc^c^ 

Ciz — \ [Ci+,C i+ \ — \ {fii — l), 



(A69a) 



which includes the same sign- factor = (— 1) as in 
Eq. (A52) to correctly represent the hopping structure in 
the Hamiltonian Eq. (A67). The total charge relative to 
half-filling is given by (up to a factor of 2) 



(A69b) 



Finally, the channel symmetry is given by the minimal set 
of two raising operators {T+,U + } = {S^,^} together 
with the z-operator {T Z ,Y} = {Zi,Z 2 } as introduced 
through Eqs. (A28) in Eq. (A30), 



U+ = J2 4<rC 3(T , Y = X)(»i<7 + h2a - 2n 3o 



(A70) 



Here the notation for the generators of SU(3) has been 
changed to another notation frequently also found in lit- 
erature, so these generators can be better distinguished 
from the generators for spin and particle-hole symme- 
try. In particular, the operators T + and T z generate an 
SU(2) subalgebra, that is linked to the full SU(3) symme- 
try through the generators U + and Y . The normalization 
of the z-operators, however, has been chosen consistent 
with Eqs. (A28), such that plain integer matrix elements 
arise. 

The spin symmetry clearly commutes with the particle- 
hole symmetry in each channel, which follows from the 
previous one-channel discussion in Eqs. (A57). There- 
fore it remains to analyze the compatibility of the SU(3) 
channel symmetry. All z-operators clearly commute. For 
the SU(3) raising operators, it follows with respect to the 
spin symmetry, 
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with a similar expression for [/+ instead of T + with a 
shift in the channel indices. Note that in order for the 
r.h.s. to vanish, the sum over the spin a is essential 



which shows the importance of the summation over a 
in Eqs. (A70). As a consequence, the SU(3) channel 
symmetry in Eq. (A70) commutes with the SU(2) spin 
symmetry, indeed. 

The compatibility of the SU(3) channel symmetry with 
the SU(2) particle-hole symmetry, however, cannot be 
established, since 



[T+ , 



s5 i2 {& n c. 
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(A72) 



cannot be made to vanish for all channels i at the same 
time. Therefore the non-abelian channel and particle- 
hole symmetries cannot coexist independently of each 
other. Nevertheless, the generators of each individual 
symmetry do commute with the Hamiltonian, which thus 
suggests a larger symmetry, with Eq. (A72) already indi- 
cating one of the additional generators. As it turns out, 
this symmetry is Sp(2m) with m the number of channels. 
This symmetry will be introduced and discussed in the 
next section. 

By reducing the non-abelian particle-hole symmetry 
to its abelian conservation of total charge, however, this 
abelian symmetry does commute with the SU(3) channel 
symmetry, 

i a.i<j ! 



In order to get a commuting abelian charge symmetry, 
the z-operators for the channel-specific particle-hole sym- 
metry must be summed over all channels i. With all com- 
muting symmetries combined, this leads to the overall 
symmetry SU(2) spin ® U(l) ch argc ® SU(3) cha nnei, consist- 
ing of the SU(2) total spin symmetry in Eq. (A68), the 
abelian total charge of the system in Eq. (A69b), and the 
channel SU(3) symmetry in Eq. (A70). 

A more conventional symmetry setup can be obtained 
by giving up the channel SU(3) symmetry. Bearing in 
mind that the channel-specific SU(2) particle-hole sym- 
metries commute with total spin, this also allows the 
symmetry setup SU(2) spin ® SU(2)f h 3 argc . 

The symmetry combinations above can be motivated 
also by a simple counting argument with respect to con- 
served abelian quantum numbers. Note that the pre- 
served abelian quantum numbers in the Hamiltonian 
Eq. (A67) are the particle number in each of the three 
channels together with the total spin S z . This results in a 
total of four z-operators, and thus four z-labels. Now, by 
including non-abelian flavors, the number of z-operators 
clearly cannot increase, but will remain the same. To- 
tal spin has one z-operator, the channel SU(3) symme- 
try has two z-operators, and the channel-specific particle- 
hole symmetries have three z-operators, which combined 
results in 1+2+3 = 6 z-operators. This set of z-operators 
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therefore cannot be independent of each other, as already 
seen in the earlier discussion. Yet, in fact, both of the al- 
ternative symmetry setups above do have a total of four 
z-operators. For SU(2) spin <g> U(l) cha rge ® SU(3) channel, 
these are 1 + 1 + 2 from spin, charge, and channel symme- 
try, respectively, while for SU(2) spin <g> SU(2)gf arge these 
are 1 + 3 from spin and each channel. 

For the symmetry setting SU(2) sp i n ® U(l) c har g c ® 
SU(3) c hannei then, the hopping term in the Hamiltonian 
in Eq. (A67) is described by a 6-component IROP 
(annihilation operators for spin-up and spin-down com- 
bined), that can be obtained, for example, numerically 
as described in Sec. B2. This leads to 

h k , k+1 = ■ 4 61 + H.c. (A73a) 

In contrast, for the second symmetry setting SU(2) sp i n ® 
SU(2)f h 3 argc , the IROPs $J required for the hopping 
term are already exactly the 4-component spinors in 
Eq. (A59), i.e. one for each individual channel, i = 
1,...,3. The hopping in the Hamiltonian is thus de- 
scribed by 

m— 3 

hk,k+i = J2 4tM+i,i> (A73b) 

i=l 

i.e. without the hermitian conjugate part as this is al- 
ready included through particle-hole symmetry. Fur- 
thermore, note that particle-hole symmetry also acquires 
even-odd alternations for the spinors along a chain [see 
App. A 9b]. 

10. The symplectic group Sp(2m) 

All Hamiltonians considered in this paper are time- 
independent, hence obey time-reversal symmetry. Time- 
reversal symmetry then is described by an anti-unitary 
operator T — T, y K, 40 that includes a standard unitary 
operation E y together with the operator K, which stands 
for complex conjugation [the notation of E y has been 
chosen for latter convenience; see Eq. (A77) below]. The 
time-reversal operator obeys T 2 = ±1, where for spin- 
half particles, such as electrons as considered throughout 
in this paper, it holds T 2 = — 1. The latter is important 
for the symmetry Sp(2m), since it implies that the uni- 
tary E y must be antisymmetric. This follows simply by 
looking at the matrix elements of the time-reversal op- 
erator for arbitrary states \a) and \b) in some real basis 
i, 

(a\fb) =J2 a i(^y)ij b *ji 

yet it also holds, 

(a\fb) = (f^b\fa) =-Y f b*(Z y ) ji a*. 

-1 hj 



As this applies for arbitrary states \a) and \b), this shows 
that, given T 2 — — 1, the unitary E y must be antisym- 
metric, indeed. 

Since a time-independent Hamiltonian obviously com- 
mutes with the time-reversal operator, it follows that 
all eigenstates of the Hamiltonian can also be written 
as eigenstates of the the time-reversal operator T. As 
a consequence, all unitary symmetry operations G — 
exp(i cia-So-) can be constrained to unitaries which 
also leave the time- reversal operator invariant. That is, 

f = GTG- 1 = G± y KG^ 
=g t k 

=>■ ± y = G± y G T . (A74) 

For the generators So- of a symmetry group this implies 
(e.g. by expansion of the exponential in G to first order 
in a a-), that 

Soty + tySl = 0, (A75) 

This exactly corresponds to the definition of the Lie al- 
gebra Sp(2m). Having a unitary, i.e. non-singular, yet 
also antisymmetric E y , this requires a global Hilbert 
space of even dimension N, since det(E a ) = det(E^) = 
(— l) Jv det(E y ) 7^ 0. While this argument holds on the en- 
tire Hilbert space, for a specific symmetry subspace (car- 
rier space) of an irreducible representation of Sp(2m) this 
is not necessarily the case. Specifically, there are IREPs 
with odd dimensions, a simple example being the scalar 
representation with dimension 1. Within such an irre- 
ducible representation, a non-singular antisymmetric E y 
does not exist. This is not a problem, however, since 
the existence of E y is required only globally, and also in 
the defining representation, which thus has to be of even 
dimension. 

Consider such a matrix representation of Sp(2m) of 
even dimension, which allows to explicitly construct the 
non-singular antisymmetric E y . In this case, an arbitrary 
matrix S( a ) within the space of the generators of the 
symmetry can be written as a tensor-product with a two- 
dimensional space, which itself can be expanded in terms 
of the Pauli matrices r CT [cf. Eq. (A6)], 

3 

S^E^r^sW, (A76) 

:e=0 

where x € {0, 1, 2, 3} = {0, x, y, z} and t = I*- 2 ) the two- 
dimensional identity matrix. Here the same letter S is 
used left and right in Eq. (A76), as their interpretation is 
related. Nevertheless, they refer to different objects. So 
in order to distinguish them, the generators on the l.h.s. 
are written with Greek-letter subscripts (a), while their 
decomposition is denoted in roman font with roman 
or numeric subscripts. Moreover, for readability, the in- 
dex a referring to a specific generator will be skipped in 
the following where not explicitly required (hence the a 
has been put in brackets). 
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Now, with representations of a symmetry unique up 
to similarity transformation, one is free to choose the 
form of the matrix representation of the operator T, y 
in Eq. (A75). In the two-dimensional (block) space 
described by the Pauli matrices then, Yi y is chosen as 
follows, 38 - 39 



1 



(rn) _ 



0^ r 



(A77) 



with (Z 



(m)\SU(m) 



the standard to x m dimensional z- 



operators for SU(m). By construction, all of these z- 
operators can be considered diagonal, as they form a 
mutually commuting set of matrices. 

Leaving the space of strictly hermitian generators, the 
canonical RLOs from the SU(m) subalgebra are given by 



Sij — 







(i * 5) 



(A81) 



where the last term explicitly denotes the tensor block- 
decomposition of to x to matrices, with O^™ 1 ) \l^ m '] an tox 
to dimensional zero [identity] matrix, respectively. This 
Ti y fulfills the minimal requirement that it is (i) unitary 
and (ii) antisymmetric. Using the E y in Eq. (A77) in 
the defining equation for Sp(2m), Eq. (A75), and the 
fact that the generators Si a \ in Eq. (A76) shall refer to 
hermitian operators to start with, this implies for the 
decomposition S x , that So = zA is a purely imaginary and 
antisymmetric matrix, while the remaining S x for x — 
(1,2,3) must be real symmetric matrices. In summary, 
this allows to rewrite the matrix block-decomposition in 
Eq. (A76) in the form, 38 ' 39 



S = 




S 3 Si - iS 2 
iS 2 iA - S 3 , 




where C = iA + S3 (D = Si + iS 2 ) is an arbitrary hermi- 
tian (symmetric) to x to matrix, respectively. The result- 
ing number of free parameters is m 2 + m for the matrix D 
(where the +m comes from the fact that the diagonal can 
be fully complex), and to 2 for the hermitian matrix C. 
The total number of free parameters of the (N = 2m)- 
dimensional matrices therefore is, 



with Sij = Eij given by the non-symmetric matrices in 
Eq. (A27b). This encodes both, raising and lowering op- 
erators, depending on i < j or i > j, respectively. Having 



(to 2 - 1) + 1 



generators from the SU(to) subalge- 



bra together with Z m , the remaining to (to + 1) operators 
are split equally into complimentary raising and lower- 
ing operators. The corresponding canonical RLOs can 
be chosen as follows, 38 ' 39 



(all 



(A82) 



with the symmetric matrices Sij = \(Eij + E 3 i). Here 
the tilde serves to differentiate the RLOs from the SU(to) 

i.e. 



(A78) subalgebra in Eq. (A81). Having symmetric Si 



Sij = Sji, this describes a total of |m(m+ 1) raising 
operators. Complemented by |m (m + 1) lowering oper- 
ators, indeed, this completes the group of generators for 
the Lie algebra Sp(2m). 

Using the canonical representation for SU(to) together 
with above extension to Sp(2m), this provides the canon- 
ical representation for Sp(2to) as in Eq. (A20). For ex- 
ample, with 



g = to (2m + 1) = f (N + 1) 



(A79) 



In case of the defining representation, by construction, 
this also corresponds to the dimension of the symmetry 
group Sp(2m). For comparison, for example, the orthog- 
onal group O(N) has dimension f (N - 1). 

Setting the off-diagonal block-matrix D in Eq. (A78) to 
zero, and using arbitrary hermitian yet also traceless ma- 
trices C, this directly demonstrates that SU(m) is con- 
tained as a subalgebra within Sp(2m). This subalgebra 
SU(m) has rank to — 1, i.e. has to — 1 z-operators. Now, 
the full Sp(2m) symmetry also includes the tracefull her- 
mitian matrix C. This introduces the remaining TO-th z- 
operator, Z m = t z ® l^ m K With a total of to z-operators, 
Sp(2m) therefore has rank to, with the z-operators given 
by 



,(m) 



where 



1 m 



k = 



(A80a) 



(A80b) 



(A83) 



referring to the i-th diagonal matrix element of the diag- 



onal matrices Zj^ , it follows 



7W5 



5 y( m ) 
D ij /J k 



(Jm) , lm)\ a.. 
\ z k.i ^ z k,j ) °m 



j 



and thus 



[Z k , fig] = [t z (8 Z k m) , \{t x ± iT v ) ® Sij] 

= \[t Zi t x ± iT y ] ® {(zi + Zj) k Sij) 
v v ' 

= ±(r x ±ir„) 

= ±(z i + z j ) k ■ Sfj, 



(A84a) 



(no summation over i or j). Similarly, for the RLOs Sij 
from the SU(to) subalgebra, with 



(rn) q I 

k ' 0l iJ 



{ z i z j)k ' Si 
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it follows, 



[Z k ,Sij] - 







-s? 



{ Z i Z j)k ' Si; 



(A84b) 



since [Z, S T ] = — [Z,S] T . This confirms that the z- 
operators together with the raising and lowering oper- 
ators are in the expected canonical form, indeed. 



a. Internal multiplet ordering 

The block-decomposition of Eq. (A76) is not yet or- 
dered w.r.t. to the RLOs, i.e. the z-labels [here, by def- 
inition, it is assumed that a raising (lowering) operator 
leads to a larger (smaller) z-label in root space which di- 
rectly links to the underlying sorting implemented in root 
space]. The starting point, however, is correct: (i) The 
(D = 2m) dimensional first state |ei) [cf. Eq. (A27a)] 
does represent the maximum weight state, indeed, and 
(ii) by applying the m — 1 lowering operators from the 
SU(m) subalgebra, this iteratively demotes the MW- 
state through the states |e2), . . . , \e m ). So far the state 
order is correct. 

However, the next lower state is obtained by the m- 
th lowering operator, i.e. the one that links to the full 
Sp(2m) symmetry. This will generate the state \erj), 
which thus is not in order. Through another sequence 
of lowering operators from the SU(m) subalgebra, finally 
this proceeds through the states |eD_i), \eD-m+i) 
with additional alternating signs. The full sequence of 
normalized states thus obtained starting from the MW- 
state, can be collected as columns into a unitary matrix 
U, 




u 



with the to x to dimensional matrix £( m ) 



(A85a) 



^(m) 



/ 



V+ 1 



(A85b) 



to be distinguished from in Eq. (A77) associated with 
time-reversal symmetry. The unitary U in Eq. (A85a) 
maps the basis into the correct order w.r.t. to sorted z- 
labels, as is assumed throughout this paper. Therefore 
this basis convention will be used henceforth, which re- 
quires U to be applied to all generators. 



The transformation of an arbitrary symmetry opera- 
tion S in Eq. (A78) then leads to S -> WSU, that is 




c (s^c) 1 



(A86) 



In E T D, E T flips the order of the rows in D with alter- 
nating signs, starting with +1 on the new first row. The 
transformation C 4 = S T C T E in the lower right block, fi- 
nally, corresponds to inversion of C w.r.t. its center with 
alternating checker-board like minus signs applied, start- 
ing with plus signs along the regular matrix diagonal. 
With C hermitian, when taken real, C is equivalent to 
transposition w.r.t. the minor diagonal, 39 thus indicated 
by superscript lowercase t [this is in contrast to the stan- 
dard transposition (-) T around the regular diagonal]. 

All generators inherited from the SU(to) subalgebra 
thus become 



(A87) 



In particular, all z-operators have the diagonal in the 
lower-right diagonal flipped to reverse order. The sim- 
ple RLOs from the SU(m) subalgebra now have two 
strictly positive entries +1 at the first upper subdiag- 
onal at symmetric positions w.r.t. the center of the ma- 
trix. The remaining simple raising operator completing 
the Sp(2m) algebra (see below) is given by the matrix 
S mm = E mm -> E mm E = +E ml in the upper right 
block, thus naturally completing the set of simple raising 
operators of the type 
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(A88) 



with a = 1, . . . , to indicating the position of the entries 
of 1 moving towards the center of the first upper off- 
diagonal. 



b. Multiplet labels for Sp(2m) for m = 3 

With the RLOs defined to have at most two matrix ele- 
ments exactly equal to 1, the canonical commutator rela- 
tions in Eqs. (A84) directly depict the diagonal elements 
of the z- matrices. As already indicated in Eq. (A83), 
these diagonals can be combined as rows into an r x D 
matrix Zk,i, to be referred to as z-matrix, with r = m 
being the rank of the symmetry and D = 2m the dimen- 
sion of the defining matrix representation. The vectors 
Zi in Eq. (A83) thus refer to the columns in the z-matrix, 
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and therefore directly reflect the q 2 -labels, i.e. the root 
space. 

For Sp(6), this 3x6 dimensional z-matrix reads 



(\ 
1 

V 1 




(A89) 



By construction, all matrix elements are integers, for sim- 
plicity. The z-labels of the defining representation are di- 
rectly specified by the columns z\ of the z-matrix. More- 
over, since the z-labels are additive for tensor-product 
spaces, this implies that the z-labels for arbitrary IREPs 
also contain integers only. 

Consequently, the root space is fully spanned by sim- 
ple linear integer combinations of the vectors z\. Fur- 
thermore, also the action of the RLOs themselves can be 
expressed as simple shifts in root space [cf. Eqs. (A84)]. 
While in the defining representation, the z-labels in the 
carrier space are clearly unique, they are not linearly in- 
dependent. In particular, it is sufficient to focus the dis- 
cussion on the linearly independent subset of the vectors 
z*i in terms of the leading 3x3 block z of the z-matrix 
in Eq. (A89). 

In terms of the three column vectors z\ in z, the simple 
roots are given (i) by the simple roots of SU(m), which 
(ii) is complemented by one further root involving z 3 , 
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where the correspondence with the raising operators indi- 
cated in the last column follows from Eq. ( A84) . Having 
di ■ dj < for i 7^ j together with taking smallest in- 
teger combinations derived from the action of RLOs in 
Eq. (A84), this suggests simple roots. 38 ' 39 

Similar to SU(A^), the convention on the sorting of 
the z-labels is chosen lexicographic, yet as always, start- 
ing from the last z-label. In this sense, the vectors a, 
in Eq. (A90) are greater than (0,0, 0) T , hence positive. 
The corresponding operators thus increase the z-labels, 
i.e. correspond to raising operators, indeed. Moreover, 
having reduced the symmetry to its simple roots, equiv- 
alcntly, this also defines the set of simple RLOs that are 
sufficient to fully explore multiplet spaces. Note that 
above convention on the sorting of the z-labels is already 
also consistent with the state order in the defining rep- 
resentation in Eq. (A89): the z-labels strictly decrease, 
starting from the MW-state (the very left column) all 
the way to the last state represented by the very right 
column. 

In principle, the z-labels of the MW-state already could 
be used as labels for the entire multiplet. However, us- 
ing the vectors Zi as (non-)orthogonal basis that spans 
the root space, also q = z -1 max {if} could be used as 
multiplet label, instead. The latter has the advantage 



that it guarantees that the multiplet labels are strictly 
positive integers or zero. For consistency with literature, 
however, the multiplet labels for Sp(2m) are still modi- 
fied somewhat further, and thus finally derived from the 
MW-state as follows, 



q — Alz \ ■ max {z} 

=Q 

where the matrix M . 



M 



(A91a) 




(A91b) 



has been added as a further minor modification for con- 
sistency with standard literature 41 which further en- 
sures that the multiplet labels lie dense, i.e. with q = 
((71,(72,93) any qi > will result in a valid multiplet. 
Overall, 



Q 



(A91c) 



For example, when applied from the right to the z-matrix 
in Eq. (A89), all resulting matrix elements (z-labels) are 
either ±1 or 0. In particular, the MW-state of the defin- 
ing representation of Sp(2 • 3) has the q-labels (1, 0, 0). 



c. Construction of Sp(2m) for m-channel setup 

Given the three-channel setup in the previous section 
with m = 3, the resulting defining representation for 
Sp(2m) is (2m = 6)-dimensional. As seen from the ear- 
lier introduction of this model in Sec. A 9 d, this contains 
an SU(3) subalgebra, together with a third z-operator, 
namely total particle conservation. This subalgebra of a 
total of 9 generators can now be completed by 6 raising 
operators together with their hermitian conjugates, i.e. 
their corresponding lowering operators. This leads to a 
total of 21 generators, consistent with the dimension of 
the group Sp(2 • 3). 

Using a sorted z-label space, this requires that the 
unitary U in Eq. (A85a) is applied to all generators 
of the defining representation, as well as to the initial 
spinor = (c it , . . . , c mT , c\ v ... , c ] mi ) T derived from 

Eq. (A76). In case of m = 3, the properly sorted 6- 
dimensional spinor (IROP) spinor becomes, 



/ 
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(A92) 



43 



This naturally generalizes particle-hole symmetry in the 
presence of channel symmetry. The symmetry preserv- 
ing hopping term in Eq. (A67), for example, can now be 

written as scalar contraction hk,k+i = J^ai^klY '^k+i v 
Note that if, in addition, also SU(2) sp i n is present, this 
would further double the dimension of the IROP in 
Eq. (A92) to a set of 12 operators, such that the hop- 
ping term in Eq. (A67) can be written as single scalar 

contraction h k ,k+i = (^i 12 ') 1 ^ • V'i+i- 

All generators are given in second quantization by the 
quadratic form S a = ifi S a ij}, with S a a 2m-dimensional 
generator from the defining representation. Specifically, 
the remaining ^m(m + 1) raising operators for the m- 
channel setup in Eq. (A82) that complete Sp(2m) are 
given by 

S± = \{r x + l r y ) ® (SyE) = f ° S f^j , (A93) 

which leads to 

S+ ee = |(cl t c^ + c] f clj. (all (A94) 

This generates a pair of particles, the nature of which 
originates from the underlying general particle-hole sym- 
metry. With {^j, = Sij for v = 1, . . . , 2m, and there- 
fore 

= ft[S,,Mfa (A95) 

the commutator relations within the matrix representa- 
tions of the defining representation earlier directly carry 
over to the quadratic second-quantized operators as in 
Eq. (A94). 

Appendix B: Numerical implementation 

Tensor-product spaces are an essential ingredient to 
numerical renormalization group techniques such as NRG 
or DMRG. State spaces are enlarged iteratively by adding 
a small local state space at a time, i. e. a physical site with 
a few degrees of freedom. With respect to the descrip- 
tion of strongly-correlated entangled quantum many- 
body states, this leads to a description which is well- 
known as matrix product states (MPS). Both, the ex- 
isting state space (iteratively constructed itself) as well 
as the newly added state-space, have finite dimension 
and well-defined symmetry labels. New representations 
can therefore only emerge through the tensor product 
of the two spaces. In particular, all iteratively con- 
structed quantum many body states strictly derive from 
the IREPs of the elementary sites. With operators usu- 
ally acting locally, these are also expressed in the sym- 
metries of the local basis. Furthermore, the local state 
space of a site is usually small. For example, a fermionic 
site has the four states described in Eq. (A62). Therefore 



the IREPs present within the local state space are usu- 
ally just the smallest non-trivial IREPs, often just the 
defining representation itself. For identical sites, the lo- 
cal symmetry space can be setup once and for all at the 
beginning of the calculation. 

Having identified and labeled all symmetries on the 
local site level, this sets the stage for generic iterative 
algorithms such as NRG or DMRG. The remainder is 
a large exercise on tensor-product spaces. By construc- 
tion, the iteratively combined spaces are finite, yet as 
they grow rapidly, they arc eventually truncated on the 
multiplct level while leaving the symmetry content of the 
individual multiplets, i.e. the CGC spaces, fully intact. 

1. Tensor product decomposition of symmetry 
spaces 

The decomposition of the tensor-product space of two 
IREPs into irreducible multiplets has already been dis- 
cussed more generally in Sees. A 5 and A 6. In the ac- 
tual numerical implementation, however, in particular 
the presence of inner and out multiplicity must be taken 
care of meticulously for overall consistency. This will be 
discussed in the following. 

Similar to Sec. A 6, consider a specific arbitrary non- 
abelian symmetry group S whose Clebsch-Gordan coef- 
ficients may not necessarily be easily accessible analyti- 
cally for arbitrary multiplets. Assume two of its IREPs, 
qi and 52, with dimensions d qi and d q2 , respectively, are 
known together with their irreducible representations of 
the generators I„ and I^ 2 \ specifically the z-operators 
(Cartan subalgebra) and the simple RLOs (simple roots). 
In practice, these representations either refer to small 
IREPs such as the defining representation, or have been 
generated through prior iterative calculations. As in 
Eq. (A36), consider their tensor-product, 

R tot = jfoj g 1 [< Z2 ] + g / [ ?2 ] ! ( B1 ) 

resulting in matrices of dimension D — d qi d q2 . Clearly 
the commutator relations are preserved, and the z-labels 
are additive under this operation [cf. Sec. A 6]. 

In order to determine the decomposition into IREPs, 
a tempting route may be through the construction of the 
group's Casimir operators in the combined state space 
and their simultaneous diagonalization together with the 
z-operators. However, in the presence of outer or inner 
multiplicity, subspaces exist that are fully degenerate in 
Casimir operators as well as in the z-operators. In this 
case, for overall consistency a unique deterministic algo- 
rithm must be constructed that (i) separates multiplets 
in the presence of outer multiplicity, and (ii) fixes a choice 
of basis for degenerate spaces within a multiplet in the 
presence of inner multiplicity. Moreover, the explicit con- 
struction of the Casimir operators bears some efforts of 
its own. In practice, therefore a more straightforward 
approach has been adopted, instead, as will be explained 
in the following. 
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D dimensional carrier space 
with matrix representation {R a } 



|MW>i 

71= 1 



\e)i 



|MW)„ 



generate multiplet using {RJ; 



get z-labels, sort states; 
get q-label from MW state; 
generate matrix repre- 
sentation |W from {R ] 



get |MW> n+ 

n — > n + 1 




yes 



Figure 11. Schematic procedure of state space decomposi- 
tion of given D dimensional vector space with known set of 
generators {-R^ 01 '}. 



The main hurdles in the decomposition of the tensor- 
product in Eq. (Bl) into IREPs is the possible occurrence 
of outer and inner multiplicity. The strategy employed 
here to deal with this situation is based on the uniqueness 
and accessibility of the MW (maximum- weight) states as 
introduced in Sec. A 3 c. For this, throughout the pro- 
cedure below, the same lexicographic sorting scheme of 
the z-labels, used to obtain the MW-state in Eq. (A26), 
is employed to order all states within an IREP. The sort- 
ing is descending, such that the MW-state appears first 
within a multiplet. 

Since the z-labels are additive, it also follows for a 
tensor product of two such representations that the first 
state automatically also represents a MW-state of some 
multiplet, 



|MW)i = |ex), 



(B2) 



where the vectors \eu) [cf. Eq. (A27a)] form the cartesian 
basis for the D-dimensional space of the representation 
R 1 " 1 in Eq. (Bl). Given that the MW-state of a represen- 
tation is guaranteed to be unique, 25 ' 38,39 the state |MW)i 
is already a proper symmetry eigenstate, i.e. an eigen- 
state of all z-operators. This was always double checked, 
in practice, as a safety measure. The further procedure 
then is schematically depicted in the work flow diagram 
in Fig. 11: starting with |MW) n =i, 

1. the symmetry eigenstate |MW)„ is used as the seed 
state to sequence its complete IREP (the current 
multiplet). This is done by repeatedly applying 



an arbitrary but fixed order of simple lowering op- 
erators only to the current set of vectors in the 
multiplet. Therefore starting with the MW-state 
|MW)„ and adding the newly acquired symmetry 
states one at a time, this introduces a well-defined 
state order, independent of whether their z-labels 
are degenerate or not. In the presence of inner mul- 
tiplicity, it is important to notice, however, that 
it is not guaranteed that a newly acquired state 
is automatically orthogonal to the already existing 
states within the current multiplet. Therefore, a 
newly acquired state, if it represents a new vector 
space component, must be orthonormalized with 
respect to the existing states. This is repeated, un- 
til the current multiplet space is exhausted. 

2. The states in the multiplet thus generated, by con- 
struction, already have well-defined z-labels (this 
again was double-checked, in practice); the states 
are sorted with respect to these labels in descending 
lexicographic order while keeping subspaces that 
are degenerate in the z-labels in their original or- 
der in order to remain deterministic. Within this 
order, the first state defines the label for the gen- 
erated multiplet, i.e. q — (?mw- In addition, the 
matrix representation in Eq. (Bl), when cast into 
the current IREP space results in the newly gener- 
ated irreducible matrix representation 1$ . 

3. If the D-dimensional vector space is not fully ex- 
hausted yet, a new seed state is determined by find- 
ing the smallest k for which |e^) exhibits a new vec- 
tor component w.r.t. the symmetry states already 
collected. Having started with k = 1 above, it fol- 
lows k > 1. After proper orthonormalization with 
respect to the previously explored space, this state 
becomes the next seed state. If it already does rep- 
resent a MW-state, which is typically the case in 
that it is destroyed by all raising operators, then 
|MW)„ + i has been found. Otherwise repeatedly 
apply simple raising operators on the current seed 
state until the unique new maximum weight state 
|MW)„ + i is reached. Continue with (1), setting 
n — > n + 1. 

4. If on the other hand, the D-dimensional vector 
space in Eq. (Bl) is already fully exhausted, the 
decomposition of the tensor-product space into n 
irreducible representations is completed, and the 
procedure terminates. 

Note that no explicit reference to z-labels has been 
made, except for step (2). That step, however, is actually 
not required right away for the decomposition, with its 
results only relevant for subsequent calculations. By con- 
struction, therefore this procedure is deterministic and 
does not dependent on dealing with degeneracies in the 
z-labels or inner and outer multiplicities. The MW-states 
are accessible by keeping IREPs sorted in their z-labels 
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throughout. They represent the entry point in sequenc- 
ing its IREP, which guarantees that inner and outer mul- 
tiplicities are dealt with in a consistent fashion. Finally, 
note that the choice of the seed states |MW)„, i.e. start- 
ing with +|e)fe, also provides the sign convention. 

The resulting unitary transformation into the ir- 
reducible symmetry subspaces directly determines (i) 
the Clebsch Gordan coefficients, and (ii) the matrix- 
representations of the newly generated IREPs. With only 
a few Clebsch- Gordan coefficients usually unequal zero 
and of order 1, small numbers below a numerical noise 
threshold for double-precision (10~ 12 ) are neglected, i.e. 
set to zero. Moreover, a non-zero Clebsch- Gordan coeffi- 
cient can typically be expressed as a rational number, or 
the square root of a rational number, an efficient approxi- 
mation of which can be found through continued fraction 
techniques. Therefore if an excellent fractional approx- 
imation was found within the same accuracy of 10~ , 
this rational approximation also was used, instead. 



2. State space initialization and operator 
compact ificat ion 

In the presence of several symmetries, a given state 
space is represented by a certain set of multiplet combina- 
tions. For a single fermionic site in the presence of spin- 
symmetry and particle-hole symmetry, this still can be 
easily characterized by hand [cf. Eq. (A62)]. The situa- 
tion, however, can quickly become more involved. For ex- 
ample, for a spinful three-channel calculation with SU(3) 
channel symmetry as in Eq. (25a), a site is represented by 
4 3 = 64 states (4 fermionic states for each of the 3 chan- 
nels). If for example, particle number, spin symmetry 
and channel symmetry is preserved, then this system ex- 
hibits SU(2) spin (gi U(l) ch arge ® SU (3) channel symmetry, as 
discussed in Sec. A 9 d. Given these symmetries, the 64- 
dimensional Hilbert space of a site cannot be decomposed 
into a tensor product of convenient smaller units with al- 
ready well-defined SU(2) spin ® U(l) ch ar g e ® SU(3) cha nnci 
symmetry labels themselves. For the channel symme- 
try it is essential, of course, that all three channels are 
present, while it is also essential for the spin symme- 
try that both spin species are present. Therefore in the 
example above, the 64-dimensional space of site already 
appears as the smallest building block. It can be reduced 
to a set of irreducible multiplet combinations, of course, 
but explicit determination can quickly become tedious if 
done by hand, while the problem can be tackled com- 
pletely generally and straightforwardly on a numerical 
level. 

In order to get started numerically, a simple and natu- 
ral starting point is the Fock space representation. While 
this usually does not represent the symmetry eigenba- 
sis, of course, nevertheless all generators of the symme- 
tries present, in particular its raising, lowering, and z- 
operators, are known in second-quantized form and can 
be equally constructed in Fock space. 



The z-operators typically have a simple form. In 
particular, for 3-channel setup mentioned above, the z- 
operators are already all diagonal in the Fock space, cf. 
Eqs. (A68), (A69b), or (A70). This thus already pro- 
vides the z-labels. Next, note that the order of the states 
w.r.t. to their z-labels is important for consistency with 
the Clebsch Gordan coefficients later, which suggests us- 
ing the same lexicographic order as for the determination 
of the MW-states. Sorting the states in this order and 
applying the same prescription for state space decomposi- 
tion as explained in Sec. B 1, this suffices to fully identify 
all symmetry multiplets within the given D-dimensional 
Hilbert space. 



a. Compactifying operators using Wigner-Eckart theorem 

Irreducible operator sets can be equally constructed 
starting from the Fock space representation of a seed op- 
erator that is part of some irreducible operator set. This 
seed operator is typically known, yet can be completed 
to an IROP set, by using the RLOs in Fock space repre- 
sentation and numerically evaluating the commutators in 
Eq. (A42) (see also subsequent discussion in Sec. A 7 a). 
Using the same unitary transformation that brings the 
Fock space into the correct symmetry eigenbasis as de- 
scribed above, the IROP set is rotated into the space of 
symmetry eigenstates. With this, however, this IROP set 
is still represented in the fully expanded multiplet space, 
i.e. this space still references both multiplet labels and 
their corresponding z-labels on the same flat level. How- 
ever, through the Wigner-Eckart theorem, Eq. (A43), 

<M»|f?i 19292.) = (q\\F« \\q 2 ) ■ dZ'Lw ^ 

many of the matrix elements can be related to each other 
through Clebsch Gordan coefficients. The IROP set can 
therefore be compactified as a tensor-product of reduced 
matrix-elements \\q2) in the multiplet space times 

the CGC space C [q2 ' ] , >. 

The CGC spaces are known from a separate numerical 
calculation, e.g. they can be generated by several itera- 
tions of tensor-product decompositions starting from the 
defining representation. Therefore, the final compactifi- 
cation in Eq. (B3) of the fully expanded matrix elements 
of the IROP also serves as a major consistency check. 
The first non-zero matrix-element (qq z \Fql, \q2q2z) for the 
already known multiplet spaces (q; 91,(72) can be used 
to determine the reduced matrix element (q\\F qi \\q2), 
with its corresponding Clebsch Gordan coefficient known. 
This, however, immediately predicts the existence of a set 
of other non-zero matrix elements within the same multi- 
plet spaces (q; 91,(72)- These matrix elements must exist 
and agree within numerical noise. The matched matrix 
elements are marked and considered taken care of. If 
the same value of a matrix element occurs several times 
within the multiplets (g; (71, (72) for the same z-labels, the 
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first one that matches is taken. This check is thus not en- 
tirely unique, but nevertheless a strong one, and sufficient 
to obtain the space of reduced matrix elements. Finally, 
having identified all non-zero matrix elements, the multi- 
plet matrix element space (q\\F qi \\q2) are stored together 
with their referenced CGC space in terms of a QSpace as 
discussed in the main text. 



becomes trivial (1.) with internal multiplet dimension of 
1. In order to stress the difference between the origi- 
nal z-label which can become negative, and the SU(2) 
q-labels of multiplets which are positive, by definition, 
the q-labels are therefore written as (+C Z ,S), empha- 
sizing the origin of the q-label C z being derived from a 
z-operator. 



Appendix C: Example QSpaces 

QSpaces represent an efficient numerical description of 
tensors of arbitrary rank in the presence of arbitrary 
quantum symmetries [cf. Eq. (5)]. This includes both 
abelian and non-abelian symmetries, with the extension 
to further symmetries such as point symmetries being 
straightforward. The QSpaces are decomposed into a set 
of reduced multiplet spaces together with their respective 
CGC (Clebsch Gordan coefficient) spaces. In the fol- 
lowing several elementary examples of QSpaces are given 
as they appeared in practice. Elementary QSpaces typi- 
cally have rank-2 (such as scalar operators with identity 
CGC spaces) or rank-3 (IREPs and IROPs with reference 
to standard rank-3 CGC spaces), while combinations of 
these through subsequent algebraic operations can easily 
result in higher-rank intermediate objects. 

The notation regarding the elementary data arrays will 
be as follows. Plain matrices of dimension m x n will be 
written as a = [an, . . . ,a in ; . . . ; a ml , . . . ,a mn ], i.e. m 
rows of equal length n separated by semicolons. The 
commas within a row are considered optional. In order 
to deal with mxnxk dimensional rank-3 objects, the no- 
tation {ax, a 2 , . . . , afc} is used, which shall indicate that 
the matrices ai, . . . , au, all of the same dimension mxn, 
are concatenated along the third dimension. Trailing sin- 
gleton dimensions will be considered implicit if required, 
e.g. a scalar such as 1. can stand for an arbitrary rank-r 
object in that a number also represents a 1 x 1 x . . . x 1 
object. Identity matrices of dimension n will be denoted 
by lK 

1. Fermionic site with U(l) c har g o <S> SU(2) sp i n 
symmetry 

Consider the state space of a single fermionic site with 
the four states: empty |0), singly occupied | f) and | \), 
and double occupied | tl)- The symmetries considered 
are particle conservation U(l) c h a rgc, and full spin symme- 
try SU(2) sp i n . The z-operators are C z = \{n^ + — 1) 
and S z = \ (jif — hi), with the corresponding quantum 
labels C z for charge and S for total spin. For consistency 
with later, here the charge is treated as the reduction 
of the non-abelian particle-hole symmetry to its abelian 
part, which also reduces the set of symmetry operations 
to the z-operator C z only [hence the factor |]. Conse- 
quently, the z-label of the underlying non-abelian sym- 
metry is promoted into a q-label, while the CGC space 



a. Symmetry space and operators of one site 

The states |0), |t), \\), and [fj,) already represent the 
correct symmetry eigenstates [cf. Eq. (A62)], 



multiplet space 


dimension 


\C,;S) 


dcz xd s = d tot 


|-i;0) =|0) 


1x1 = 1 


|+|; o) =\U) 


1x1 = 1 


I 0;|>={|t),|;» 


1x2 = 2 



The matrix elements of a generic Hamiltonian in this 
basis can be written as QSpace [see definition in Eq. (5)], 





\{C Z -S) (C' Z ;S') 


\\H\\ 


CGC 


spaces 


H = < 


+ i;0 +i;0 
0;| 0;I 


h i 
h i 

+ 2-0 

h i 

0,2" 


1. 
1. 
1. 


1. 
1. 

1(2) 

> 



(C2) 



The Hamiltonian is a scalar operator, hence its rank as 
an IROP can be reduced from three to two, as it is 
the only operator in its irreducible set. Consequently, 
all CGC spaces reduce to the identity, as reflected in 
the last two columns of the QSpace (C2). Each of 
the remaining two indices explicitly refers to symmetry 
states, hence the QSpace requires the two sets of q-labels 
q = (C z ; S) and q' = (C' z ; S') referring to the first (sec- 
ond) index shown in the first (second) column, respec- 
tively. With the Hamiltonian preserving the symmetries, 
it must be block-diagonal, i.e. q = q' for all records 
in (C2). Both of the symmetry spaces (±|;0J have a 
single state only, therefore the corresponding entries in 
the multiplet space /i±i/2,o are lxl dimensional blocks, 
i.e. numbers. The last symmetry multiplet (0; ^) has 
two states owing to the SU(2) symmetry, see (CI). By 
means of the Wigner Eckart theorem, the space of re- 
duced matrix elements, ^o,i/2i is therefore again a num- 
ber while the CGC space becomes a 2-dimensional iden- 
tity matrix. Therefore the most general representation of 
a scalar operator for a single fermionic level in the pres- 
ence of U(l) c harge <8> SU(2) sp ; n symmetry is given by the 
three numbers {/i-i/2,0) ^+i/2,o> ^o.i/a} m the multiplet 
space. The remaining matrix elements are constrained 
due to symmetry. 

As an example for a non-scalar IROP, consider the 
spinor of particle creation operators ^r s = {c|,c|} that 
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Table I. Example QSpaces in the presence of U(l) c har g c ® SU(2) sp i n symmetry for a single fermionic site (tpg and V>s), and for 
the combination of the state space of two sites (A-tensor). Having the two symmetries of abelian U(l) c har g c and non-abclian 
SU(2) S pi n , the respective CGC spaces C (trivial) and S appear in the right columns. The record index v in the first column, as 
well as the explicit specification of the dimensions of the reduced multiplet space and the combined CGC spaces are just added 
for better clarity. For comparison, Tbl. II shows how the QSpaces (C3a) and (C4) are modified for the case that the abelian 
charge symmetry also becomes a non-abelian SU(2) c har gc particle-hole symmetry. 



4 = < 



record 








reduced matrix elements 


CGC spaces 




index v 


(C Z ;S) (C' Z ;S') 


(C'J 


;S") 




dimension 


C S 


dimension 


I. 


0;i 


+1. 

2 ' 


i 

2 


1. 


lxlxf 


1-, {[1;0],[0;1]} 


2x1x2 


2. 


±i;0 0;| 


+ 1. 

2 ' 


1 
2 


72 


lxlxl 


1, ^{[0 1],[-10]} 


1x2x2 



i>s = < 



record 








reduced matrix elements 


CGC spaces 




index v 


(C Z ;S) (C' Z ;S') 


(C'J 


;S") 


HlM, 


dimension 


C S 


dimension 


1. 


0;| ±i;0 


-l . 

2 ' 


i 

2 


1. 


lxlxl 


1-, {[1;0],[0;1]} 


2x1x2 


2. 


^;0 0;| 


-1 . 

2 ' 


1 

2 


-72 


lxlxl 


1., ^{[0 1],[-10]} 


1x2x2 



QSpaces ip' s = {cj, cj} and %j)s = { — &i; Cf} representing IROPs for a single spinful fermionic level [c/. Eqs. (A49)]. 
Note that the IROP xj) s is interpreted differently compared to the IROP (i/>s)S hence tp s ^ (i'sY [e.g. note the 
sign in multiplet space in the second record of tps or the reverted signs in the q-labels for (C'J; S") associated 

with the IROP in the third column; see text]. 
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record 
index v 


site 1 site 2 combined 
(C Z ;S) (C' Z ;S') (C'J;S") 


multiplet space 
T dimension 


CGC Spaces 

C S dimension 


1. 


^i;0 ^;0 -1;0 


1. lxlxl 


1. 1. lxlxl 


2. 
3. 


=1-q o-- —■- 

2 ' u u ' 2 2 ' 2 
u, 2 2 ' 2 ' 2 


{[1],[0]} 1x1x2 
{[0],[1]} 1x1x2 


1. {[10], [0 1]} 1x2x2 
1. {[1;0],[0;1]} 2x1x2 


4. 
5. 
6. 


^i;0 |;0 0;0 
0;| 0;| 0;0 
|;0 ^-;0 0;0 


{[1],[0],[0]} 1x1x3 
{[0], [1], [0]} 1x1x3 
{[0],[0],[1]} 1x1x3 


1. 1. lxlxl 

L [° 71' 75 °1 2x2x1 
1. 1. lxlxl 


7. 


0;| 0;| 0;1 


1. lxlxl 


{[1 0;0 0], 

L [° 75 - 75 °]' 2x2x3 
[0 0;0 1]} 


8. 
9. 


0^ i-i 

u, j 2 ' u 2 ' 2 

1-0 O' 1 i-i 

2 , u w, j 2 » 2 


{[1],[0]} 1x1x2 
{[0],[1]} 1x1x2 


1. {[1;0],[0;1]} 2x1x2 
1. {[10], [0 1]} 1x2x2 


10. 


|;0 |;0 +1;0 


1. lxlxl 


1. 1. lxlxl 



QSpace of identity A-tensor combining two fermionic sites. Site 1 with symmetries (C z ; S) and site 2 with symme- 
tries (C' z ; S') are combined into the global symmetry (C'J; S"). The specific order of the records is irrelevant and 
hence arbitrary. Here, the records have been sorted with respect to the combined quantum labels q" = (C'J; S"), 
where groups with the same q" are indicated by horizontal lines for clarity. The dimensions in the last (third) 
index are therefore the same within a group that shares the same (C'J; S"). 



encodes SU(2) spin symmetry [cf. Eq. (A49) ], with its 
QSpace representation shown in (C3a). The z-labels of 
the IROP set tpg are determined through the z-operators 



C z and S z acting on the components of ip s , 
[C'z, c\] = • c CT 

with a = = {+1, -1}. The IROP $\ is therefore 
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identified with the multiplet q" = (C'J;S") = (^5), 
as indicated in the third column of (C3a). The QSpace 
representation of ipg derives from the matrix-elements 

4^(C Z S\-((^ S )(^^\C' Z S')) 

using the Wigner Eckart theorem as in Eq. (A43). 

The operator index in the QSpace (C3a) is listed third, 
by convention. The two non-zero matrix elements of each 
c\ within the four-dimensional space of single fermionic 
site implies a total of four non-zero matrix elements in 
tpg, all having norm 1, with one matrix-element being 
negative. These matrix elements can be directly identi- 
fied in QSpace (C3a). Since the reduced matrix elements 
\\ip s \\ and the CGC spaces are to be interpreted as ten- 
sor product, the \Jl factors in the last line cancel. With 
tyg representing non-hermitian operators, the first col- 
umn q = (C Z ;S) is in general different form the second 
column q' = (C' Z ;S'). Moreover, since Tp\ creates one 
particle, the first column, for example, cannot contain 
the empty state (^', 0), while the second column cannot 
contain the double occupied state (±^;0). 

In contrast, the QSpace representation of the IROP 
-0s, i.e. without the dagger, is shown in (C3b). Note 
that for ips = {~ cy,tf) to be an irreducible operator 
as compared to ipg — {cLcj^}, the reverse order in spin 
and the minus sign in the first component is essential 
[see discussion along with Eqs. (A49)]. In terms of the 
QSpace (C3b), this leads to the extra minus signs in the 
multiplet space of the second row. Moreover, the z-labels 
of the operator ips itself flipped sign w.r.t. ipg as expected 
as it removes a particle rather than adding one [see the 
multiplet labels q" = (C Z ;S") in the third column of 
(C3b)]. This is to emphasize that the application of the 
Wigner Eckart theorem must be performed consistently, 
i.e. switching sides in the application of an operator as 
in (C Z S\ ■ {^\C' Z S')) = C0|C z S))t • \C' Z S') must be dealt 
with carefully. 



b. Identity A-tensor for two fermionic sites 

Consider the combination of two fermionic sites. Al- 
luding to Fig. 1, let site 1 (2) be described by \i) (|cr)), re- 
spectively, both representing a 4-dimensional state space 
{|0) , |t) , ||) , It-i)} OI their own. The decomposition of 
the combined space in terms of the overall symmetry 
U(l)char g e <S> SU(2) sp i n is fully described by the rank-3 
QSpace (C4). 

Given the U(l) c har g c ® SU(2) sp i n symmetries, the 
abelian charge quantum number C z simply adds up, 
while for the SU(2) spin symmetry, the usual SU(2) ad- 
dition algebra applies. The overall number of multiplets 
in the combined space q" is given by the last number 
(index 3) in the dimensions specified with the multiplet 
space. The specific input combinations entering a certain 



combined space q" are easily verified. The q" = (— i; ^) 
sector, for example, derives from the two configurations 
{q,q'} = {(-|;0),(0;i)} and {(0; |),(-|;0)}. There- 
fore the dimension of the reduced multiplet space for 
this q" is 2. Each of these multiplets has an internal 
z-space which is itself of dimension 2 [last column] . The 
combined total dimension of the q" = (~; ~) sector is 
therefor given by the product 2-2 = 4. Consistently, 
the dimension of the two 4-dimensional sites combined 
add up correctly to 16 states total. That is, multiplying 
the last dimension in the reduced multiplet space with 
the last dimension in the combined CGC spaces for each 
block separated by horizontal lines, bearing in mind that 
the multiplet space and the CGC spaces are to be com- 
bined in a tensor-product, yields the overall dimension of 
the combined space, M + 2-2 + 3-1 + 1-3 + 2-2 + M = 16. 

The A-tensor in (C4) is an identity A-tensor, in that 
up to permutations, plain identity matrices are split-up 
on the reduced multiplet level. By considering, for exam- 
ple, the q" = (±|; \) symmetry sector in records 2 — 3 
or 8 — 9 of the QSpace (C4), the multiplet space when 
viewed together, i.e. ignoring all brackets, resemble the 
structure of a 2-dimensional identity matrix. Similar so 
for the q" = (0; 0) space in records 4 — 6, having essen- 
tially a 3-dimensional identity matrix in the multiplet 
space. Allowing for arbitrary unitaries in the multiplet 
space in QSpace (C4), this then becomes the most gen- 
eral unitary transformation of the product space of two 
fermionic sites that also respects the symmetries consid- 
ered. 



2. Fermionic sites in the presence of particle-hole 
symmetry 

The tensors introduced in the previous section for 
U(l) c hargo ® SU(2) sp i n symmetry will now be written 
more compactly still by assuming the stronger particle- 
hole SU(2) c har g c symmetry instead of the plain abelian 
U(l) c hargo- The symmetry considered in the following 
is therefore SU(2)f g = SU(2) spin ® SU(2) chargc . The z- 

operator for charge, C z , is now complemented by the rais- 
ing operator C+ for charge SU(2) [cf. Eq. (A52)]. The 
combined symmetries are given by the multiplet label for 
both charge and spin SU(2), i.e. the non-negative labels 
q = (C, S) with the z-labels of the charge symmetry now 
also taken care of by the CGC spaces. 

The basis for a single fermionic level given given 
SU(2)gQ symmetry has been introduced in Eq. (A62). 
Therefore the full space of the four states {0, t, -Utl} 
can be reduced to the two symmetry multiplets 



multiplet space 


dimension 


\C,S) 


dc x d s = d tot 


||;0> ={|0>,»m>} 


2x1 = 2 


|0;|> ={|t),|l)} 


1x2 = 2 



(C7) 



The most general scalar operator such as the Hamiltonian 
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Table II. Example QSpaces in the presence of SU(2)gQ = SU(2) sp i n <g> SU(2) c h ar go symmetry for a single fermionic site (ipcs), 
and for the combination of the state space of two sites (^-tensor). The CGC spaces for SU(2) c h arge and SU(2) sp i n are indicated 
by the C and S, respectively. The record index v as well as the explicit specification of the dimensions are just added for clarity. 
For comparison, Tbl. I shows the same QSpaces for the case where the particle- hole symmetry is reduced to abelian charge 
conservation. 





record 
index v 


(C;S) (C';S') 
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red. matrix elements 
HV'csll dimension 


CGC spaces 
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dimension 


V>cs = < 
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0;i |;0 
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2 ' 2 


sy/2 


lxlxl 


^{[o i];[-io]} 


{[1;0];[0;1]} 


2x2x4 




2. 


§;0 0;| 


1 . 1 

2 ' 2 


V2 


lxlxl 


{[1;0];[0;1]} 


7s{[o i];[-io]} 


2x2x4 



(C5) 



QSpace of spinor i/+js defined in Eq. (A59) with the reduced matrix elements for given symmetries already calculated 
in Eq. (A63). The operator index within the IREP is listed third, as usual. The alternating sign s required with 
particle-hole symmetry appears with the reduced matrix elements in the first record only [cf. Eq. (A63); note that 
the same sign s is also picked up by the double occupied state, cf. Eq. (A62) or (C7)]. 



A= { 



record 
index v 


site 1 site 2 combined 
(C;S) (C";5") (C";S") 


multiplet space 
|| A || dimension 


CGC spaces combined 
C S dimension 


1. 
2. 


0;| 0;| 0;0 
i;0 §;0 0;0 


{[1],[0]} 1x1x2 
{[0],[1]} 1x1x2 


1. ^j[0 - 1; 1] 2 x 2 x 1 
4g [0 - 1; 1] 1. 2x2x1 


3. 


0;§ 0;| 0;1 


1. lxlxl 


{[1 0; 0], 
1. ^[0 1; 1], 2x2x3 ( 
[0 0; 1]} 


4. 
5. 


0- i ! -0 i-A 

1- 0-- -■- 

2 , u u, j 2 ' 2 


{[1],[0]} 1x1x2 
{[0],[1]} 1x1x2 


{[10], [0 1]} {[1;0],[0;1]} 2x2x4 
{[1;0],[0;1]} {[10], [0 1]} 2x2x4 


6. 


|;0 |;0 1;0 


1. lxlxl 


{[1 0; 0], 

^=[0 1; 1], 1. 2x2x3 
[0 0; 1]} 



QSpace of identity ^-tensor combining two fermionic sites - the symmetry records are sorted with respect the 
combined symmetry q" = (C", S"), where groups with the same (C", S") are separated by horizontal lines for 
clarity. Considering the tensor-product of multiplet and CGC spaces, the combined space has total dimension 
of2-l + l- 3 + 2- 4+ l- 3 = 16as expected for two spinful fermionic levels. Compared to the A-tensor in 
(C4) with abelian charge conservation, the number of combined symmetry sectors has been further reduced 
from 6 to 4 [i, e. number of horizontally separated groups sharing the same q"] , with an overall reduction in 
the number of multiplets present in the QSpace reduced from 10 to 6 [having 1 + 2 + 3 + 1 + 2 + 1 = 10 in 
(C4), and here 2 + 1 + 2 + 1 = 6]. 



is given by the QSpace, 





f (C;S) (C';S') 
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CGC spaces 


-1 


|;0 i;0 
0;i 0;i 


hi 

h i 

o,i 


1. 
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Thus only the two reduced matrix elements h 1 / 2 .o and 
ho, 1/2 are l e ft f ree to choose without compromising 
SU(2)fg symmetry. 

A non-scalar IROP is given in terms of the 4-element 
spinor ^cg in Eq. (A59), which combines two creation 
and two annihilation operators. Its symmetries have been 
identified in Eq. (A61), leading to the IROP ipes with 



the QSpace presented in (C5). Note that the size of the 
third dimension is 4 [see combined CGC dimension in 
the last column] , consistent with the four operators that 
constitute the IROP. 

The QSpace representation for the A-tensor combining 
two fermionic levels is given in (C6). The standard SU(2) 
addition rules are quickly confirmed. For example, the 
combined symmetry q" — (0, 0) [records 1-2] can result 
from two combinations, namely (| , 0) (g> ( \ , 0), or (0, |) ® 
(0, 1), leading to a two dimensional multiplet space. All 
of this is transparently encoded in the QSpace. Here, 
QSpace (C6) again shows an identity A-tensor, as seen, 
for example, in the combined reduced multiplet space of 
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CG dimension 
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1. 


0; 


0; 


0; 


1 
2 


1 


1x1x1x2= 


2 


2. 


0; 


0; 


1 . 
2 ■ 





1 


1x1x2x1= 


2 


3. 


0; 


i . 

2 ' 


0; 





1 


1x2x1x1= 


2 


4. 


0; 


1 . 

2 ' 


1 . 
2 ■ 


1 
2 


1 


1x2x2x2= 


8 


5. 


||; 


0; 


0: 





2 


2x1x1x1= 


2 


6. 


||; 


0; 


1 . 

2 ■ 


1 
2 


1 


2x1x2x2= 


8 


7. 


II; 


l . 

2 ' 


0: 


1 

2 


1 


2x2x1x2= 


8 


8. 


II; 


1 . 

2 ' 


1 . 

2 ' 





1 


2x2x2x1= 


8 


9. 


i; 


0; 


0; 


1 
2 


1 


3x1x1x2= 


6 


10. 


i; 


0; 


1 . 

2 ■ 





1 


3x1x2x1= 


6 


11. 


i; 


i . 

2 ' 


0: 





1 


3x2x1x1= 


6 


12. 


If; 


0; 


0; 





1 


4x1x1x1= 
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Table III. State space of 3-channel site with SU(2)f c 
SU(2) spin ® SU(2)® 3 Mgc symmetry. 
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Table IV. State space of 3-channel site with SU(2) sp i n 
U(l)char gc ® SU(3)channoi symmetry. 
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CG dimension 




index 


\s 


Sp(6) } dim. 


ds ■ dc z ■ d T = 


dtot 
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|0; 


1) 1 


1 x 1 x 14 = 
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Table V. State space of 3-channel site with SU(2) sp i n <g>Sp(6) 
symmetry. 



records 1 — 2, or 4 — 5). 



3. Three channels with SU(3) channel symmetry 

Consider a system with three spinful particle-hole sym- 
metric channels, as introduced in Sec. A9d. A single site 
then has a full Hilbert space of dimension 4 3 = 64. Three 



symmetry settings are analyzed: a set of plain SU(2) 
symmetries, a combination with the SU(3) channel sym- 
metry, and finally the largest symmetry present which in- 
cludes the enveloping symplectic symmetry Sp(G). Since 
the QSpaces in given context are extensive, a more com- 
pact comparison of these symmetry settings in the nu- 
merical context is given, instead. 

The first setting, SU(2)fg = SU(2) spin <g> SU(2)f h 3 arge , 
is based on four independent SU(2) symmetries. The 64- 
dimensional state space of single site decomposes into 13 
multiplets in 12 symmetry sectors, as listed in Tbl. III. 
All of these contain a single representative multiplet, ex- 
cept for the space 1 1 ; 0; 0; 0} in row 5, which contains two 
multiplets. In contrast, using the SU(2) sp j n (g)U(l) c hargc® 
SU(3) c hannci symmetry, instead, the 64 states of the three 
fermionic levels decomposes into the 10 multiplet spaces 
listed in Tbl. IV. Thus compared to the SU(2)§| sym- 
metry setting in Tbl. Ill, the number of multiplet spaces 
is further reduced with all multiplet spaces containing a 
single multiplet only. This suggests that the latter sym- 
metry including the channel SU(3) is somewhat more ef- 
ficient as it allows to compactify multiplet spaces more 
strongly. Given the multiplet space in Tbl. IV, for ex- 
ample, the most general Hamiltonian in the 64 x 64 di- 
mensional Hilbert space compatible with given symmetry 
consists of the 10 reduced matrix elements appearing in 
the multiplet space only. 

A further strong boost in numerical efficiency can be 
obtained, if the Hamiltonian supports it, by combin- 
ing the particle-hole symmetry of SU(2) spin <g)SU(2)f h 3 arge 
with the channel symmetry in SU(2) spin (g) U(l) c h a r g c ® 
SU(3) c hannci to their enveloping Sp(6) symmetry. The 
resulting state space for the state space of a single 3- 
channel site is given in Tbl. V. The 64 x 64 dimensional 
Hilbert space has thus been reduced to a total of four 
multiplets only. 

All three symmetry settings have been successfully im- 
plemented within the NRG framework. By starting with 
a single site [i.e. the basic fermionic three-level unit as 
introduced in Sec. A9d], and iteratively adding a site 
within the NRG, new multiplet spaces are quickly ex- 
plored and built up within the first few NRG iterations. 
In practice, the CGC spaces of newly generated multi- 
plets are also stored for latter retrieval. Once truncation 
of the state space within NRG sets in, however, the gen- 
eration of IREPs eventually saturates to within a finite 
range of multiplet spaces. 

The resulting space for adding up to three further 
sites without truncation is indicated in Tbl. VI. With 
reasonable numerical resources, it is feasible within 
the SU(2) spin <g) U(l) ch ar g c ® SU(3) ch annoi (second) or 
SU(2) S pi n <E) Sp(6) (third) setting, to keep all states up 
to three sites total within the NRG, first truncating only 
when a fourth site is added. This leads to Nk — 9, 086 
[Nk = 1,232] kept multiplets for the second and third 
symmetry setting, respectively. The corresponding mem- 
ory requirements for a general basis transformation for 
adding another site (A-tensor) then amounts to about 
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sites 
n 


abelian dim. 
D* = 64™ 


SU(2) spin ® SU(2)f h 3 rge 
N s D D'/fl memory MEM* 


SU(2) spin ® U(l) ch argc ® SU(3) cha nncl 

N s D D* /D memory MEM* 


SU(2) spin ® Sp(6) 
N s D D*/D memory MEM* 


1 

2 
3 
4 


64 

4,096 
262, 144 
16,777,216 


12 13 4.9 <18K 

61 388 10.6 528K >8.7M 
192 14,229 18.4 27 M >23G 
469 590,856 28.4 24 G >65T 


10 10 6.4 < 13K 
69 260 15.8 359 K >12M 
226 9,086 28.9 11 M >31G 
565 366,744 45.7 6.8 G >85T 


4 4 16 <6K 
23 61 67.1 162K >34M 
60 1,232 213 7M >112G 
132 31,640 530 334 M >355T 



Table VI. Comparison of different symmetry scenarios for the same underlying physical system of a symmetric 3-channcl 
setup, analyzing the product spaces of up to n — 4 sites. Each site represents a Hilbert space of dimension 64, thus n sites 
amounts to an overall Hilbert space of dimension D* — 64" [second column]. This state space can be decomposed into D 
multiplets in Ns symmetry sectors using an ^4-tensor for the addition of every new site. These ^-tensors are encoded in terms 
of QSpaces. The total memory requirement for each ^-tensor is listed, given sparse CGC representation. In addition, as a 
comparison to a fully abelian setting, MEM* indicates the memory that had been required if the tensor products between 
reduced multiplets and CGC spaces was carried out explicitly [K,M,G,T for kilo-, mega-, giga-, and tera-bytes, respectively]. 



sites 
n 


SU(2) spin ® SU(2)f h 3 rge 
ds dc dc dc sparsity CGS/A 


SU(2) spin 
ds dc/(i) 


® U(l)chargo t 

dc sparsity 


3 SU(3)channol 
CGS/A 


SU(2) spin ® Sp(6) 
ds rfs P (6) sparsity CGS/A 


1 


4 2 2 2 


0.36 


0.8 


4 


1 


8 0.12 


0.8 


4 14 0.027 0.90 


2 


7 3 3 3 


0.36 


0.8 


7 


1 


27 0.12 


0.8 


7 126 0.027 0.90 


3 


10 4 4 4 


0.28 


0.11 


10 


1 


64 0.064 


0.33 


10 616 0.011 0.94 


4 


13 5 5 5 


0.22 


< 10" 3 


13 


1 


125 0.039 


0.003 


13 2457 0.006 0.55 



Table VII. Comparison of different symmetry scenarios as in Tbl. VI in terms of (i) largest multiplet dimension d for each 
individual symmetry, and (ii) overall average sparsity of the CGC spaces, i.e. the number of non-zero elements divided by the 
total number of matrix elements. The last columns for each symmetry (CGS/A) shows the memory requirement of all sparse 
CGC spaces in a given QSpace A n relative to the entire QSpace. 



7G [0.3 G]. The corresponding full NRG iteration for 
adding another site then takes several hours on a state- 
of-the-art 8-core workstation. 

The same calculation, however, gets quickly impossible 
as fewer symmetries are available or used in the actual 
computation as can be seen from Tbl. VI. For example, 
if only the abelian part of the symmetry had been ac- 
counted for in the computation, the corresponding mem- 
ory requirement can be estimated by considering the ex- 
plicit tensor product of the multiplet space with the CGC 
spaces, leading in terms of the SU(2)^ setting to about 
23 G and 65 T(!) for n = 3 and n = 4, respectively, 
the latter being completely hopeless in practice. The 
explicit treatment of non-abelian symmetries, however, 
clearly makes the latter case feasible with a reasonable 
amount of numerical resources. 



a. SU(3) symmetry 

The irreducible SU(3) multiplets generated in the ac- 
tual NRG run using SU(2) spin ® U(l) cha rgc ® SU(3) ch annci 
symmetry as presented in the main text [cf. Sec. V], 
are shown in terms of their weight diagrams in Fig. 12. 
For comparison, the SU(3) IREPs present in the de- 
scription of a single site are g su ( 3 ) = (91,(72) G 
{(0,0), (0,1), (0,2), (1,0), (2,0)}, cf. Tbl. IV. The ap- 
parent symmetry w.r.t. to flipping the quantum numbers 
in (31,92) is also reflected in the overall set of multiplets 



generated within the NRG. As seen in Fig. 12, all IREPs 
(qi, 92) with qi+q2 < 8 are present, except for (0, 8) and 
(8,0). 

Inner multiplicity, as expected for SU(3), is clearly 
present and depicted in the weight diagrams of Fig. 12 by 
the encircled set of points. There the number of points 
inside a circle stands for the multiplicity of the corre- 
sponding z-labels in the multiplet. Inner multiplicity 
decreases in shells as one moves outward the multiplet, 
which is seen particularly well for the multiplets q\ = 92 • 
The states on the outer circumference have no multiplic- 
ity, i.e. have unique z-labels, as expected. This demon- 
strates the uniqueness of the maximum weight state [cf. 
Sec. A 3 c], which was required for the numerical state 
space decomposition in Sec. B 1. 

Due to the two-dimensional label structure of the 
SU(3) multiplets together with inner multiplicity, their 
internal dimensions can get significantly larger as com- 
pared to SU(2) multiplets. The largest SU(3) multiplet 
(4, 4) encountered in the actual NRG run presented in 
Fig. 12, for example, has an internal irreducible dimen- 
sion of d — 125 [see also Tbl. VII]. This implies, for 
example, that with respect to the diagonalization of a 
Hamiltonian, a 125-fold degeneracy has been reduced to 
a single multiplet. In contrast, the multiplet structure 
for SU(2)®q is clearly weaker as it only includes SU(2) 
symmetries. There the largest quantum numbers encoun- 
tered in an NRG run with comparable number of kept 
states includes S < 6 in the spin sector, leading to an 



■52 



10 
J, 
-10 

10 
-10 

10 
J, 
-10 

10 

«? 

-10 

10 
-10 

10 

o 5 
-10 

10 
^ 
-10 

10 

«? 

-10 



(0.0) : 


(0,1) : 







• 




d=1 


d=3 



(0,2) 



••• 



d=6 



(0,3) 



-A- 



d=l0 



(0,4) 



d=l5 



(0,5) 



••• 
. •••• . 
••••• 



d=21 



(0,6) 



-Aft- 
••••••• 



d=28 



(0,7) 



•V 



•••• 

.*»». 
••«••••• 



d=36 



(1.0) 



d=3 



(1.1) 



■ w 



d=8 



(1,2) 



•ML 



d=15 



(1.3) 



d=24 



(1,4) 



■ *©©©•■ 
• ®®.®®. 



d=35 



(1,5) 



JOSft. 



d=48 



(1,6) 



■ ■ ■ ■•<!)• ■ ■ ■ ■ 

•3XD* 
• ®®®« 
. ..®®<1)®. . 
•®®®©®. 
m^j^j^)^ 



d=63 



(1,7) 



mm 

■ ' ■ ■ 

•©©©• 
•©©■©©• 
•©©©©©•■ 
•®ffl®.®®®. 

V.v.v.v.* 

d=80 



(2,0) 



d=6 



(2,1) 



••• 

■ »<e®«. 



d=l5 



(2,2) 



d=27 



(2,3) 



•®@®« 

■*®©®®. ' 



d=42 



(2,4) 



. . .... . . . 

•©<©• 
•©«©• 

..®.®-®®.. . 
•<»JK®J®J3>» 



d=60 



(2,5) 



••• 



•©©■©©• 

■ •<£©©©©•■ ■ 



d=81 



(2,6) ... 
.... ..©©.. . . 

•©©©• 
•©©.©©. 
.®®®@®. 
■»®®®©©©. 
.©®®®®®®. 
m^j^j^j^jm 

d=105 



10 



10 



(3.0) : 


(3,1) : 






• 


.... ■ ■ ■ ■ 

mm 


d=10 


d=24 



(3.2) : 


(3,3) : 


•••• 


mlS^m 


•©©©• 

■■■■■vmr---> 


•®®®®. 


. *»• ; 




d=42 


d=64 



(3,4) 



. . ••«.. . . 
•©©©• 

•®®.®©. 
•®@®®®. 



m^xi^m 



d=90 



(3,5). 



.©©©>©• 
.©©©®®. 

«®®.@©®©. . 



•®@®®©®. 

-WW--. 

: d=120 



10 



10 



(4,0) 



*•*•*•*•* 
••• 

V 



d=15 



(4,1) 



JMMMML 

WW 



% 



d=35 



(4,2) 



••••• 

*®®.®®. 
•®®®@®. 

...*5fr... 



d=60 



(4,3) 



•©©<©©• 
.©@®@®. 

.©©.©■©©©• 



d=90 



(4,4) 

■.®®®®. ' ' 
•®@®@®. 

•®©©'@®®. 

•©©>©©©©©• 

•®®®®®. 



d=125 



10 



10 



(5,0) 



...... 

••••• 

w 



d=21 



(5,1) 



•®®®®®« 

mjmgm 
mm>m 



mm 



d=48 



(5,2) 



■ »®©©®®.' 



d=81 



(5,3) : 

' iStSffiSSSm' ' 
.©©®.@®®. 
•®®©©©©®. 
, .®© , ©@©®. 1 
.®@@©®. 

■ «®®®« ■ ■ ■ ■ 
•••• 

d=120 



10 



10 



(6,0) 



w.y.v 



.... 

•A* 



d=28 



(6,1) 



■ ■ ■ 

d=63 



(6,2) 



m^^$<^*m 

..®®.®®®®. . 
•®®@®®. 

•«? . 
••*d=105 



10 







10 



(7,0) 



•••••••• 

••••••• 

. . •••••• . . 

••••• 

••» # 

* d=36 



(7,1) 



•••••••• 

•®®®®®®®« 
•®®®JJ>®®« 

•®®®« 



d=80 



10 



10 



-10 



10-10 



10 



Figure 12. Weight diagrams of SU(3) multiplets generated in a typical NRG run for the symmetric 3-channel system 
including SU(2) sp i n <g> U(l) c har g e ® SU(3) c hannci symmetries (A = 4, St runc = 7). The multiplet label (91,172), as defined in 
Eq. (A34), is specified with each multiplet in the upper left corner of its panel. For a weight diagram of a specific IREP 
(91,92), the corresponding z-labels (921,922) of SU(3) for each individual state within the multiplet are depicted as points in a 
two-dimensional plot. In case of inner multiplicity, i.e. that several states within the same IREP share exactly the same z-labels, 
these states are shown as an encircled group of smaller points. The dimension for every multiplet (number of points drawn 
within a panel) is indicated to the lower right of each panel. The first panel [multiplet (0, 0)] represents the scalar representation 
with multiplet dimension d = 1. Multiplet (1,0) represents the defining three-dimensional representation [cf. Sec. A4c], and 
(1, 1) the regular representation of dimension 8 equal to the dimension of SU(3), i.e. the number of its generators. The largest 
multiplet encountered in given NRG run is the multiplet (4, 4) with an irreducible dimension of d = 125. 
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Table VIII. Sp(6) multiplets generated in a fully con- 
verged NRG run for the symmetric 3-channel system using 
SU(2) apin ® Sp(6) (A = 4, £ trunc = 7). Multiplet (0,0,0) 
represents the scalar representation of dimension 1, multi- 
plet (1,0,0) the defining representation of dimension 6, and 
multiplet (2, 0, 0) the regular representation of dimension 21 
which is also equal to the number of generators for Sp(6). 
The largest tensor-product decomposition was between the 
product spaces of IREPs of dimension 14 and 512, yielding 
a combined product space dimension of 7, 168. Run time of 
the bare NRG run was about 2 hrs on a state-of-the-art 8-core 
workstation with moderate memory requirements of < 4.5 G. 

individual multiplet dimension of at most 13 [see also 
Tbl. VII]. In the overall combination of the symmetries, 
this implies that for comparable number of states, i.e. 
for a comparable accuracy within the NRG, on average 
about 50% more multiplets need to be kept within the 
SU(2)i?g setting as compared to the case when SU(3) is 
included [see Fig. 9 in the main text]. 

Finally, the individual weight diagrams in Fig. 12 show 
well-known symmetries, such as a reflection symmetry 
of each diagram around the vertical y-axes, or the re- 
flection symmetry between the multiplets (91,92) and 
(<72,<Zi) around the horizontal axis. These Weyl sym- 
metries may be used to evaluate or encode CGC spaces 
more efficiently. 25 For the purpose of this paper, however, 
these symmetries were not exploited, given also that the 
pure numerical evaluation of the CGC spaces as outlined 
earlier was already sufficiently fast. 



b. Sp(6) symmetry 

The complete set of Sp(6) symmetries generated in the 
fully converged NRG run (using A = 4 and -Etrunc = 7 as 
used in the results in the main text), is listed in Tbl. VIII. 
All multiplets had been generated within the first four 
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Wilson shells. The fact that the symmetry Sp(6) fully 
incorporates non-abelian particle-hole and channel sym- 
metry, manifests itself by observing that all eigenenergies 
in the multiplet spaces are now strictly non- degenerate 
throughout an entire NRG calculation. Huge degenera- 
cies of several thousands can be split off in terms of tensor 
products with Sp(6) multiplets. 

For given model, the symmetry Sp(6) in fact also al- 
lowed to reduce the rather coarse discretization of A = 4 
in the NRG calculation underlying Tbl. VIII. For com- 
parison, if A = 2 is used, instead, while keeping the same 
-E'trunc = 7, it turns out, the largest multiplet generated 
is (2,1,2) of dimension 5720. The largest intermediate 
product space to be decomposed into IREPs becomes as 
large as 14 x 1386 = 19, 404. Having A = 2, this required 
twice the Wilson chain length for the same range in en- 
ergy scales, leading to an overall run time of the entire 
NRG run of about 32 hrs with still reasonably manage- 
able memory requirements of < 20 G. 

As a rough general estimate, typical multiplet dimen- 
sions, as they occurred in practice, scale like 10 r where 
r is the rank of the symmetry. For SU(2) sp j n , this im- 
plies multiplets of dimension < 10, for the SU(4) c hannci 
symmetry, indeed, one had multiplets of dimension of 
< 100, while now for Sp(6), a symmetry of rank 3, one 
easily reaches multiplet dimensions on the order of a few 
1000 (cf. Tbl. VIII). Therefore with increasing rank of the 
symmetry, the numerical effort strongly shifts from the 
multiplet space to the CGC spaces. For sets of smaller 
symmetries with rank r < 2 this leads to a strong gain 
in numerical efficiency, while the numerical overhead for 
the CGC spaces remains negligible. Reaching symme- 
tries of rank 3, such as Sp(6), the numerical effort within 
the CGC spaces can now become comparable to or even 
larger than the operations on the higher multiplet level. 

Table VII summarizes the situation by comparing the 
maximal multiplet spaces with the corresponding spar- 
sity and memory requirements of the CGC spaces for the 
first few A-tensors, when combining up to n = 4 sites 
without truncation. As the internal multiplet dimensions 
quickly grow for higher rank symmetries, nevertheless 
only an ever smaller fraction of Clebsch-Gordan coeffi- 
cients are non-zero. As seen from Tbl. VII, the sparsity 
roughly grows exponentially with the rank of the symme- 
try. Nevertheless, with the memory requirement of the 
sparse CGC for Sp(6) comparable or even larger than 
the storage of the reduced matrix elements on the higher 
multiplet level (see last column in Tbl. VII), full storage 
including also the zero CGC spaces would have extremely 
inflated overall storage requirement. In this sense, here 
sparse storage of CGC spaces becomes mandatory. 

In order to distinguish numerical noise, i.e. negligible 
CGC matrix elements, from actual matrix elements then, 
this requires an accurate evaluation of the CGC matrix 
elements. Double precision accuracy as compared to the 
exact theoretical CGC coefficients was sufficient, in prac- 
tice. In particular, this implies, that also the matrix ele- 
ments of the generators for given IREPs of the symmetry 
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are known essentially exactly at any step. In the iterative 
approach, however, when new multiplets are generated 
through tensor products with smaller entities, numeri- 
cal error can easily pile up. In this case, the accuracy 
of the matrix elements of the generators must be better 
than double precision. For practical purposes, quad pre- 
cision on matrix elements of the generators turned out 
sufficient. 

Finally, for example, with standard CGC spaces of 



rank-3, the sparse storage of CGC spaces also requires 
the generalization of standard sparse storage and sparse 
operations to arbitrary rank tensors. All of these is 
achieved, in general, by proper efficient permutations in 
sparse index space which requires a fast sorting scheme, 
together with reshaping of higher rank-objects to stan- 
dard two-dimensional sparse objects with fused indices, 
since this allows to employ fast routines for standard 
sparse multiplication. 



